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FOREWORD 


A  tMo-day  seminar  on  Water  Quality  Data  Interpretation  was  held  In 
Atlanta,  Georgia  on  8-9  February  1978.  Xrhe  purpose  of seminar  was 
to  provide  a  forum  for  Corps  of  Engineers  personnel  who  are  routinely 
Involved  In  water  quality  data  collection,  analysis,  and  Interpretation. 
Topics  addressed  during  the  serntnar^ Include  data  collection,  data 
processing,  laboratory  quality  control,  and  Interpretation  of  physical- 
chemical,  biological,  and  chemical  contaminant  data.  Ten  papers  pre- 

,  I 

sen  ted  during  the  seminar  are  contalned-bere^nr- 

The  views  and  conclusions  expressed  In  these  proceedings  are  those 
of  the  authors  and  are  not  Intended  to  modify  or  replace  official 
guidance  or  directives  such  as  engineer  regulations,  manuals,  circulars, 
or  technical  letters  Issued  by  iJtxe  Office  of  the  Chief  of  Engineers. 


R.  G.  WILLEY 
Edi tor 
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THE  COMMITTEE  ON  WATER  QUALITY 


By 


Janice  E.  Rasgus 


1 


The  Committee  on  Water  Quality  was  formally  established  by  ER  15-2-10, 
dated  15  December  1975.  The  members  of  the  Committee  are  appointed  by 
the  Chief,  Engineering  Division,  Directorate  of  Civil  Works,  Office, 

Chief  of  Engineers  and  represent  each  Division  and  several  of  the  laboratories. 
The  Committee  Is  charged  with  providing  guidance  for  developing  a 
comprehensive,  coordinated  Corps-wide  water  quality  managment  program. 

The  Committee  met  three  times  during  1977  and  accomplished  the 
following: 

1)  Recommended  three  water  quality  training  courses  which  will  be 
held  during  FY  78.  These  are: 

a.  Selective  Withdrawal  Design. 

b.  Water  Quality  and  Ecosystem  Modeling. 

c.  Water  Quality  Aspects  of  Water  Control. 

2)  Provided  a  preliminary  listing  of  water  quality/environmental 
training  courses  available  to  Corps  personnel  through  other  agencies. 

3)  Developed  a  workshop  for  laboratory  chemists  at  WES  to  be  held 
on  22  and  23  February  1978.  The  purpose  of  the  workshop  is  to  discuss 
the  technical  aspects  of  laboratory  quality  control  and  analytical 
procedures . 

4)  Sponsored  a  water  quality  seminar,  "Water  Quality  Data  Collection 
and  Management",  which  was  held  In  Denver,  Colorado,  on  25-26  January 
1977.  A  proceedings  was  published  and  distributed  to  each  participant. 

5)  Revised  and  Issued  ER  1130-2-334,  dated  16  December  1977.  It 
requires  the  reporting  of  water  quality  activities. 

6)  Drafted  an  ER  which  addresses  quality  assurance  procedures  for 
Corps  and  contracted  laboratories  and  laboratory  personnel  needs. 

7)  Issued  ER  15-2-11,  dated  1  November  1977.  It  describes  the 
Committee's  consulting  service  and  the  procedures  for  obtaining  these 
services.  This  service  Is  one  of  the  Committee's  Important  functions 
and  Is  desirable  for  disseminating  Information  and  providing  a  means  for 
effective  Corps-wide  communication. 


North  Central  Division;  Vice  Chairman,  Committee  on  Water  Quality 

Paper  1 


8)  Consulted  with  the  Louisville  District  on  the  Cave  Run  reservoir 
project  In  Morehead,  Kentucky.  During  the  summer  months  the  reservoir 
experiences  low  dissolved  oxygen  concentrations  and  consequent  releases 
of  slugs  of  high  concentrations  of  Iron,  maganese  and  hydrogen  sulfide. 
The  meeting  helped  generate  a  better  understanding  of  the  problems  and 
developed  some  practical  solutions. 

During  the  coming  year  the  Committee  will  address  the  following 
Items: 

1)  Manpower  requirements  for  properly  staffing  a  water  quality 
management  function. 

2)  Data  management  systems. 

3)  Water  quality  laboratory  quality  control  and  personnel  needs. 

4)  Technology  transfer. 

5)  Training  needs. 

Members  of  the  Committee  should  be  contacted  with  any  questions, 
concerns  or  suggestions.  These  will  be  brought  to  the  attention  of  the 
Committee  for  discussion.  The  members  of  the  Committee  on  Water  Quality 
are: 

Dr.  Mark  Anthony  -  ORD 
Mr.  Charles  Bradshaw  -  LMVD 
Mr.  John  Bushman  -  OCE 
Mr.  Richard  DlBuono  -  NED 
Mr.  Harry  Dotson  -  SPD 
Mr.  Earl  Elker  -  OCE 
Dr.  Rex  Eley  -  WES 
Mr.  John  Grace  -  WES 
Mr.  Alfred  Harrison  -  MRD 
Mr.  David  Legg  -NPD 
Mr.  Thomas  Malsano  -  NAD 
Dr.  Harlan  McKlm  -  CRREL 
Mr.  Milton  Millard  -  OCE 
Ms.  Janice  Rasgus  -  NCD 
Mr.  Julian  Raynes  -  SAD 
Mr.  Charles  Sullivan  -  SWD 
Mr.  Jerry  Willey  -  HEC 
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ENVIRCMHENTAL  AND  WATER  QUALITY 
OPERATIONAL  STUDIES 
by 

Milton  Millard 


The  Civil  Works  Program  of  Che  Corps  is  unique  as  it  encompasses  the  en> 
tire  spectrum  of  Water  Resource  Development.  Traditionally,  the  role  has 
involved  the  planning,  design,  construction,  and  operations  of  water  re¬ 
sources,  primarily  for  navigation,  flood  control,  hydropower  and  associated 
activities.  The  recent  concern  for  the  preservation  and  protection  of  the 
environment  has  resulted  in  legislation  that  makes  environment  and  water 
quality  additional  considerations  in  the  management  of  these  water  resources 
and  has  added  new  dimensions  to  our  basic  management  requirements.  For  pro¬ 
jects  in  the  planning  phase,  this  means  that  the  environmental  and  water  quality 
objectives  must  be  addressed  along  with  the  basic  project  purposes.  For  pro¬ 
jects  Chat  are  operational,  environmental  and  water  quality  objectives  mean 
additional  pressures  to  meet  new  standards  and  reduce  adverse  environmental 
impacts.  To  help  meet  these  new  responsibilities,  the  Corps  has  undertaken 
another  major  research  effort. 

Without  going  into  detail  on  the  background  of  the  study,  the  sequence  of 
events  depicted  (Fig.  1)  followed  an  OCE  staff  study  made  in  1975.  Problem 
identification  and  assessment  was  accomplished  in  1976,  detailed  program 
planning  was  done  in  1977  and  the  research  and  study  program  is  being  ini¬ 
tiated  in  FY  1978.  For  budget  purposes,  the  program  was  Identified  as  "The 
Environmental-Water  Quality  Operational  Studies  Program:  and  was  Inmedlately 
referred  to  by  the  acronym,  "EWQOS." 

Environmental  quality  problem  Identification  and  assessment  was  the  first 
step  in  formulation  of  the  EWQOS  program.  A  field  office  survey  in  which 
sll  divisions  were  Involved  was  conducted  by  a  team  from  OCE  and  WES.  The 
major  purpose  of  this  survey  was  to  obtain  information  on  Che  impacts  of  the 
new  requirements  on  the  various  functional  areas  of  the  Corps*  organiaatibn, 
namely,  the  Planning,  Engineering,  and  Construction-Operations  Divisions. 

Some  of  Che  impacts  of  Environmental  Quality  problems  on  Civil  Works  activi¬ 
ties  are  shown  in  Figure  2. 

Environmental  and  Water  Quality  Problems  were  documented  at  152  Corps  Pro¬ 
jects  representing  200  megawatts  of  power  generating  capacity,  300  billion 
gallons  of  water  supply,  and  over  125  million  recreation  days  annually. 

A  major  finding  was  Chat  Che  environmental  and  water  quality  problems  most 
significantly  Impact  on  Project  Operations  and  are  not  limited  to  any  one  part 
of  Che  country. 


Milton  Millard,  Chief  Western  Section,  Operations  Branch,  Con-Ops  Division, 
Directorate  of  Civil  Works,  Office  Chief  of  Engineers 
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Figure  3  give*  eome  id««  of  how  the  Corps'  projects  sre  distributed  sround 
the  country*  but  does  not  necesssrily  indicste  the  distribution  of  the 
■sjor  environaentsl  probleu. 

The  end  product  of  the  first  phsse  of  this  prograi  wss  s  report  thst 
served  ss  Justificstion  for  the  progrsa  to  OMB  end  the  Congress.  A  few 
of  the  p rob lens  identified  end  questions  raised  by  the  report  sre: 

s.  vniat  effect  do  Corps  projects  have  on  water  quality*  is 

the  effect  significant*  and  is  there  a  need  to  find  a  solu¬ 
tion? 

b.  Under  the  present  energy  situation  and  the  increased  ia- 
portsnce  placed  on  hydropower*  what  are  the  effects  of 
hydropower  operations  at  punpback  facilities? 

c.  What  do  we  know  of  the  quality  of  reservoir  releases  and 
their  Impact  downstream? 

d.  What  are  the  impacts  of  flow  variations  from  reservoir 
projects?  What  are  the  minirniSB  flew  requirements?  What 
are  the  impacts  of  hydropower  peaking  operations  and  other 
flow  variations? 

e.  What  are  the  Impacts  of  fluctuating  reservoir  pool  ele¬ 
vations*  i.e.*  exposed  mud  flats*  aesthetics*  the  effects 
on  fisheries  and  other  habitats? 

f.  What  about  reservoir  eutrophication  and  problems  of  algal 
blooms?  Uow  does  this  affect  project  purposes  such  as 
water  supply  or  recreation? 

The  report  also  indicates  a  need  to  improve  procedures  in  the  preparation 
of  environmental  assessments*  a  need  for  definition  of  data  requirements 
and  Improved  data  management  techniques.  In  addition*  the  report  points 
out  the  need  for  a  better  understanding  of  the  impacts  of  the  various 
activities  in  and  on  the  waterways  on  the  environment.  What*  if  any* 
adverse  impacts  are  there  from  navigation*  from  the  construction  of  cut¬ 
offs  and  training  structures*  bank  erosion  control  methods*  and  the  nunt- 
erous  other  activities? 

The  formulation  of  the  EWQOS  Program  considered  the  problems*  the  possible 
resources*  the  time  to  be  available  for  the  work*  and  the  following 
concepts : 


A.  The  probability  of  Success 

B.  Technical  constraints 

C.  Cost  of  the  Research 

D.  Potential  solutions 
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The  program,  as  developed.  Is  organized  Into  six  research  projects  plus 
associated  field  studies  aimed  at  meeting  the  program  objectives.  These 
field  studies  form  a  major  part  of  the  program.  Examples  of  major  research 
areas  to  be  addressed  are: 

a.  The  Improved  understanding  of  reservoir  hydrodynamics 
through  the  use  of  scaled  models  to  get  a  better  feeling  for 
reservoir  project  effects  on  water  quality  and  ecology. 

b.  Control  of  nuisance  algal  blooms  -  what  are  the-major  factors 
causing  them  and  what  steps  may  be  taken  to  control  them? 

c.  Reservoir  Release  Criteria  -  what  Is  required  to  maintain  a 
suitable  downstream  habitat  how  may  the  effects  of  algal  re¬ 
leases  be  reduced? 

The  success  of  this  program  Is  dependent  on  the  Interactions  of  all  elements 
within  the  Corps  and  the  many  outside  Interests  Involved  In  environmental 
and  water  quality  problems. 

In  the  area  of  technology  transfer  and  application,  we  will  try  to  be 
Innovative  and  expand  from  the  usual  scientific  reports  that  usually  do 
not  get  to  the  field  operating  people.  In  each  report,  when  practicable, 
we  will  Include  a  "How  to  Apply"  section  Or  appendix  which  will  also 
be  published  as  a  technical  letter  or  manual.  In  cases  where  extensive 
research  work  units  are  Involved,  we  will  develop  synthesized  reports  and 
"How  To"  manuals.  In  addition,  where  practical,  we  will  provide  field 
deoionstrdtions  of  new  techniques  and  procedures. 

To  assist  In  this,  each  District  and  Division  will  be  expected  to  provide 
continuous  feedback  on  Its  research  need  and  the  results  of  the  program's 
outputs.  To  facilitate  this  feedback  and  assure  that  your  needs  are  fed 
Into  the  program,  a  Field  Review  Group  has  been  established  to  work  with 
the  technical  monitors  In  management  of  the  program.  You  should  all  become 
acquainted  with  your  division  representative  on  the  Field  Review  Group 
and  work  closely  with  him  to  assure  that  your  needs  are  considered. 

To  get  the  most  out  of  this  program  and  avoid  duplication  of  efforts, 
all  activities  have  been  and  will  continue  to  be  coordinated  closely  with 
Federal  and  State  agencies,  consultants  from  Industry  and  academia,  and 
with  private  Interest  groups.  Early  In  the  program  we  held  an  Interagency 
briefing  on  our  proposed  activities  and  received  valuable  feedback.  We 
propose  another  such  briefing  In  the  near  future  and  then  on  a  regular 
schedule  as  the  program  progresses. 

The  EWQOS  Program  Is  based  on  a  six-year  effort  at  an  anticipated  cost 
of  $30  million.  This  Is,  of  course,  a  planning  figure  which  will  be 
reviewed  annually  by  the  Administration  and  the  Congress. 

The  Initial  planning  for  the  program  has  been  completed  -  we  are  funded 
for  our  first  year's  efforts.  Field  sites  have  been  selected  and 
work  will  be  Initiated  after  the  selections  are  formally  approved  by 
the  Division  Engineers  Involved.  In  a  short  time,  the  first  RFP  will 
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Over  the  last  ten  years  there  has  been  a  dramatic  increase  In 
the  Interest  In  water  quality  problems  and  In  finding  solutions  to  these 
problems.  At  the  same  time  traditional  Corps  roles  In  water  resources 
development  and  management  have  magnified  these  concerns  relative  to 
our  activities.  Water  quality  studies  are  now  conducted  at  every  phase  of 
project  development*  survey  reports,  general  and  feature  design 
memoranda,  and  for  projects  that  are  completed  and  those  under  constrqctlon. 
Each  water  quality  study  Is  different*  yet  at  the  same  time  there  Is  a 
coomon  denominator,  and  the  common  denominator  Is  data.  You  simply 
can't  do  a  water  quality  study  without  data.  At  the  Denver  Seminar  last 
yearicwo  days  were  spent  talking  about  data  collection,  and  little  more 
was  done  than  scratching  tiie  surface. 

Someone  might  ask  the  question,  why  should  the  Corps  collect  water 
quality  data?  The  states  are  collecting  data,  ERA  Is  collecting  data, 
uses,  USFWS,  and  universities  are  all  collecting  data  for  one  purpose  or 
another.  But  that's  just  the  point,  their  purposes  for  collecting  data 
are  considerably  different  In  the  vast  majority  of  cases  from  ours. 

For  example,  when  EIA  collects  data  to  monitor  Industrial  discharges, 

99%  of  the  time  that  data  will  not  be  of  any  worth  to  a  person  sitting 
in  a  District  office  trying  to  regulate  a  reservoir.  This  leaves  us  In  a 
position  where  for  Corps  unique  purposes,  at  best  we  must  collect  data  to 
supplement  that  available  from  others  and  at  worst  we  must  conduct 
complete  data  collection  programs  of  our  own. 

While  looking  back  over  approximately  10  years  of  notes  related  to 
data  collection,  I  came  across  an  Interesting  piece  from  a  1969  training 
course  that  was  held  at  HEC.  I  think  It  points  out  some  problems  that 
still  exist  with  data  collection.  It  Is  not  always  possible  to  anticipate 
future  needs,  hence  there  exists  a  tendency  for  the  creation  of  data 
programs  for  data  sake  alone.  Ignoring  the  purposes  to  which  the  Information 
gathered  should  be  directed.  Regretably,  there  too  often  seems  to  be 
money  to  gather  more  data,  but  relatively  little  to  evaluate  It.  Thus, 
the  process  goes  on  at  ever  expanding  levels  until  we  become  virtually 
Inundated  with  data.  At  this  juncture,  we  may  have  to  give  some  attention 
to  data  disposal.  At  the  very  least  we  should  look  toward  optimizing 
our  efforts  In  data  management  so  that  we  get  maximum  value  for  each 
dollar  expended.  While  we  have  made  Improvements,  too  often  we  wind  up 
In  the  same  situation  today  as  we  were  In  10  years  ago.  We  collect  too 
much  data  and  often  not  the  right  data.  We  need  to  get  only  the  Information 
we  need  and  then  move  on  to  Interpret  It  and  come  up  with  the  solutions 
to  the  problems. 

1 

Hydrologic  Engineering  Section,  Office  Chief  of  Engineers 
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The  challenge  of  developing  and  managing  data  collection  programs 
Is  great.  It  taxes  the  Ingenuity  of  every  water  quality  manager.  A 
major  difficulty  In  developing  efficient  data  collection  efforts  Is  that 
"cookbook"  approaches  to  sampling  program  design  simply  can't  be  set  up. 

A  quick  look  at  why  we  collect  water  quality  data  will  Illustrate  this 
point.  Data  presently  collected  to  provide  project  baseline  Infoimatlon, 

for  early  Identification  of  problems,  to  provide  guidance  to  reservoir 
regulation  elements,  as  a  basis  for  development  of  design  criteria,  as 
input  to  ElS's  and  for  numerous  other  reasons.  The  needs  are  clearly 
different  and  the  programs  to  meet  these  needs  must  be  just  as  diverse. 

Under  all  this> however,  there  Is  a  logical  general  approach  to  sampling 
program  design  and  management  which* if  followed* will  produce  good  results. 

For  a  water  quality  data  collection  program  to  be  sucdessful,  many 
Important  factors  must  be  brought  together.  The  program  must  first 
be  established  to  meet  an  Identified  need.  It  must  address  only  very  well 
defined  objectives  which  have  been  developed  based  on  the  need.  It  has 
to  be  a  team  effort  In  the  sense  that  It  must  include  Input  from  both 
field  and  laboratory  personnel  In  order  to  reflect  reasonable  workloads  and 
minimize  the  possibility  of  technical  problems  arising  during  the  implementa¬ 
tion  phase.  Finally,  the  output  from  the  effort  must  satisfy  the  need 
that  the  program  was  designed  to  address. 

Figure  1  Illustrates  the  interactions  required  to  set  up  a  data 
collection  program.  The  data  collection  program  is  in  the  center.  We 
have  needs  considered,  a  program  nianager,  and  input  from  field  and  lab 
personnel.  Also,  output,  feedback  and  constraints  are  shown. 

Now  let's  take  these  Items  one  at  a  time.  There  are  probably  an 
Infinite  number  of  needs  for  which  to  collect  water  quality  data.  Some  are  to 
develop  baseline  Information,  to  assess  our  projects,  as  a  basis  for 
project  modification,  for  early  problem  Identification,  for  operations 
and  to  provide  reservoir  regulation  guidance.  We  collect  data  for 
criteria  development,  for  design,  for  EIS  preparation.  We  collect  data 
for  use  In  coordination  and  cooperation  with  other  agencies,  Federal 
and  state.  We  use  their  data  whenever  we  can  and  they  use  our  data 
whenever  they  can.  We  need  data  In  order  to  administer  the  regulatory 
program.  We  need  to  set  up  the  program  to  match  the  need. 

Next,  there  Is  the  program  manager.  What's  his  job?  He  first 
must  assimilate  the  needs  and  establish  the  objectives.  He  Is  responsible 
for  setting  up  and  conducting  the  program.  To  do  this  he  must  assess 
laboratory  and  field  capability  (with  suitable  Input  from  these  people) 

In  order  to  deteirmlne  whether  he  uses  Inhouse  people  or  contracts  the 
effort.  After  the  program  Is  underway  he  receives  the  data  and  Is 
responsible  for  data  management.  He  adapts  the  program  to  meet  changing 
objectives  and  needs.  He's  responsible  for  Interpreting  the  data.  He 
Is  responsible  for  the  output  which  Is  a  continuous  process.  He  has 
to  encourage  feedback.  This  Is  the  only  way  that  we  are  going  to  be 
able  to  keep  up  the  Interest  of  the  laboratory  and  field  people.  He 
must  let  them  know  what  the  data  Is  used  for  and  why  It  Is  needed. 
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what  are  the  constraints?  The  most  critical  one  is  manpower.  There 
never  seems  to  be  enough  people  to  do  all  the  work.  Another  is  time. 
Invariably  the  need  for  the  data  and  the  need  for  the  evaluation  of  the 
data  was  "yesterday"  so  we  are  faced  with  a  situation  where  we  must  set 
up  a  program  that  is  greatly  compressed  in  time.  The  final  constraint 
is  money.  But  keep  in  mind  we  need  a  good  balance  between  data  collection 
and  evaluation.  Both  are  needed  to  make  a  study  successful.  Data 
collected  but  not  evaluated  should  not  have  been  collected  in  the  first 
place. 


What  about  lab  personnel?  One  of  the  most  important  considerations 
here  is  the  need  for  laboratory  people  to  input  certain  technical  con¬ 
siderations  to  the  program  manager  so  he  can  be  fully  aware  of  the  problems 
the  laboratory  may  encounter  over  the  life  of  the  program.  Input  such 
as  required  level  of  detection  for  given  parameters  and  the  capability 
of  the  lab  as  it  relates  to  volume,  workload  and  what  analyses  the 
lab  is  capable  of  performing.  There  is  no  point  in  setting  up  a  lab 
program  if  your  laboratory  people  don't  have  the  capability  of  doing  it. 
Finally  the  lab  personnel  can  give  you  some  insight  as  to  "turnaround 
time"  you  can  reasonably  expect. 

What  about  field  personnel?  The  same  considerations  apply  to  field 
personnel.  If  you've  got  a  field  crew  to  collect  the  data,  their  technical 
input  and  knowledge  of  the  area  should  be  a  guide  in  setting  up  a  program; 
particularly  when  it  comes  to  frequency  of  sampling,  location  of  sampling 
stations,  etc.  You  can't  get  away  from  the  fact  that  the  man  who  knows 
the  area  can  provide  you  with  a  lot  of  information.  This  aspect  is  too 
often  neglected  in  setting  up  a  program. 

Finally  let's  consider  the  output.  This  item  not  only  represents  the 
final  report  but  also  the  continual  interpretation  and  feedback  of  data. 

We  must  continually  evaluate  our  data  so  that  we  can  make  improvements  and 
adjustments  to  the  program.  Output  also  Includes  data  management  whether 
it's  a  file  drawer,  a  computer  system  or  whatever. 

Now  how  do  we  pull  all  this  together?  Figure  2  shows  five  basic 
steps  needed  to  establish  an  effective  data  collection  program.  The 
first  step,  and  by  far  the  most  critical,  is  to  define  the  scope  and 
objectives  to  meet  the  identified  needs.  After  the  objectives  are  set 
we  must  assemble  all  the  available  data.  This  includes  any  pertinent  data 
collected  by  others  which  will  often  give  some  indication  of  what  we 
need  to  collect  and  how  to  collect  it.  Infortbatlon  such  as  available 
water  quality  data,  land  use  information,  point  and  non-point  sources 
of  pollution  in  the  watershed  and  maps  of  the  area  is  pertinent.  All  this 
information  should  be  assembled  and  evaluated  before  the  plan  is  put 
together.  Thirdly,  a  preliminary  program  is  established.  This  is  the 
point  where  the  laboratory  and  the  field  personnel  should  be  brought 
into  the  process  to  make  it  a  coordinated  effort.  A  good  idea  is  to 
keep  the  preliminary  plan  unconstrained  for  the  first  cut.  Since  we 
never  have  enough  money,  manpower  and  time  to  do  all  we  would  like  to 
do,  developing  an  unconstrained  program  first  will  give  us  the  base  to 
evaluate  the  effect  of  reductions  when  the  constraints  are  applied. 


3 


Paper  3 


In  other  words,  what  we  are  losing  as  we  reduce  the  scope  of  the  program 
to  meet  funding,  personnel  and  time  limitations.  The  fourth  item  is  to 
review  the  preliminary  plan,  possibly  collecting  some  water  quality  data 
to  support  out  initial  judgements.  This  review  enables  us  to  then  develop 
the  final  program.  Ask  yourself  again.  Does  the  program  address  the 
objectives?  Is  the  program  going  to  meet  the  needs?  Is  it  going  to  produce 
the  output  required?  Is  the  program  cost  effective?  If  the  answer  to 
these  questions  is  yes,  we  now  have  our  final  program. 

Now  for  the  data  collection  program  itself.  Figure  3  shows  seven 
items  common  to  any  water  quality  data  collection  program.  The  first 
one  is  the  objectives,  Ifeike  them  specific  -  define  them  well  and  you'll 
save  yourself  a  lot  of  problems  as  you  get  into  the  program.  The  second 
thing  is  parameters.  What  do  we  really  need  to  know?  It  may  be  all 
right  in  the  first  cut  to  take  a  large  number  of  samples  and  analyze  many 
parameters.  However,  you  better  come  back  to  the  office,  do  some 
interpretation  and  think  about  cutting  down.  Don't  be  drawn  into  the 
situation  where  you  are  collecting  a  lot  of  information  because  someone 
five  years  from  now  may  have  a  use  for  it.  The  third  item  is  sampling 
sites.  These  have  to  be  representative.  There's  no  use  sampling  unless 
it  gives  you  a  good  picture  of  what  the  water  quality  is.  The  fourth 
item  is  sampling  frequency.  Get  enough  samples  to  answer  the  questions 
but  not  too  many.  This  item  has  the  biggest  potential  for  wasted  effort. 
Often  data  is  collected  more  frequently  than  we  really  need  to  answer 
the  questions.  Next  is  sample  handling  and  preservation.  What  are  the 
requirements?  This  is  where  the  experience  of  the  laboratory  and  field 
people  can  be  of  great  assistance.  The  sixth  item  is  the  lab  and  field 
analyses.  What's  needed?  Should  the  effort  be  contracted  or  done  in 
house.  This  should  be  clearly  laid  out.  Along  these  same  lines,  resist 
the  temptation  to  put  water  quality  monitors  here,  there  and  everywhere. 
They're  often  expensive  to  run,  have  operation  and  maintenance  problems 
and  yield  more  data  than  we  need.  The  final  item  is  flexibility.  We 
should  be  constantly  reviewing  the  program  and  making  changes.  If  the 
program  can't  be  modified  with  relative  ease  our  task  will  be  much  more 
difficult. 

Factors  to  be  considered  under  each  of  the  seven  items  discussed 
above  are  shown  on  Figures  4-10.  While  the  figures  may  not  be  all 
inclusive,  they  should  nevertheless  provide  a  guide  to  setting  up  most 
water  quality  sampling  programs. 

In  summary,  developing  and  managing  water  quality  data  collection 
programs  is  not  an  easy  task.  It  requires  co-ordination  and  communication 
with  many  interested  parties.  Furthermore,  there  is  no  standard  procedure 
to  follow.  However,  if  the  ideas  presented  in  this  paper  are  followed, 
we  should  be  able  to  raeasureably  improve  our  efforts  in  this  important 
area. 


4 


Paper  3 


References 


1.  Committee  on  Water  Quality,  "Water  Quality  Data  Collection  and 
Management,"  Seminar  Proceedings,  25-26  January  1977,  Denver,  Colorado, 
U.S.  Army  Corps  of  Engineers,  Washington,  D.C.,  1977. 

2.  Training  Course,  Number  27,  Water  Quality  Management,  HEC,  Davis, 
California,  March  1969. 


Paper  3 


DATA  COLLECTION 
PROGRAM 
FLOW  CHART 


EXTERNAL 

constraints 


DATA  COLLECTION 


PR06RAf4  DEVELOPMENT 

1.  ASSIMILATE  NEEDS  &  ESTABLISH  SCOPE  AND  OBJECTIVES 

2.  ASSEMBLE  ALL  AVAILABLE  INFORMATION 

3.  ESTABLISH  PRELIMINARY  PLAN 

4.  REVIEW  PLAN 

5.  ESTABLISH  FINAL  PLAN 


Paper  3  figure  2 


DATA  COLLECTION 

pmm 


1.  OBJEQIVES 

2.  PARAMETERS 

3.  SAIWNG  SITES 

4.  SA^PLING  FREQUENCY 

5.  SA^PLE  HANDLING  &  PRESERVATION 

6.  LAB  &  FIELD  ANALYSES 

7.  FLEXIBILITY 


Paper  3  figure  3 


FAaORS  TO  CONSIDER 
IN  SETTING  OBJEaiVES 


1.  NEED  FOR  DATA 

2.  PHASE  OF  PROJEQ 

3.  TYPE  OF  PROJECT 

4.  CONSTRAINTS 

5.  PROBABILITY  OF  SUCCESS 


Paper  3  figure  4 


FAQORS  TO  CONSIDER 
IN  PAR/VtTER  SELEaiON 


1.  REUTIONSHIP  TO  OBJECTIVES 

2.  VOTIS®  LM)  USES 

3.  PROJECT  PURPOSES 

4.  KNOWN  OR /WICIPATED  PROBLEMS 

5.  COST  OF  ANALYSES 

6.  SAll^  PRESERVATION  REQUIREMENTS 

7.  ty^ORATORY  SUPPORT  CAPABILITY 

8.  SECONDARY  DATA  REQUIREMENTS 
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FACTORS  TO  CONSIDER  IN 
SATfLING  SITE  SELECTION 


1.  RELATIONSHIP  TO  OBJECTIVES 

2.  Flow.  CURRENT  AND  MIXING  PATTERNS 

3.  PHVSIOV.  CHAIWnERISTICS  (F  WATERBtW 
A.  POLUJTANT  SOURCES 

5.  ACCESSABILITY  OF  SITE 

6.  ECONOMICS 


I 
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FACTORS  TO  CONSIDER 
IN  Sfm.m  FREQUENCY 


1.  RELATIONSHIP  TO  (UECTIVES 

2.  NORMAL  VARIATIONS  IN  CONCENTRATION  OF  PAR/VEIERS 

3.  HYDimiGIC  I^FL1JENCES 
ECONOMICS 


Paper  3  figure  7 


1 


F/CrORS  TO  CONSIDER 

IN  SmE  HANDLING,  PRESERVATION  AND  TRANSPORTATION 


1.  RELATIONSHIP  TO  OBJECTIVES 

2.  SATflE  CONTAINERS 

3.  Pf^SERVATION  AND  FILTRATION  REQUIRErerS 
A.  SAMPlf  STORAGE 

5.  TRANSPORTATION  AVAILABLE 
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FACTORS  TO  CONSIDER 
IN  lAB  AND  FIELD  ANALYSES 


1.  RELATIONSHIP  TO  CEJECTIVES 

2.  CAPABILITIES  -  PAR/WERS  AND  VOUU^ 

3.  ACCURACY  REQUIREreirS  OF  STUDY 
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OHIO  RIVER  DIVISION'S 
LABORATORY  MASTER  FILE  SYSTEM 

1 

Henry  C.  Jackson 

1.  General.  Before  the  creation  of  the  LABMASTER  system,  the  division 
water  quality  laboratory  received  san^les  from  the  field,  tested  them, 
and  prepared  reports  of  test  results  manually.  LABMASTER  trats  created  as 
a  means  of  systematically  recording  stemdard  information  about  sanqples 
entering  the  laboratory,  preparing  bench  sheets  describing  all  sauries 
logged  into  the  system,  recording  test  results,  billing  districts  for 
testing,  preparing  reports  of  completed  and  pending  results,  and  pro¬ 
viding  sumagement  tools  for  efficient  operation  of  the  laboratory. 

Field  results  were  stored  with  laboratory  results  for  a  given  sample  to 
insure  a  complete  information  base.  While  the  system  was  originally 
designed  to  run  on  the  G437  conputer,  the  G437  computer  can  no  longer 
provide  adequate  service  (econcxny,  turnaround  time,  special  plotting 
packages,  time  sharing,  etc.).  Therefore  the  information  retrieval 
portion,  as  well  as  the  field  data  storage  portion  of  the  system,  has 
been  moved  to  INFOMET  (Computer  Sciences  Corporation  INFONET)  aund  only 
the  division  laboratory  portions  will  remain  on  the  G437.  Programs  on 
the  G437  should  be  run  only  by  the  division  laboratory. 

2.  gyrating  Environment.  The  G437  portion  of  the  system  is  designed 
to  run  on  a  Corps  standard  64K  G437  computer  with  9  disc  drives,  6 
7-track  tape  drives,  card  punch,  card  reader,  emd  line  printer.  The 
INFONET  por~tion  of  the  system  is  designed  to  mn  on  a  Computer  Sciences 
Corporation  UNIVAC  1108  con^uter  with  associated  peripherals  including 
disc  drives  and  9-track  tape  drives.  Access  to  the  INFONET  system  is 
made  thru  COPE-1200  terminals  with  card  reader,  c£urd  punch,  line  printer, 
and  7-track  tape  drives  as  well  as  thru  a  variety  of  low-speed  inter¬ 
active  terminals. 


3.  Programs .  The  LABMASTER  system  is  currently  supported  by  12 
programs  on  the  G437  for  division  laboratory  use  and  20  programs  on 
INFONET  for  storage  and  retrieval  by  districts.  Other  programs  are 
available  but  not  supported  on  the  G437.  Detailed  descriptions  of  these 
programs  may  be  found  in  Appendix  B.  The  programs  are  as  follows: 


G437 


a. 

401H0130 

Laboratory 

Master 

File  Addition 

b. 

401H0140 

Laboratory 

Master 

File  Update 

c. 

401H0150 

Laboratory 

Master 

File  Billing 

d. 

401H0170 

Laboratory 

Master 

File  Station  List 

e. 

401H0190 

Laboratory 

Master 

File  Bench  Sheet 

f. 

401H0200 

Laboratory 

Master 

File  Station  Correction 

9- 

401H0210 

Laboratory 

Master 

File  Table 

Ohio  River  Division,  Reservoir  Control  Center,  Hater  Quality  Section 
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401H0220 

Laboratory  Master 

Pile 

Sample  Age 

i. 

401H0225 

Laboratory  Master 

File 

Management  Information 

j- 

401H0310 

Laboratory  Master 
Correction 

File 

Station- Identification 

k. 

401H0430 

Laboratory  Master 

File 

INFONET  Transfer 

1. 

401H0440 

Laboratory  Master 

File 

Interim  Bill 

INFONET 

a. 

401265 

Laboratory  Master 

File 

Sorted  General  Report 

b. 

401340 

Laboratory  Master 

File 

Field  Data  Input 

c. 

401CAT 

Laboratory  Master 

File 

Catalog 

d. 

401EDT.JA 

Laboratory  Master 

File 

Pre-edit 

e. 

401GET 

L2j3oratory  Master  File  G437  to  INFONET  Transfer 
(Translate  aind  Separate) 

f. 

401GEX 

Laboratory  Master 
Transfer 

File 

G437  to  INFONET  Correction 

g- 

401IND 

Laboratory  Master 
River  Index 

File 

Station  Selection  by 

h. 

401LOD 

Laboratory  Master 

File 

System  2000  Data  Base  Load 

i. 

401MER 

Laboratory  Master 

File 

Update  emd  Merge 

j- 

4010BY 

Laboratory  Master 

File 

Define  and  Equate  Files 

k. 

401RET 

Ledxiratory  Master 

File 

Data  Retrieval 

401RIV 

Laboratory  Master 

File 

River  Data  Report 

m. 

401SE1 

Laboratory  Master 

File 

Data  Selection  by  Type 

n. 

401SR1, 

L^Jx)ratory  Master 

File 

Data  Selection 

o. 

401STA 

Laboratory  Master 

File 

Station  Statistical  Summary 

P- 

401STD 

Ledxiratory  Master 

File 

Statistics 

q- 

401SUB 

Laboratory  Master 

File 

Station- ID  Selection 
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SYSTEM  FLOWCHRRT 
FOR  G437  update 


Fig.  1 
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r.  401SUM  Laboratory  Master  File  Station-Day  Summary 

s.  401UNP  Laboratory  Master  File  Unpunctuated  File  Creation 

t.  40ixyz  Laboratory  Master  File  XYZ  Coordinate  Append 

4.  Interaction  of  Programs  on  the  G437.  When  a  water  quality  sample 

is  received  by  the  division  laboratory,  it  is  assigned  a  leiboratory 
number.  The  tag  attached  to  the  sample  contains  information  about  the 
sample  which  is  keypunched  directly  from  the  tag  (see  figure  1  for 
system  flowchart)  run  with  401H0130  to  log  the  sample  into  the  system. 
Laboratory  bench  sheets  are  prepared  by  running  401H0190.  Management 
information  programs  such  as  401H0220  or  401H0225  provide  guidance  to 
leiboratory  personnel  in  choosing  which  tests  to  perform.  After  tests 
are  performed  and  results  recorded  on  bench  sheets,  input  to  401H0140  is 
keypunched  directly  from  the  completed  bench  sheets.  Billing  for  labo¬ 
ratory  work  is  accomplished  by  401H0150  which  is  run  on  a  quarterly 
basis.  Interim  billing  information  is  obtained  by  running  401H0440 
during  each  update  cycle.  Information  to  be  transferred  to  INFOIIET  is 
obtained  from  the  LABMASTR  tape  (LABMASTR  is  the  8-character  tape  file 
name  used  on  the  G437)  by  401H0430.  A  complete  list  of  stations  may  be 
obtained  by  running  401H0170.  Various  programs  such  as  401H0200  or 
401H0310  may  be  used  to  correct  information  after  it  has  been  placed  on 
a  LABMASTR  tape  file. 

5.  Interaction  of  Programs  on  INFONET.  Leiboratory  results  on  tape 
extracted  by  401H0430  from  the  G437  are  transmitted  to  INFONET  via  the 
COPE-1200  terminal  and  placed  on  a  tape  file  (4  backup  tapes  are  main¬ 
tained)  .  Files  referenced  to  the  current  data  in  each  of  the  district 
water  quality  INFONET  libraries  are  defined  and  equated  by  4010BY  (see 
figure  2  for  system  flowchart) .  Necessary  character  translation  and 
separation  by  district  code  is  done  by  401GET.  Special  correction  files 
in  each  district  library  are  created  by  401GEX  on  an  as-needed  basis. 

Field  data  input  to  INFONET  may  be  accomplished  by  using  401340  to 
generate  LABMASTER  format  data  records  (or  by  any  suitable  user-coded 
program).  Disc  files  created  by  401GET,  401GEX,  401340,  etc.,  may  be 
edited  by  401EDT.JA  before  placing  the  data  on  a  tape  master  file  (LABFIL) 
with  401MER  (see  figure  3  for  system  flowchart) .  A  detailed  catalog  of 

a  LABFIL  may  be  obtained  by  running  401CAT.  A  disc-resident  subset  of 
the  data  on  a  LABFIL  may  be  obtained  by  running  401SE1,  401SEL,  401SUB, 
or  401RET.  A  formatted  report  of  information  on  a  LABFIL  or  a  subfile 
may  be  obtained  by  running  401265,  401RIV,  or  401SUM.  A  system  2000 
data  base  loader  input  may  be  obtained  by  running  401LOD.  Standard 
statistical  analysis  of  the  information  in  a  sorted  subfile  may  lie 
obtained  by  running  401STD  or  401STA.  Temporary  subfiles  (one  file  per 
project)  for  binding  into  a  quick- access  unpunctuated  tape  may  be  obtained 
by  running  401UNP.  Raw  data  records  may  be  modified  for  plotting  purposes 
by  running  401XYZ  which  appends  the  river  mile,  elevation,  and  distance 
from  left  bank  to  each  record.  Stations  may  be  selected  by  401IND  which 
allows  the  user  to  specify  STORET  streeun  indices  and  river  miles  as 
selection  criteria.  Subsequent  raw.  data  retrieval  may  be  made  by  401RET. 

A  detailed  description  of  each  station  is  contained  in  a  special  random 
access  file  called  DESFIL  which  may  be  accessed  by  any  program.  Examples 
of  program  use  may  be  found  in  Appendix  C. 
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SYSTEM  FLOWCHART 
FOR  DATA  TRANSFER 


/LflB-TflPE 


COPE  1200 


SYSTEM 

* 

/^DRT.HO  / 

L  U«  T 

CURRENT 

OflTUL.HO 


4010BY 


system 

COPY 


foot  .ho 


401GET 


/YtHfOO  HI 
/UZV3!!] 


jm”*! 


f  TYnrOB  HB  I 
/UYHiOS 


REPORT 


I  YYHrOB .H«  , 

/UZ^OBY  I 


401CflT 

401CHT 

401CHT 

401CBT 

HUNTINGTON 

CRTRLOC 

LOUIS^ULE 

crtrloc 

NASHVILLE 

catalog 

PITTSBURGH 
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BRD  R^C  ERRFIL  SHOULD  BE  CORRECTED  WITH  SYSTEM  EDITOR  BEFORE  NEXT  UPOfl^'E 
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6.  Data  Elements.  A  standard  method  of  describing  the  test,  location, 
and  time  of  sem^llng  was  required.  The  following  descriptive  terms  have 
been  adopted: 

TEST  Five  character  STORET  parameter  code  for  a  st2mdard 

laboratory  test  (see  Appendix  A  for  list) . 

LABORATORY  Five-digit  nvimber  assigned  to  each  sample  by  the  lab- 
NUMBER  oratory . 


STATION 


Nine-character  location  code  described  below: 


DISTRICT 

1 

2 

3 

4 

5-9 

PROJECT 

TYPE 

1 

2 

3 

4 

5 

6 

7 

8 

0  or  9 


One-digit  code 
Huntington  District 
Louisville  District 
Nashville  District 
Pittsburgh  District 
Other 

Three-character  code  assigned  by  district 

One  digit  type  of  sampling  code  assigned  by  district 

River 

Lake 

Benthic 

Sewage 

Industrial 

Potable  Water 

Subsurface  Water 

Surface  Soil 

Other 


RIVER  MILE 


Four-digit  code  assigned  by  district 


IDENTIFICA-  Thirteen-digit  time  and  depth  code  described  below. 
TION  All  fields  are  right- justified  and  zero-filled: 


YEAR 


Two-digit  calendar  year 
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MONTH 


Two-digit  month 


DAY 

TIME 

DEPTH 

CONCENTRA¬ 

TION 


Two-digit  day  of  month 
Four-digit  time  of  day  (24-hour  clock) 
Three-digit  depth  of  sampling 
Seven-character  test  result  described  below: 


TRACE  INDICATOR 


L  Minimal  test  result 

G  Meiximal  test  result 


)i  Normal 

VALUE  Six-character  numeric  test  result  (decimal  point, 

plus  sign,  emd  minus  sign  are  permitted) . 

COST  Five-digit  test  cost  in  cents 

BILL  INDICATOR  if  billed;  space  otherwise 


EPA  INDICATOR  Three -character  flag  described  below: 

EPA  Data  on  record  previously  traxvsferred  to  STORET 

REV  Test  result  on  record  has  t>een  changed  and  not 

trcuisferred  to  STORET 


ANY 


A)  normally  described  Station- ID  appears  on  this 
record 


Data  not  trcUisferred  to  STORET 

INFONET  One-character  indicator  (with  meaning  only  on  G437 

INDICATOR  system)  as  described  below: 


I 


Data  transferred  to  INFONET 


Data  not  transferred  to  INFONET 


TRANSFER  One-character  indicator  (with  meaning  only  on  G437 

INDICATOR  system)  as  described  below: 


T  Data  tremsferred  to  district  LABMASTER 

Data  not  transferred  to  district  LABMASTER  (obsolete) 
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A  special  coding  system  has  been  devised  to  facilitate  the  naming  of 
main  stem  Ohio  River  stations.  This  coding  system  Is  as  follows: 

POSITION  INFORMATION 

1  District  code 

2-3  OR 

4  Hundreds  position  of  river  mile  (i.e.,  6  for  river 

mile  625.7  from  Pittsburgh) 

5  TYPE  CODE  (usually  1) 

6-8  Tens,  units,  emd  tenths  positions  of  river  mile 

(i.e.,  257  for  river  mile  625.7) 

9  Relative  position  in  tenths  of  cross  section  (0 

for  left  bank  to  9  for  right  b2mk) 

7.  File  Names  and  Structure  of  the  G437  System.  Tape  and  disc  files 
on  the  G437  system  must  be  eight-character  names.  The  master  tape  file 
containing  all  laboratory  results  as  well  as  work-in-progress  is  called 
LABMASTR.  The  LABMASTR  tape  has  a  standard  GE  tape  ledsel  with  a  VALUE- 
OF-ID  of  "LABMASTR."  Records  are  64  characters  long  and  blocks  are  100 
records  long  with  block  serial  m’mbers.  The  first  record  of  the  tape 
consists  of  all  zeros  except  columns  33-40  which  contain  the  last 
billing  date.  The  remaining  data,  as  described  below,  are  sorted  in 
ascending  order  on  the  first  32  coltnnns: 


POSITION 

INFORMATION 

1-5 

TEST 

6-10 

LABORATORY  NUMBER 

11-19 

STATION 

20-32 

IDENTIFICATION 

33-39 

CONCENTRATION 

40-44 

COST 

45 

BILL  INDICATOR 

46-59 

Not  Used 

60-62 

EPA  INDICATOR 

63 

INFONET  INDICATOR 

64 

TRANSFER  INDICATOR 
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The  teible  of  acceptable  STORET  parameters  and  associated  information  is 
stored  in  a  permcuient  disc  file  called  "XXXXXXXX"  on  the  "IDSEDITl" 
disc-  Records  are  64  characters  long  emd  the  blocking  factor  is  60, 
Data  are  sorted  in  ascending  order  on  the  first  5  columns.  The  record 


description  is  as 

follows : 

COLUMN 

INFORMATION 

1-5 

STORET  Parameter  Code 

6-13 

Line  1  of  Title 

14-21 

Line  2  of  Title 

22-29 

Line  3  of  Title 

30-34 

Cost  in  Cents 

35-42 

Obsolete  Number  in  401H0160 

43-46 

Shelf  Life  Code 

47 

Scunple  Type  Code  1=WATER 

3=MUD,  otherwise  blank 

48 

Minimum  Value  Flag  (L  if 
space  otherwise) 

trace  values  are  reported; 

49-54 

Minimum  Acceptable  Value 

55 

Maximal  Value  Flag  (G  if 
space  otherwise) 

maucimal  values  are  reported 

56-61 

Maximum  Acceptable  Value 

The  code  which  is  used  in  the  tremsmission  of  information  from  the  G437 
to  INFONET  is  called  "L’B-TAPE."  It  is  an  unlabeled  card  image  tape 
with  80  character  records  unblocked.  The  first  64  characters  of  each 
record  are  the  same  as  on  the  LABMASTR  tape  and  the  remaining  16  char¬ 
acters  are  blcink. 

8.  File  Names  and  Structure  on  INFONET.  Tape  euid  disc  files  on 
INFONET  have  six-character  names  (may  be  less  than  6)  emd  an  optional 
two-character  version  which  may  be  appended  to  the  end  of  the  file  name 
(separated  by  a  period) .  Special  tape  files  defined  as  "generation  data 
sets"  and  providing  backup  cap2d3ility  ("grandfather,"  "father,"  son," 
etc.)  are  described  by  affixing  a  comma  and  a  backup  number  to  the  end 
of  a  tape  file  name  (i.e.,  MYTAPE.H1,2) .  File  names  may  be  further 
qualified  by  eq>pending  a  slash  and  the  library  name  (i.e.,  401MER/ORDLAB 
would  be  in  the  shared  division  library  emd  ERRFIL.H1/UZV003  would  be  in 
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the  Huntington  District  library).  Unless  otherwise  stated,  all  files 
have  standard  keys  and  records  are  64  characters  long.  The  standard 
sort  sequence  is  ascending  order  of  columns  11-32  and  1-5.  Standard 
version  names  for  LABMASTER  files  are  as  follows: 

.HO  for  Division  Office 

.Hi  for  Hiintington  District 

.H2  for  Louisville  District 

.H3  for  Nashville  District 

.H4  for  Pittsburgh  District 

Stemdard  file  names  (with  version  represented  by  xx)  are  as  follows: 


BAD. XX  File  containing  records  rejected  by  editor  program 

(401EDT. JA/ORDLAB)  to  be  corrected  with  the  system 
editor  and  re-edited. 


DAT. HO 


Temporary  disc  file  containing  laboratory  data 
most  recently  transferred  to  INFONET. 


DATFIL.HOfO  Tape  file  containing  laboratory  data  transferred 
to  INFONET  (,0  indicates  most  recent  data). 


DESFIL.xx  File  containing  station  description  information 

in  the  EPA  STORET  format  as  described  in  the  description 
for  program  401DES,  DESFIL  may  be  accessed  randomly 
or  sequentially  because  it  has  3-word  nonstandard 
keys  corresponding  to  the  station  code  and  STORET 
record  number. 


DUPFIL.xx  File  containing  duplicate  records  detected  by  master 

file  update  program  (401MER/ORDLAB)  to  lae  checked 
and  (in  most  cases)  deleted  with  the  system  editor. 

EDTINP.xx  File  containing  input  data  for  editor  program 

(401EDT. JA/ORDLAB) . 


ERRFIL.XX 


FLDINP.XX 


FLDOUT.XX 


LABCOP.XX 


File  containing  records  rejected  by  master  file 
update  program  (401MER/0RDLAB)  to  be  corrected  with 
the  system  editor  cUJd  re-edited. 

File  containing  field  data  input  to  be  run  with 
field  data  format  program  (401340/ORDLAB) . 

File  containing  field  data  after  reformating  into 
standard  format  by  401340/ORDLAB. 

Tape  resident  duplicate  of  lednratory  master  file 
to  be  used  as  laackup. 
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LABDEF  System  2000  data  base  definition  file. 

LABLOO  System  2000  data  base  load  instruction  file. 

LABFIL.xx  Tape  resident  input  laboratory  master  file  con¬ 

taining  all  district  leUboratory  results. 

LABFIL.xx, 3  Tape  resident  output  laboratory  master  file  (3 

indicates  that  3  generations  of  backup  tapes  are  to 
be  maintained)  which  is  created  by  401MER/ORDLAB . 

LABINP.xx  File  containing  laboratory  data  after  merging  all 

ledioratory  provided  files  (i.e.,  770105. xx,  77999. xx, 
etc. ) . 

LABUNP.xx-UNP  Unpunctuated  tape  resident  file  containing  one 

subfile  for  each  project  (input  files  are  generated 
by  401UNP/ORDLAB)  which  is  used  for  quick  access  of 
data  for  an  entire  project. 

LOAODO  System  3000  data  base  load  "do  file." 

MERINP.xx  File  containing  all  data  to  be  updated  on  master 

file  by  401MER/ORDLAB. 

STAFIL  File  containing  stations  selected  by  401IND  which 

will  later  be  used  as  input  to  401RET. 

STALEG  File  containing  acceptable  stations  for  use  by 

401EDT. JA/ORDLAB.  Station  number  is  in  positions  1-9, 
maximum  depth  is  in  positions  10-12,  and  EPA  transfer 
code  is  in  positions  13-15  ("EPA"  or  blank) .  This 
file  is  sorted  on  positions  1-9. 

TABFIL  File  containing  information  on  every  acceptable 

STORET  parameter  (semie  format  as  XXXXXXXX  on  IDSEDITI 
on  G437) . 

yyirandd.xx  File  containing  laboratory  information  placed  on 

INFONET  on  year  yy,  month  mm,  and  day  dd  (i.e., 

770109. XX  would  contain  information  for  9  Jan  77) . 

yy99zz.xx  File  containing  corrections  to  laboratory  information 

on  INFONET  (i.e.,  779994. xx)  which  should  be  carefully 
inspected  by  district  personnel  before  using  in  eui 
update . 

9.  Special  INFONET  System  Routines.  There  are  a  few  system  routines 
on  INFONET  which  are  especially  useful  in  massaging  amd  organizing  data 
files.  They  are  well  documented  in  the  INFONET  manuals.  The  following 
is  a  list  of  some  of  the  more  useful  available  routines: 
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a.  System  Editor.  This  is  a  powerful  text  editor  which  is 
invoked  by  the  coIBn^uld  lEDIT  file  name.  The  user  should  refer  to  the 
manuals  or  pocket  guide  for  a  full  list  of  editor  conoiands. 

b.  System  Sort.  This  is  an  easy  to  use  sorting  routine  which  may 
be  executed  by  a  single  comnand.  A  typical  sort  coimnemd  to  sort  a  file 
in  the  standard  LABMASTER  sequence  «»ould  be  ISORT  64  inputfile  outputfile 
11,22  1,5. 


c.  System  Merge.  This  program  provides  an  easy  way  to  concatenate 
two  or  more  files  into  another  file  without  disturbing  their  order,  or 

in  the  case  of  already  sorted  files,  it  may  be  used  to  create  a  single 
sorted  file  from  two  or  more  input  files.  A  typical  merge  command  used 
with  LABMASTER  data  files  is  !MERGE,COPy  64  filel,file2,file3  outputfile. 

d.  SYSTEM  2000.  System  2000  is  em  extremely  powerful  data 
management  language  which  can  enable  a  user  to  perform  a  limitless 
variety  of  queries  on  a  properly  created  datcdsase.  It  is  essential  that 
the  user  have  a  System  2000  manual  before  attempting  to  use  this  program. 

e.  ALADIN  cmd  REPORT  II.  ALADIN  is  a  more  economical  data  manage¬ 
ment  lamguage  which  provides  cm  interface  with  the  REPORT  II  report 
%n:iter  which  enables  the  user  to  prepare  one-time  especially  formatted 
reports  quickly  2md  easily. 
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INTRODUCTION 


Vteter  quality  data  represents  infomation  about  the  state  and  behavior 
of  the  chemical  and  physical  aspects  of  a  hydrologic  system.  Analyses  of  this 
information  are  directed  tcward  separating  the  consequential  fran  the 
incaansequential  aspects  of  the  data,  and  point  to  a  better  understanding  of  the 
underlying  causal  and  randan  mechanisms  of  the  system. 

Uncertainty  in  the  water  quality  data  may  arise  fran  a  multitude  of  sources, 
including  unknown  or  uncertain  system  inputs,  uncertain  system  responses  to 
known  inputs ,  and  obser/ational  errors  resulting  from  sampling,  storaoe  and 
testing  procedures.  In  spite  of  these  uncertainties,  kncsi/ledge  often  is  available 
vhich  cdlows  a  deterministic  appraisal  of  at  least  a  portion  of  the  response  of 
the  system  being  monitored,  so  that  the  system  response  may  be  partitioned  into 
a  deterministic  portion  and  a  residual.  In  this  discussion,  techniques  are 
presented  which  enable  a  rational  analysis  of  the  chserved  data  to  be  made  with 
a  minimum  a  priori  knowledge  of  the  behavior  of  the  system.  The  procedures  cu:e 
designed  to  aid  ini  the  identification  and  characterization  of  both  the  randan 
and  non-randan  ocnponents  of  the  data. 

CkK:e  these  conponents  of  system  respc»)se  have  been  identified,  causative 
inechaniams  for  the  observed  non-randcm  behavior  nay  be  sought,  and  the  non-randon 
portion  of  the  data  used  for  calibratiai  of  the  deterministic  models  thus 
developed.  The  statistical  properties  of  the  random  conponents  of  the  data  may 
be  used  in  conjuction  with  the  detenninistic  system  response  to  predict  average 
and  extreme  system  behavior. 

The  statistical  analysis  of  data  is  an  integral  part  of  the  development 
and  calibration  of  models  of  hydrologic  systems.  Further,  because  the  process 
described  above  allows  the  observed  system  behavior  to  be  separated  into 
detenninistic  and  "random"  coiponents,  such  analyses  generally  should  allow 
realistic  estimates  to  be  made  of  "normal"  and  "extreme"  conditions  vhich  may 
occur  in  the  future. 

FITTING  OF  DETSBMINISTIC  TO  SAMPLE  DATA 

In  order  to  provide  motivation  for  the  study  of  the  abstract  concepts  of 
i^pectred  analysis  of  data  to  be  presented  later  in  this  section,  we  begin  our 
discussion  of  data  analysis  by  examining  the  fitting  of  detenninistic  models  to 
observed  water  quality  data. 

^Nuclear  Engineering  Conpai^,  Inc.,  Icuisville,  Kentucky.  Formerly  Chief 
Hydroixjgist,  law  Eiigineering  Testing  CJorpany,  Marietta,  Geor^. 
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Such  analyses  have  infxsrtant  practical  applications.  For  exanple,  veter  quality 
models  are  often  calibrated  by  adjusting  the  model  parameters  until  the  model 
output  is  "close"  to  the  observed  monitoring  data.  *1116  adequacy  .of  aujreanait 
between  observed  data  and  model  output  may  be  detonnined  subjectiv^y  or 
objectively. 

Hie  most  familiar  objective  criterion  is  the  method  of  least  squares.  In 
this  method,  the  parameters  of  the  model  are  adiusted  until  the  sun  of  the  squares 
of  the  differences,  tetween  predicted  and  observed  values  is  a  miniinum.  In  F<^y 
cases,  the  differences,  or  residuals,  obtained  by  this  procedure  have  desireable 
statistical  properties,  and  thus  the  prediction  of  the  livelihood  of  occurrence 
of  euctreme  events  through  such  an  analysis  is  enhanced. 

As  an  example  of  the  techniques  of  parameter  estiiration  by  the  method  of 
least  squares,  suppose  we  generate  a  sequence  of  temperature  measurements  Over 
a  one  year  period  in  a  standard  stream.  We  assvme  t^t  the  (unknown)  function 
relating  tenperature  to  time  is: 

T  =  68  +  5  sin  {2wt)  +  e 

Miere  t  is  time  in  years,  and  e  is  a  normally  -  distrihuted  randan  variate  from 
N  (0,  16) .  Using  the  abcye  equaticxi  and  a  table  of  normal  randan  variates,  the 
following  "observaticais"  may  be  generated. 


t 

T  =  T  -  e 

e 

T 

0 

68 

-0.12 

67.88 

,25 

73 

2.36 

75.36 

.5 

68 

1.65 

69.65 

0.75 

63 

-3.44 

59.56 

1 

68 

-1.68 

66.32 

in  a  pratical  problem,  we  will  have  available  the  observations  T,  and  the 
parametric  model 

t'  =  a  +  b  sin  (2iTt) . 

Our  objective  is  to  use  the  observed  data  to  estinete  the  parameters  a  and  b 
as  well  as  the  population  parameters  of  the  distribution  of  the  randan  corponent 
of  the  observations.  The  error  associated  with  the  ith  observation  is  defineti 
by  the  relation 

®i  “  ^i  “  (2fft3 
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Using  the  usual  prcxredures  of  least  squares  analysis,  we  obtain  the  parameter 
estimates 

k  =  67.54 

and  A 

b  =  7.9. 

An  analysis  of  the  residuals  yields  the  following  estimates  of  the  sample 
mean  and  variances 

W=  0 

and  2 

S  ^  =  5.84. 
e 

Using  our  derived  relation,  we  may  estimate  the  tenperature  at  t  =  0:5  vAiich, 
on  the  average,  will  be  exceeded  only  five  percent  of  the  time.  To  do  this, 
we  first  oonpute  the  deterministic  portion  of  the  system  response 

I  A  A  . 

from  the  relation  T  =  a  +  b  sin  (2irt)  =  67.754 

we  tlien  ccrpute  Oq 

®o..95  =  ^  ^^0.95^  “  3.981 

vtfiere  Z  is  the  .95  percentile  point  of  the  standard  normal  distribution. 

To.95  =  67.54  +  3.981  =  71.735 

Since  in  this  exanple  we  knew  the  model  frem  which  the  observations  were  generated , 
we  may  carpute  the  actual,  or  true,  value  of  T„  95  as 

T  =  68  +  4  (1.645)  =  74.58 
0.95 

We  see  by  conparison  of  the  above  estimate  of  with  the  true  value  that  our 

best  estimate  of  gg  is  too  low  by  almost  3®.  For  some  applications,  the  risk 
of  making  such  an  error  may  not  be  acceptable,  and  in  these  instances  we  are  led 
to  the  consideration  of  tolerance  limits  of  our  estimate.  By  a  tolerance  limit 
we  mean  an  estimate  of  (for  exanple)  gg  for  which  vre  can  be  assured  that  if 
we  make  such  an  estimate  in  a  large  ncn«5er  of  cases,  the  inequality 

^0.95  (estimated) 2 'f 0.95  (true) 

will  be  true  at  least  a  specified  number  of  times.  To  see  how  such  an  estinate 
may  be  constructed,  we  first  address  the  problem  of  describing  confidence  limits 
on  our  estimates  of  the  parameters  a  and  b.  A  confidence  limit  in  this  case 
means  a  region  in  tl.e  a  -  b  plane  in  which,  if  we  construct  such  regions  for  a 
large  number  of  cases,  the  true  values  of  a  and  b  will  lie  within  at  least  a 
specified  percentage  of  the  regions.  Such  a  region  may  be  constructed  as  follows 
(Grant,  1973) ;  Define  the  regression  conponent  of  the  errors  as 

N 

Reg  (e)  =  r  |a  +  b  sin  (2nti)  -  'y 


Paper  5 
3 


and  the  residue  of  the  errors  as 

2 

Res  (e)  =  E  e^  -  Reg  (e) 

The  function  Reg  (e)  has  rank  2  and  Res  (e)  has  rank  N-2  degrees  of  freedom. 

The  a  percent  confidence  region  for  a  and  b  is  then  defined  as  the  region  of 
all  pairs  of  points  (a,b)  which  satisfy  the  inequality 

Reg(e)  ^2  F  (  a  ;  2,  N--2) 

Reg (e)  N-2  100 

pirnire  1  shews  a  sketch  of  the  95%  confidence  region  for  a  and  b  in  our  exairple. 

Once  we  lave  available  the  95%  confidence  region  for  the  parameters 
a  and  b,  we  nay  define  the  95%  tolerance  limit  of  the  estimate  t'  as  the 
maxinum  of 


T  =  a  +  b  sin  (2iit) 

over  the  95%  confidence  region.  The  population  parameters  of  the  residual 
errors  (in  this  case,  the  standard  deviation) ,  are  estimated  from  the  errors 
evciluated  at  the  critical  Vcilues  of  a  and  b.  For  our  exanple,  we  find 
trial  and  error  that 


t'o.95  =  71*51 

and  thus. 


=  6.62 
e 

exar  95%  tolerance  limit  is  thus 

^0.95,0.95  T*  0.95  +  (Z0.95) 

=  81.40 

The  above  results  have  been  presented  for  the  case  of  linear  least  squares, 
that  is,  for  the  case  for  vrtiich  the  unknown  parameters  a  and  b  ajpear  linearly 
in  the  postulated  model.  The  results  can  be  generalized  in  a  natural  manner 
to  include  non-lineeir  cases.  The  details  of  this  extension  are  presented 
in  Hartley,  1964;  Grant,  1973;  and  Wallace  and  Grant,  1971. 

DRTft  ANALYSIS  WHEN  THE  UNDERLYING  DETEmtNISTIC  FACTORS  ARE  NOT  KNCT^JN 

In  the  preceeding  section,  we  examined  an  application  of  the  method  of 
least  squares  in  the  analysis  of  data.  In  that  exaitple,  we  assumed  that  the 
form  of  the  under  1>  ing  deterministic  mechanism  was  known.  Hiis  assunption 
allowed  the  analysis  to  be  performed  by  a  straightforward  separation  of  the 
<*served  responses  into  a  deterministic  and  randem  conponents.  The  statistical 
pre^^erties  of  the  parameter  estimates  and  the  residual  errors  then  allowed  the 
oenputation  of  statistical  estimates  of  the  system  response. 

In  this  section,  we  assvme  that  the.  form  of  the  deterministic  mechanism  is 
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unkncMn,  and  ve  examine  techniques  of  analysis  vrfiich  allow  us  to  infer  certain 
properties  of  the  deterministic  mechanisms  as  well  as  the  stochastic  properties 
of  the  system.  Rrom  these  inferences,  we  then  v»uld  postulate  forms  of  these 
processes  and  proceed  with  our  analysis  as  before. 

Whenever  we  are  atteupting  to  analyze  a  set  of  data,  we  generally  find 
that  the  data  set  has  two  properties  w^ich  constrain  our  techniques  of  analysis. 
First,  although  the  data  may  be  representative  of  a  continuous  system,  we  will 
tave  knowledge  of  the  state  of  the  system  only  at  discrete  points.  Second, 
we  will  find  that  our  sequence  of  data  covers  only  a  finite  time  period,  so 
that  we  are  unaware  of  the  state  of  the  system  except  during  this  interval. 

Given  these  constraints,  v/e  must  first  decide  whetlier  to  treat  the  system 
being  monitored  as  a  discrete  system  or,  by  interpolating  between  the  given 
data  points,  to  atterrpt  to  approximate  the  continuous  nature  of  the  ^stem. 

Ihe  restriction  imposed  by  the  finite  length  of  our  sanple  requires  us  to  treat 
our  system  as  periodic  with  a  pjeriod  equal  to  our  record  length,  or  to  extend 
the  record  based  upon  known  or  assumed  properties  of  the  system. 

In  this  discussion,  v/e  will  confine  ourselves  to  an  examination  of 
techniques  of  the  spectral  einalysis  of  discrete  periodic  systems.  Extensions 
of  these  techniques  to  continuous  aperiodic  systems  are  presented  in  Kisiel, 

1969,  and  to  discrete  apjeriodic systems  in  Kisiel,  1969,  and  fppenheim  and 
Schafer,  1975. 

Suppose  we  have  a  sequence  of  N  observations  Xn,  n  =  0,  1,...,  N-1. 

We  may  write  the  Fourier  transform  of  this  secjuence  (Cppenheim  and  Schafer,  1975) 
as  follows: 

X(k)  «  X-  EXP  £7-  j  (2ji)nk7 
n=o  “  N 

Where  the  X  (k)  are  the  Ftourier  coefficients  of  the  series  representation  of 
the  sequence  x^^,  that  is 

X(k)  EXP^  (2ir)nk3 
k=o  ^ 

Statistical  properties  of  the  sequence  x_  vAiich  are  of  interest  in  data  analysis 
include  the  mean,  the  variance »  and  the  autocorrelation  and  autovariance  functions. 
Hiese  quantities  are  defined  as  follows; 

Mean  =  E  (xn)  for  each  n. 

Variance  =  E  (xl  -  for  each  n, 

n 

Autocorrelation  sequence 

A  (n,m)  =  e(x  .  x„) , 

^xx  n 

Autooovariance  sequaice 

^xx  ®  O’^n  ”  *  ^^m  ~  “  ^xx  ^ 
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For  a  stationary  process,  we  may  write,  since  means  and  variances  are 
independent  of  time. 


if  =  E  (x^) 

'  X  =  E  {x„  -  X)  2 


^  (m)  =E  (V^+m) 
Yxx  =  ^xx  ^ 


For  ergodic  sequences,  we  nay  estimate  tJie  above  averages  as  the  time  averages 
of  the  single  sequence.  Various  formulas  have  been  proposed  to  ccrpuire 
estirates  of  the  averages  of  a  data  sequence.  For  example,  means  ard  variances 
nay  be  oonputed  by  the  formulas: 


X 


tJ-l 

1  I 
N  o 


Tlie  autocovariance  sequence  may  be  approxinated  by  the  relation 

k  (k)  =(N-K)“^  '!f‘^(x  -JJ)  (X  ,  -  x) 

^  i''=o  ^ 

and  the  autocorrelation  function  by  tiie  relation 

Kx(k)  =  (k)  +  X  2. 

Ihe  spectral  density  function  of  a  sample  is  defined  as  the  Fourier 
transform  of  the  autocovariance  sequence.  An  unsmoothed  estimate  of  the  one¬ 
sided  spectral  density  function  is  given  by 

G  (f)  =  2  (k  (o)  +  2  E  K  (r)  cos  (irrf)  +  K  (m)  cos  .tTrmf)  f 
M  r=l  ^  f  c  fcj 

Where  1/fc  equals  twice  the  sampling  frequency,  Tliis  critical  frequency,  known 

as  the  t/yquist  frequency,  determines  the  behavior  of  the  analysis.  Frequencies 
higher  thii  the  critical  frequency  are  not  visible  to  the  analysis;  however, 
the  variance  associated  with  these  frequencies  may  be  reflected,  or  "aliased" 
into  the  visible  spectrum. 

Itie  above  estinate  of  tlie  spectral  density  function  is  not  particularly 
suitable  for  analysis  because  it  is  unstable  witJi  respect  to  sample  size. 
Various  techniques  of  obtaining  a  smoothed,  or  stable,  estimate  are  present^ 
in  Oppesnheim  and  ScJiafer,  1975,  and  a  rath^  simple  procedure  is  presented  in 
Kisiel,  1969. 

nie  primary  purpose  of  spectral  analysis  of  monitoring  data  is  to  allcw 
seme  separation  of  the  mechanisms  of  variance  in  the  data,  and  to  simplify 
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the  statistical  properties  of  the  data.  V-tieii  aporopriate^  and  properly  used, 
spectral  analysis  can  te  a  powerful  tool  to  aid  in  the  interpretation  of  data. 
The  use  of  specttal  analysis  Ccin  be  illustrated  by  tlie  follcwing  sinple  exairple. 
Si;^pose  we  sanple  a  system  gejierated  by  tlie  process 

x(t)  =5  sin  (2iTt) 


at  intervals  of  0.25.  Wc  obtain  the  sequence  0,5, 0,-5. - .  Suppose  that  we 

have  a  sequence  of  13  such  data  points  which  we  use  to  estimate  the  spectral 
density  function  of  tl^e  system.  By  the  above  equations,  we  find 

K  =  11.54 
o 

Kj^  =  0 

K2  =-11.32 


K4  =  11.1],,. 


and 


e  (f)  =  0.5(11.54  -  22.72  TTf  +  11.11  cos  2tt3. 

The  spectral  density  function  is  to  be  evaluated  at  the  discrete  frequencies 
f  =  0.5k,  and  we  find 


G  (0)  =  0.04 
G  (0.5)  =  0.22 
G  (1)  =  22.69 
G  (1.5)  =  0.22 
G  (2)  =  -0.04. 

A  plot  of  the  \'ariance  spectra  is  shown  in  Figure  2.  The  preponderance  of 
the  variance  associated  v/ith  f  =  1  indicates  that  our  data  saiple  was  obtained 
from  a  periodic  process  with  period  equal  to  one  year.  Of  course,  spectral  density 
functions  of  real  data  usually  are  more  coiplicated  than  the  above  example,  and 
the  interpretations  of  tne' variance  spectra  is  not  so  straightfcx^^ard.  Examples 
of  interpretation  of  more  complicated  spectra  are  presented  in  Kisiel,  1969. 

Once  the  deteministic  contributions  of  a  data  sequence  cure  identified,  the 
contributions  of  these  deterministic  components  can  be  ranovod  frem  the  data. 

The  residuals  obtained  by  this  process  represent  the  random  conponents  of  the 
system,  and  may  be  cinalyzed  by  standard  satistical  methods.  The  results  of  such 
analysis  allow  the  syntiiesis  of  the  physical  system  and  a  prediction  of  extreme 
conditions  which  my  occur  within  the  system  in  the  future. 
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OOWIENTS  ON  STJ-IPLING  FREQUENCIES 


The  above  discussion  illustrates  the  requirements  of  a  rational  water 
quality  sairpiing  prxigram.  The  frequency  of  sanpling  must  satisify  the  Nyquist 
relation  for  the  highest  frequency  vhich  the  subsequent  data  analysis  will 
consider.  Because  natural  systems  rcurely  are  governed  by  a  pure  harmonic 
process,  sairpiing  at  a  multiple  of  two  or  three  times  the  Nyquist  frequency  is 
reoomnended  (Kieso.l,  1969). 

Frequencies  too  high  to  be  visible  to  the  analysis  of  a  data  set  either 
will  be  folded  back  into  the  visible  frequency  range,  a  phenomenon  known  as 
aliasing,  or  will  appear  as  randan  noise  to  the  analysis.  Figure  3  shows  a 
conceptual  diagram  of  the  aliasing  of  a  high  frequency  process  as  a  result  of 
widely-spaced  sampling. 

Requirements  upon  sample  size  and  the  lejngth  of  the  autocovariance  sequence 
are  dictated  hy  the  statistical  properties  of  the  estimates  of  the  averages  and 
tlie  spectral  density  function  of  the  data.  The  maximum  number  of  serial  lags, 
m,  used  in  computations  generally  should  be  10%  of  the  sanple  size  or  less. 

For  a  specified  standard  error  of  estimate  e  of  the  spectral  density  function, 
the  sanple  size  required  is  N  =  m/e^,  and  the  length  of  record  required  is 
NAT  (Kiesel,  1969). 


SU^IMAPY 

Statistical  analysis  of  \^ater  quality  data  usually  is  performed  in  order 
to  utilize  the  information  contained  within  the  data  set  to  obtain  an 
understanding  of  the  system  from  vhich  the  data  were  obtained.  Natural 
systems  generally  exhibit  characteristics  which  are  in  part  deterministic  and 
in  part  random.  A  rational  analysis  of  water  quality  data  requires  that  the 
deterministic  elements  of  the  system  be  identified  and  their  influences 
deleted  from  the  data.  The  stochastic  and  random  conponents  of  the  system  then 
can  be  characterized  statistically. 

Techniques  for  the  analysis  of  data  exhibiting  deterministic  and  stochastic 
properties  have  been  discussed,  and  sinple  examples  of  the  use  of  the  methods 
have  been  presented.  The  implications  of  the  analysis  techniques  presaited 
are  discussed  in  terms  of  the  requirements  which  these  techniques  impose  upon 
data  collection  programs. 
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SPECTRAL  DENSITY  FUNCTION 


f  {CYCLES/ YEAR) 


ALIASING  OF 

HIGH  FREQUENCY  SIGNALS 


DATA  VALIDATION 

By 

James  D.  Ashworth*^ 


Introduction 


This  paper  is  a  general  presentation  of  the  topic  of  environmental  data 
validation.  It  is  subdivided  into  theoretical  considerations,  methodo¬ 
logy,  and  procedures  presently  in  use.  Neither  computer  software  nor 
hardware  are  discussed.  The  reader  is  advised  that  terms  used  are 
defined  in  the  Appendix. 

Theoretical  Considerations 


Intuitively,  data  validation  must  rely  upon  the  organization  of  records 
into  groups  based  upon  common  attributes  and  Independent  components.  The 
validating  process  then  requires  the  ability  to  discern  known  Intra-group 
properties  and  relationships.  Group  boundaries  are  not  crossed  during  any  test. 

In  general,  records  resulting  from  environmental  monitoring  contain  five 
independent  components;  i.e.,  parameter  name,  time,  and  3  dimensional 
location.  Consequently,  the  theoretical  number  of  group  types  or 
classes  is  2^  ■  32.  Figure  1.  is  a  summarization  of  the  theoretical 
grouping  and  testing  process.  Each  row  represents  an  organization  for 
all  of  the  data  to  be  validated. 

Confining  this  discussion  to  non-redundant  records,  the  Class  1  groups 
contain  one  record  each.  Intra-record  tests  of  syntax,  component  logic, 
and  reference  comparison  are  possible.  The  Class  II  groups  contain  a 
complete  sample  permitting  tests  for  balances  and  inter-parameter  re¬ 
lationships.  Class  III  groups  can  be  chosen  by  varying  any  one  of  the 
four  space-time  components  %diile  holding  the  other  three  as  well  as 
parameter  name  constant.  This  class  enables  testing  for  continuity. 

Classes  I  thru  III  constitute  a  comprehensive  validating  effort.  Higher 
classes  are  confounded  by  multiple  variables  rendering  relationships 
difficult  to  define  and  cumbersome  for  validating  purposes.  In  fact 
these  data  groupings  are  used  in  sophisticated  modeling  efforts,  a  major 
purpose  for  the  creation  of  the  data  base. 

Testing  is  normally  performed  in  class  order  with  concurrent  editing. 

No  edit  can  be  made  without  a  manual  effort  involving  a  multiple  path 
decision.  Ultimately,  substantial  evidence  of  a  mistake  must  be  found 
or  exceptions  must  be  accepted.  Class  I  exceptions  are  most  easily 
corrected  since  testing  concentrates  on  the  independent  portion  of  the 
record.  All  higher  classes  rely  on  value  testing.  While  Classes  I  and 
II  exceptions  tend  to  be  distributed  randomly  throughout  data  with 


1  Water  Quality  Section,  HunClngtoo  District 


Paper  6 


respect  to  all  components,  Class  III  exceptions  tend  to  be  outliers. 
Consequently,  incorrect  Class  III  editing  decisions  can  bias  the  data 
base.  Finally,  the  value  range  test  which  is  often  associated  with 
Class  I  is  legitimately  a  part  of  classes  III  or  higher  inasmuch  as  the 
range  is  determined  by  experience. 

Methodology 

Tests  of  record  syntax  search  for  illegal  characters  based  on  the  allow¬ 
able  syntax.  Examples  of  component  logic  include  determination  that 
recorded  year,  month,  day,  hour,  and  minute  do  not  exceed  current  year, 
twelve,  the  maximum  for  the  month,  23,  and  59,  respectively.  Reference 
comparisons  can  be  made  for  station  number,  date,  time,  depth,  para¬ 
meter,  and  value.  Within  Class  I,  however,  the  value  range  is  based  on 
physical  constraints  and  not  experience. 

Intra-sample  balances  for  charge  and  total  and  dissolved  constituents 
are  appropriate  Class  II  tests  as  well  as  alkalinity-acidity-pH, 
specific  conductance-ionic  strength,  and  dissolved  oxygen-temperature 
relationships. 

Continuity  can  be  detected  computationally  or  graphically.  Such  testing 
is  based  upon  the  assumption  that  nature  rarely  produces  a  discontinuous 
function.  The  major  prerequisite  is  high  data  density  with  respect  to 
the  fixed  independent  components.  A  common  application  is  the  plot  of 
depth  versus  concentration  holding  horizontal  location  and  time 
constant.  Computational  approaches  rely  upon  changes  in  the  slope  of 
the  chosen  function.  It  should  be  stressed  that  since  Class  III  tests 
seek  outliers,  incorrect  application  can  bias  the  data  base. 

Procedures  In  Use 


The  validating  procedures  used  in  the  maintenance  of  four  large  data 
bases  are  outlined  in  Figure  2.  It  is  seen  that  data  collected  with  the 
aid  of  robot  monitors  receives  somewhat  different  attention  than  that 
collected  by  manual  methods.  This  is  primarily  the  result  of  three 
factors,  higher  data  but  lower  parameter  densities  and  decreased 
opportunity  for  transcription  errors,  where  data  is  collected 
mechanically.  Consequently,  the  Geological  Survey  often  tests  for  con¬ 
tinuity  of  Class  II  relationships  corresponding  to  line  10  of  Figure  1 
at  monitored  stations,  but  for  Ionic  Balance,  Class  II,  of  grab  samples 
assumed  to  measure  all  constituents. 

The  proposed  validating  system  of  the  National  Climatic  Center  will  use 
asymptotic  singular  decomposition  of  a  matrix  to  streamline  Classes  II 
and  III  relational  testing  as  well  as  provide  data  compaction.  Value 
range  limits  based  on  physical  constraints  are  contrasted  with  those 
based  on  experience  by  testing  procedures  for  records  from  Orsanco 
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THEORETICAL  MATRIX  OF  RECORD  GROUPING  CHOICES  FOR  DATA 
VALIDATION  TESTING 
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The  Storet  System  is  unique  as  a  central  repository  for  data  from  many 
sources.  It  is  rightly  termed  a  collection  of  data  bases.  The  valida¬ 
ting  tests  listed  are  optional.  The  value  range  tests  became  part  of 
the  system  in  1975.  They  consist  of  user  supplied  station  or  run 
specific  value  ranges.  Application  of  this  logic  with  a  single  range  to 
all  data  accumulated  prior  to  1975,  approximately  30  million  data 
points,  identified  approximately  45,000  exceptions.  To  date,  these  have 
not  been  edited. 


VALIDATING  TESTS  IN  USE 
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FIG 


HIERARCHAL  GLOSSARY 


1.  Parameter  -  Name  of  any  measured  quantity;  i.e.,  dissolved  oxygen 
pH,  specific  conductance,  etc.;  an  independent  component  of  a 
record. 

2.  Location  -  Designation  of  sampling  point  in  3  dimensional  space, 
often  times  an  arbitrarily  assigned  "station  number"  and  a  depth. 
The  location  consists  of  3  independent  components  of  a  record. 

3.  Time  -  Date  and  time  of  sample  collection,  an  independent  com¬ 
ponent  of  a  record. 

4.  Value  -  Digital  results  of  a  measurement,  the  dependent  component 
of  a  record. 

5.  Record  Identity  -  The  set  of  5  independent  components  of  a  record 
parameter,  3  dimensional  location  and  time. 

6.  Record  -  A  value  and  associated  identity. 

7.  Sample  -  Any  complete  record  set  with  common  location  and  time, 
parameter  is  varied.  Complete  is  interpreted  as  "all  which  were 
selected  for  determination". 

8.  Validation  -  Act  of  testing  for  correctness. 

9.  Exception  Record  -  Record  which  failed  some  test  for  validity. 

10.  Edit  -  Act  of  changing  a  record. 

11.  Data  -  Any  set  of  records. 

12.  Data  Base  -  The  set  of  all  historical  records. 
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TEMPERATURE  AND  DISSOLVED  OXYGEN 
DATA  INTERPRETATION 


By 


Dennis  E.  Ford* 


INTRODUCTION 


The  data  Interpretation  process  is  viewed  differently  by  people 
depending  on  their  particular  needs  and  interests.  A  statistician 
may  view  data  interpretation  only  in  terms  of  statistics  while  a  policy 
maker  may  view  it  as  a  definite  answer  to  the  question:  Does  a  par¬ 
ticular  water-quality  parameter  violate  standards?  In  this  paper, 
data  interpretation  is  considered  synonymous  with  a  basic  understanding 
of  all  factors  and  processes  responsible  for  or  Influencing  the  par¬ 
ticular  phenomenon  or  parameter  of  Interest.  A  logical  explanation 
therefore  exists  for  all  observed  phenomena. 

The  general  theme  of  this  paper  will  be  to  identify  some  of  the 
Important  factors  and  processes  Influencing  temperature  and  dissolved 
oxygen  (DO)  phenomena  and  to  illustrate  them  with  field  data.  This 
can  be  a  formidable  task  for  water-quality  data  because  of  the  numerous 
interactions  between  the  morphometry,  physics,  chemistry,  and  biology 
of  a  river  or  reservoir  system.  All  of  these  interactions  must  be  con¬ 
sidered  in  the  interpretation  process  for  all  water-quality  variables 
including  temperature  and  DO. 

Only  reservoirs  will  be  considered  in  this  paper.  Since  similar 
processes  occur  in  rivers,  some  of  the  information  and  techniques 
presented  herein  will  also  be  applicable  for  riverine  systems.  The 
general  order  of  presentation  will  be  reservoir  characteristics, 
temperature  data,  DO  data,  and  oxygen  demands. 


RESERVOIR  CHARACTERISTICS 


Morphometry  -  According  to  recent  estimates,  there  are  over  400 
authorized  CE  reservoirs.  These  reservoirs  vary  in  surface  area  from 
approximately  1  km^  to  10^  km^  and  in  volume  from  approximately  10^  m^ 
to  lOll  m^.  The  general  shape  of  most  reservoirs  is  probably  best 
described  as  dendritic.  In  such  a  heteorgeneous  group,  it  is  unlikely 


*  Hydrologist,  Environmental  Laboratory,  U.  S.  Army  Engineer  Waterways 
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that  similar  processes  dominate  In  reservoirs,  but  relationships 
do  appear  to  exist  between  basic  morphometric  parameters  such  as 
surface  area,  volume,  depth,  drainage  basin  size,  and  Inflow  (e.g.. 
Figures  1-3) .  These  relationships  are  encouraging  because  they  bring 
a  semblance  of  order  to  a  diverse  group  of  reservoirs.  These  morpho¬ 
metric  relationships  Indicate  that  reservoirs  with  similar  morphology 
may  exhibit  other  similarities  and  that  relationships  with  other 
parameters  such  as  thermocline  depth  and  mixing  coefficients  may 
exist  (Figures  4  and  5). 

An  important  factor  of  reservoir  morphometry  is  the  depth  to 
length  ratio  (i.e,  H/L) .  A  typical  value  Is  10  This  ratio  Implies 
that  problems  involving  vertical  variations  may  be  treated  Independently 
from  problems  involving  horizontal  variations.  It  also  explains  why 
one-dimensional  mathematical  reservoir  models  have  been  relatively 
successful  In  describing  reservoir  dynamics. 

Flow  Regime  -  The  flow  regime  within  a  reservoir  is  generally 
characterized  as  turbulent  and  unsteady.  Typical  Reynolds  numbers 
are  always  greater  than  10^,  clearly  indicating  turbulent  flow.  In 
stratified  reservoirs,  where  buoyancy  forces  dominant,  the  Reynolds 
number  Is  not  a  valid  criterion  for  turbulence.  Instead,  a  Richardson 
number  or  densimetric  Froude  number  Is  used.  These  criteria  indicate 
sporadic  turbulence  in  all  but  the  most  stable  density  gradients. 
Generally,  the  turbulence  is  nonisotropic  because  of  the  small  H/L. 

This  Implies  that  the  reservoir  is  not  homogeneous,  and  that  obtaining 
representative  water-quality  samples  may  be  difficult. 

The  flow  regime  is  unsteady  because  the  controlling  or  driving 
forces  are  unsteady.  For  example,  daily  variations  in  wind,  solar 
radiation,  and  inflow  are  shown  in  Figures  6  and  7.  Diurnal,  nocturnal, 
and  hourly  variations  also  occur.  The  response  of  a  lake's  water 
surface  temperature  and  mixed  layer*  depth  to  daily  variations  in 
solar  radiation  and  wind  speed  is  shown  in  Figure  6.  Under  constant 
forcing,  it  would  take  several  days  for  a  medium  size  reservoir  to 
attain  steady  state  conditions.  Such  steady  state  conditions  are 
unlikely  in  reservoirs. 

Reservoir  Natural  Lake  -  There  are  several  reasons  why  a  reser¬ 
voir  may  differ  from  a  natural  lake.  Several  of  these  reasons  are 
discussed  by  Baxter  (1977).  One  obvious  reason  is  some  reservoirs  have 
the  capability  to  release  water  from  different  depths  (i.e.,  selective 
or  bottom  withdrawal)  rather  than  only  at  the  water  surface,  as  in  a 
natural  lake.  Reservoirs  also  have  shorter  residence  times,  larger 
drainage  basins,  larger  shoreline  development  ratios,  etc.  The  impact 
of  each  of  these  differences  on  the  water  quality  of  a  reservoir  is 
beyond  the  scope  of  this  paper;  however,  at  this  time  it  is  only 
necessary  to  understand  that  differences  do  exist  and  that  extrapolating 


^  The  mixed  layer  is  defined  later  in  the  text  and  in  Figure  8. 
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results  from  natural  lakes  to  reservoirs  may  be  subject  to  error. 

A  natural  lake  may,  however,  be  a  simplified  subset  of  a  reservoir. 
In  many  instances.  Important  processes  are  much  easier  to  identify 
in  natural  lakes  since  there  are  fewer  complicating  factors.  Some 
processes  from  natural  lakes  will  therefore  be  used  to  Illustrate 
important  points  in  this  paper. 


TEMPERATURE  DATA 


Of  all  the  water-quality  parameters,  temperature  data  is  perhaps 
the  easiest  to  obtain  and  the  most  useful.  Since  all  chemical  and 
biological  processes  are  temperature  dependent  to  varying  degrees  and 
since  these  processes  are  also  governed  by  the  degree  of  stratification 
and  mixing,  temperature  measurements  are  essential  to  any  water-quality 
study.  A  reservoir  is  characterized  as  being  of  poor  or  good  water 
quality  based  on  many  other  variables  in  addition  to  temperature. 

Thermal  Stratification  -  Before  discussing  thermal  stratification, 
it  is  important  to  review  some  basic  terminology.  Several  basic 
limnology  books  describe  stratification  in  detail  (e.g.,  Wetzel,  1975; 
Ruttner,  1963;  Hutchinson,  1957).  The  important  terms  are  illustrated 
in  Figure  8.  The  epllimnion  is  the  upper  strata  of  warm,  turbulent 
water.  It  is  usually  characterized  by  relatively  uniform  temperatures. 
The  deep,  cold,  relatively  undisturbed  region  is  termed  the  hypolimnion. 
Between  the  epllimnion  and  hypolimnion  is  the  metalimnion  which  is 
characterized  by  a  strong  temperature  (density)  gradient.  The  plane 
of  inflection  or  of  maximum  temperature  gradient  is  termed  the  thermo- 
cline.  Other  definitions  have  been  proposed  for  the  thermocline 
(e.g.,  the  l°C/m  criterion),  but  they  are  not  in  wide  use  and  should 
be  avoided.  Another  term  illustrated  in  Figure  8  is  the  "mixed  layer." 
The  mixed  layer  is,  as  Implied  by  the  name,  the  overlying  isotropic 
layer.  Since  the  mixed  layer  refers  to  the  instantaneous  depth  of  the 
overlying  isotropic  layer,  it  differs  from  the  epllimnion,  which  is 
an  averaged  mixed  layer,  in  two  respects.  First,  its  depth  is  usually 
1  ss  than  the  depth  of  the  epllimnion.  Second,  it  is  much  more  dynamic 
than  the  epllimnion. 

There  are  many  other  terms  used  to  describe  the  thermal  structure 
of  a  lake.  These  can  be  found  in  the  above  references  and  are  related 
to  the  number  of  turnovers  (i.e.,  periods  of  complete  vertical  mixing) 
occurring  within  a  lake  and  to  the  strength  of  the  stratification. 

For  the  most  part,  these  terms  are  of  minor  importance  compared  to  the 
difference  between  holomictic  and  meromictic  lakes.  In  holomictic 
lakes,  the  entire  water  column  completely  circulates  or  turns  over. 

Lakes  that  cannot  circulate  completely,  and  exhibit  a  deep  strata 
that  is  perennially  stagnant  in  the  water  column,  are  termed  meromictic 
lakes.  The  reason  for  this  condition  may  be  physical,  chemical,  or 
biological.  In  any  case,  it  can  have  a  profound  effect  on  the 
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temperature  and  mixing  structure  and  consequentially  on  the  water 
quality  of  a  lake.  Bottom  withdrawal  is  an  effective  means  to  elimi¬ 
nate  this  undesirable  condition. 


Factors  Affecting  Stratification  -  Since  all  lakes  exhibit  some 
degree  of  thermal  stratification,  it  is  important  to  Investigate  the 
factors  responsible  for  stratification.  The  principal  factors  influencing 
the  formation,  strength,  and  extent  of  stratification  are  the  density 
of  water,  solar  radiation  and  the  heat  budget  at  the  air-water  interface, 
and  the  mixing  resulting  from  advection  and  wind-induced  phenomena. 

It  is  well  known  that  the  density  of  water  varies  with  temperature 
(Figure  9).  The  importance  of  this  variation  in  determining  the 
distribution  of  heat  within  a  lake  was  originally  documented  by 
Blrge  (1910).  Two  factors  are  important.  First,  the  maximum  density 
of  water  occurs  at  4^0.  Colder  or  warmer  water  will  be  less  dense  and 
tend  to  float  on  top  of  the  heavier  4°C  water.  Second  and  more  important, 
the  density  of  water  decreases  with  an  escalating  rate  with  increasing 
temperatures  above  4°C.  Therefore,  the  buoyant  forces  and  resistance 
to  mixing  increases  with  temperature.  There  is  over  an  order  of 
magnitude  difference  in  the  energy  requirements  to  mix  a  1°C  temperature 
difference  at  25°C  than  at  5°C. 

The  energy  available  to  warm  the  waters  of  a  reservoir  ultimately 
comes  from  solar  radiation  which  varies  seasonally.  The  seasonal 
variation  of  solar  radiation  follows  a  sinusoidal  curve  with  a  maximum 
in  late  June.  In  addition,  diurnal  cycles  also  occur.  Water  temperatures 
respond  to  both  of  these  cycles  with  a  slight  delay. 

Solar  radiation  is  absorbed  at  the  water  surface  and  selectively 
with  depth  depending  on  the  wavelength  of  the  light,  properties  of  the 
water,  and  the  matter  suspended  in  the  water.  This  absorption  is 
usually  assumed  to  be  exponential  with  depth  (i.e..  Beer's  Law),  but 
surface  effects  result  in  minor  discrepancies  in  the  top  meter  or  so 
of  a  lake  (Figure  10).  To  avoid  these  discrepancies,  it  is  sometimes 
assumed,  especially  by  mathematical  modelers,  that  a  certain  fraction 
(g)  of  the  solar  radiation  is  absorbed  in  a  surface  layer  of  a  speci¬ 
fied  thickness  (e.g.,  2  ft  (0.6  m))  and  the  remaining  fraction  is 
absorbed  exponentially  with  depth  with  an  extinction  coefficient,  q. 

The  specific  magnitudes  of  6  and  q  depend  on  the  properties  of  the 
water,  but  there  is  a  definite  relationship  between  6  and  n  (Figure  11). 
For  physical  significance,  the  3  and  n  used  in  mathematical  model 
studies  should  reflect  this  relationship.  If  field  measurements  are 
not  available,  the  extinction  coefficient  can  be  estimated  directly 
from  Secchi  disk  measurements  (Figure  12) .  Considering  all  the  un¬ 
certainties  that  go  into  a  Secchi  disk  measurement,  the  relationship 
depicted  in  Figure  12  is  excellent. 

Field  measurements  are  also  remarkably  consistent  in  predicting 
the  Secchi  disk  depth  to  occur  at  the  10-15  percent  light  level 
(Ford,  1976).  This  fact  is  useful  in  calculating  extinction  coefficients 

4 


Paper  7 


I 


directly  from  Secchl  disk  measurements.  For  example,  assuming  exponential 
absorption  with  depth,  the  vertical  distribution  of  light  with  depth  Is 


where  Is  the  light  Intensity  at  depth  z 

Is  the  light  Intensity  at  the  water  surface. 

Now,  If  the  Secchl  disk  depth  Is  assumed  to  occur  at  the  12.5  percent 
light  level,  then  I  /I  Is  0.125  and  the  extinction  coefficient  Is  given 
by  ^  ° 

n  =  2.1/SD 


where 


SD  Is  the  Secchl  disk  depth  In  meters 
n  Is  the  extinction  coefficient  In  1/m. 

The  problem  with  using  these  estimates  of  light  extinction  Is  the 
penetration  of  light  Is  dependent  on  algae  concentrations,  suspended 
solids,  etc.,  and  highly  variable  (Figure  13).  Some  Investigators 
have  been  successful  In  relating  Secchl  disk  depths  to  a  suspended 
solids  concentration  (F.  Schlebe,  personal  communication,  1978)  and 
chlorophyll  a  (Carlson,  1977).  Comparison  of  Figures  13  and  14  shows 
how  Important  light  penetration  Is  In  determining  the  depth  of  the 
thermocline.  The  thermocllne  depth  Is  deeper  for  the  years  1968-72 
when  the  Secchl  disk  depths  are  larger. 

In  contrast  to  heating,  cooling  of  a  water  body  can  only  occur 
at  the  water  surface.  It  Is  therefore  possible  for  the  surface  water 
to  decrease  In  temperature  while  deeper  water  Increases  In  temperature. 
If  the  temperature  of  the  surface  water  drops  below  the  temperature  of 
the  deeper  water  and  still  remains  above  A^C,  the  water  column  becomes 
thermally  unstable  and  natural  convection  and  mixing  commences.  The 
convection  mixing  Is  termed  penetrative  convective  mixing  If  It  results 
In  the  entrainment  of  water  from  an  underlying  region  of  stable  density 
(temperature)  gradient.  Otherwise,  It  Is  termed  nonpenetratlve  con¬ 
vective  mixing  and  mixing  occurs  only  down  to  a  depth  to  eliminate 
density  Instabilities.  It  Is  generally  believed  that  convective  mixing 
In  lakes  and  reservoirs  Is  penetrative  although  most  mathematical 
models  consider  It  to  be  nonpenetratlve. 

Mixing  In  lakes  and  reservoirs  results  from  the  cumulative  effects 
of  Inflows,  outflows,  wind  generated  currents,  surface  and  Internal 
waves,  Langmuir  circulation,  natural  convection,  etc.  Since  the 
contribution  of  these  phenomena  and  their  Interactions  to  mixing 
are  not  completely  understood,  their  cumulative  effects  are  usually 
lumped  Into  an  eddy  or  effective  dlffuslvlty  term.  This  approach 
Is  commonly  used  In  mathematical  modeling. 
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The  eddy  diffusion  approach  Is  based  on  an  analogy  with  molecular 
diffusion.  In  molecular  diffusion,  a  property  Is  transferred  down 
a  concentration  gradient  by  the  random  motion  of  molecules  with  no 
overall  transport  of  the  fluid.  The  flux  or  transport  of  a  property 
Is  therefore  equal  to  a  molecular  diffusion  coefficient  times  a  concen¬ 
tration  gradient.  That  Is 

Flux  =  K 

mdx 

where 


K  Is  the  molecular  diffusion  coefficient 
m 


■K—  is  the  concentration  gradient. 


The  eddy  diffusion  approach  also  assumes  that  the  turbulent  flux 
of  a  constituent  is  proportional  to  an  eddy  diffusion  coefficient  times 
a  gradient.  The  eddy  diffusion  coefficients  are  usually  assumed  to  be 
much  larger  than  the  molecular  diffusion  coefficient.  One  important 
distinction  of  the  molecular  diffusion  coefficient  is  a  property  of  the 
fluid  while  the  eddy  diffusion  coefficient  depends  on  the  flow,  density 
stability,  etc.,  of  the  fluid.  A  simple  example  will  illustrate  the 
implications  of  this  fact.  It  is  commonly  recognized  that  the  strong 
gradients  in  the  metalimnion  of  a  lake  inhibit  turbulence  and  mixing. 
Since  the  eddy  diffusion  approach  assumes  the  transport  of  a  constituent 
is  proportional  to  the  concentration  gradient,  mixing  or  transport 
will  be  greatest  in  the  metalimnion  where  the  gradients  are  the  largest. 
This  conclusion  is  contrary  to  field  observations.  Extreme  care  must 
therefore  be  used  in  analyzing  data  based  on  this  approach. 

Eddy  diffusion  coefficients  can  be  back  calculated  from  tempera¬ 
ture  data  and  the  one-dimensional,  unsteady,  thermal  energy  equation 
(l.e.,  diffusion  equation) 

^  ^  3  K  ^ 

9t  ^  ^  3z  3z 

where  T  is  the  mean  temperature 

z  is  the  depth 

is  the  eddy  diffusion  coefficient 

H  is  the  internal  flux  of  solar  radiation 
z 


If  this  equation  is  integrated  down  the  water  column  from  depth  z  to 
the  bottom  depth  H,  then 


K 


z 


This  equation  can  be  solved  graphically  or  with  the  aid  of  a  computer. 
Typical  results  are  shown  in  Figure  15.  Variations  with  time  and 
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directly  from  Secchl  disk  measurements.  For  example,  assuming  exponential 
absorption  with  depth,  the  vertical  distribution  of  light  with  depth  Is 


where  Is  the  light  Intensity  at  depth  z 

Is  the  light  Intensity  at  the  water  surface. 

Now,  If  the  Secchl  disk  depth  Is  assumed  to  occur  at  the  12.5  percent 
light  level,  then  I  /I  la  0.125  and  the  extinction  coefficient  Is  given 
by 

n  -  2.1/SD 


where 


SD  Is  the  Secchl  disk  depth  In  meters 
n  Is  the  extinction  coefficient  In  1/m. 

The  problem  with  using  these  estimates  of  light  extinction  Is  the 
penetration  of  light  Is  dependent  on  algae  concentrations,  suspended 
solids,  etc.,  and  highly  variable  (Figure  13).  Some  Investigators 
have  been  successful  In  relating  Secchl  disk  depths  to  a  suspended 
solids  concentration  (F.  Schlebe,  personal  communication,  1978)  and 
chlorophyll  a  (Carlson,  1977) .  Comparison  of  Figures  13  and  14  shows 
how  Important  light  penetration  Is  in  determining  the  depth  of  the 
thermocline.  The  thermocllne  depth  is  deeper  for  the  years  1968-72 
when  the  Secchl  disk  depths  are  larger. 

In  contrast  to  heating,  cooling  of  a  water  body  can  only  occur 
at  the  water  surface.  It  is  therefore  possible  for  the  surface  water 
to  decrease  In  temperature  while  deeper  water  Increases  In  temperature. 
If  the  temperature  of  the  surface  water  drops  below  the  temperature  of 
the  deeper  water  and  still  remains  above  40C,  the  water  column  becomes 
thermally  unstable  and  natural  convection  and  mixing  commences.  The 
convection  mixing  is  termed  penetrative  convective  mixing  If  It  results 
In  the  entrainment  of  water  from  an  underlying  region  of  stable  density 
(temperature)  gradient.  Otherwise,  It  is  termed  nonpenetratlve  con¬ 
vective  mixing  and  mixing  occurs  only  down  to  a  depth  to  eliminate 
density  Instabilities.  It  Is  generally  believed  that  convective  mixing 
in  lakes  and  reservoirs  is  penetrative  although  most  mathematical 
models  consider  It  to  be  nonpenetratlve. 

Mixing  In  lakes  and  reservoirs  results  from  the  cumulative  effects 
of  Inflows,  outflows,  wind  generated  currents,  surface  and  Internal 
waves,  Langmuir  circulation,  natural  convection,  etc.  Since  the 
contribution  of  these  phenomena  and  their  Interactions  to  mixing 
are  not  completely  understood,  their  cumulative  effects  are  usually 
lumped  Into  an  eddy  or  effective  diffusivity  term.  This  approach 
Is  commonly  used  In  mathematical  modeling. 
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The  eddy  diffusion  approach  Is  based  on  an  analogy  with  molecular 
diffusion.  In  molecular  diffusion,  a  property  Is  transferred  down 
a  concentration  gradient  by  the  random  motion  of  molecules  with  no 
overall  transport  of  the  fluid.  The  flux  or  transport  of  a  property 
Is  therefore  equal  to  a  molecular  diffusion  coefficient  times  a  concen¬ 
tration  gradient.  That  Is 


where 


Flux 


K  Is  the  molecular  diffusion  coefficient 
in 


3c 


Is  the  concentration  gradient. 


The  eddy  diffusion  approach  also  assumes  that  the  turbulent  flux 
of  a  constituent  Is  proportional  to  an  eddy  diffusion  coefficient  times 
a  gradient.  The  eddy  diffusion  coefficients  are  usually  assumed  to  be 
much  larger  than  the  molecular  diffusion  coefficient.  One  Important 
distinction  of  the  molecular  diffusion  coefficient  Is  a  property  of  the 
fluid  while  the  eddy  diffusion  coefficient  depends  on  the  flow,  density 
stability,  etc.,  of  the  fluid.  A  simple  example  will  Illustrate  the 
Implications  of  this  fact.  It  Is  commonly  recognized  that  the  strong 
gradients  In  the  metallmnlon  of  a  lake  Inhibit  turbulence  and  mixing. 
Since  the  eddy  diffusion  approach  assumes  the  transport  of  a  constituent 
Is  proportional  to  the  concentration  gradient,  mixing  or  transport 
will  be  greatest  In  the  metallmnlon  where  the  gradients  are  the  largest. 
This  conclusion  Is  contrary  to  field  observations.  Extreme  care  must 
therefore  be  used  In  analyzing  data  based  on  this  approach. 


Eddy  diffusion  coefficients  can  be  back  calculated  from  tempera¬ 
ture  data  and  the  one-dlmenslonal,  unsteady,  thermal  energy  equation 
(l.e.,  diffusion  equation) 

il  =  3  K  W  -  ^^z 
9t  ^  ^  9z  9z 

where  T  Is  the  mean  temperature 

z  Is  the  depth 

Is  the  eddy  diffusion  coefficient 

H  Is  the  Internal  flux  of  solar  radiation 
z 

If  this  equation  Is  Integrated  down  the  water  column  from  depth  z  to 
the  bottom  depth  H,  then 


This  equation  can  be  solved  graphically  or  with  the  aid  of  a  computer. 
Typical  results  are  shown  In  Figure  15.  Variations  with  time  and 
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depth  (i.e.,  stratification)  are  evident.  If  daily  temperature  data 
are  used,  then  the  calculated  eddy  diffusion  coefficients  can  vary 
by  several  orders  of  magnitude  from  one  day  to  the  next.  It  should 
be  noted  that  this  approach  does  not  consider  horizontal  advectlon. 

Since  some  reservoirs  are  dominated  by  horizontal  advectlon,  the  use 
of  this  approach  may  result  in  unrealistic  estimates  for  eddy  diffusion 
coefficients. 

One  application  of  eddy  diffusion  coefficients  is  to  calculate 
the  flux  of  a  nutrient  (e.g.,  phosphorus)  from  the  metallmnion  to  the 
epilimnion.  The  transport  of  mass  and  momentum  are  not  necessarily 
equal  to  the  transport  of  thermal  energy.  Since  most  dlf fusivlties 
are  based  on  the  transport  of  thermal  energy  (i.e. ,  back,  calculated 
from  temperature  data) ,  their  application  to  other  constituents  such  as 
dissolved  oxygen  or  nutrients  may  be  questionable. 

Recent  sensitivity  analyses  on  an  ecological  model  at  the 
Waterways  Experiment  Station  (WES)  (Thornton  and  Lessem,  1976)  showed 
the  diffusion  coefficient  to  be  a  highly  sensitive  variable.  Pertur¬ 
bations  of  the  diffusion  coefficient  that  had  minimal  effect  on  the 
temperature  structure  had  dramatic  effects  on  the  algae  concentrations. 
One  reason  for  this  sensitivity  is  differences  between  concentration 
gradients.  Although  the  physical  location  (depth)  of  maximum  gradients 
may  be  similar  for  different  constituents  (e.g.,  the  metallmnion), 
the  magnitudes  of  the  gradients  are  not  similar.  For  example,  in  the 
metallmnion  of  a  lake,  temperatures  may  vary  between  10°  and  20^0 
while  nutrient  concentrations  may  vary  over  two  orders  of  magnitude. 

The  use  of  a  single  constant  eddy  diffusion  coefficient  to  determine 
the  transport  of  temperature  and  nutrient  in  this  example  would  result 
in  an  excessive  transport  of  the  nutrient. 

Seasonal  Temperature  Cycle  -  To  see  how  the  above  factors  interact 
to  form  a  seasonal  temperature  cycle,  a  natural  lake  will  be  used  as 
a  simple  example  (Figure  16)  to  minimize  the  effects  of  inflow  and 
outflow.  In  Figure  16,  the  ice  cover  on  the  lake  went  out  on  17  April. 
Thereafter,  the  entire  water  column  mixed  for  several  days  (i.e., 
period  of  spring  turnover).  The  density  differences  at  the  low  tempera¬ 
ture  were  not  sufficient  to  prevent  complete  mixing.  As  the  water 
column  warmed,  density  differences  Increased  and  it  became  more  difficult 
to  mix  the  entire  water  column.  Stratification  started  to  form  at  the 
bottom  of  the  lake  because  the  density  differences  and  resulting 
buoyancy  forces  were  small  compared  to  the  kinetic  energy  input  (i.e., 
wind).  As  the  solar  radiation  increased;  water  temperatures  increased; 
density  differences  increased;  and  the  thermocline  moved  upward  because 
the  kinetic  energy  input  could  not  overcome  the  ever  increasing 
buoyancy  forces.  The  minimum  thermocline  depth  was  usually  achieved 
sometime  around  the  time  of  summer  solstice  or  time  of  maximum  heat 
input.  Once  stratification  formed,  the  hypolimnlon  temperature  remained 
relatively  constant  until  fall  overturn.  Fall  overturn  occurred  when 
the  surface  temperature  approached  the  hypolimnlon  temperature. 
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Generally,  the  slope  of  the  Isotherms  indicates  the  degree  of 
mixing;  however,  this  assumes  no  internal  heating  from  solar  radiation 
occurs.  The  relatively  flat  Isotherms  in  the  metallmnlon  during  mid¬ 
summer  indicated  very  little  mixing  in  this  region  (Figure  16) . 
Similarly,  the  steep  slope  of  the  14°  and  16°C  isotherms  during  late 
May  and  early  June  should  indicate  intensive  mixing.  During  this 
period,  however,  the  Secchi  disk  depth  increased  from  2  m  to  6  m, 
indicating  the  heating  was  due  to  internal  absorption  of  solar 
radiation  and  not  mixing  (Figure  17). 

The  increased  mixing  resulting  from  inflows  and  outflows  in 
reservoirs  can  result  in  deviations  from  the  above  example.  For 
instance,  periods  of  overturn  may  be  extended  and  become  more  frequent, 
and  the  slope  of  the  isotherms  may  be  increased.  If  bottom  withdrawal 
is  used,  hypolimnetlc  temperatures  may  increase.  For  example,  in 
Figure  18,  the  hypolimnetlc  temperatures  in  Beltzville  Lake  in  1976 
remained  relatively  constant  when  only  small  amounts  of  water  were 
released  through  the  lower  ports.  In  contrast,  in  1972,  when  large 
quantities  of  water  were  released  through  the  flood  gates  in  response 
to  Hurricane  Agnes,  the  hypolimnetlc  temperatures  increased.  Man- 
induced  perturbations  such  as  power  operations  can  also  increase  mixing 
and  hypollmnion  temperatures  (Figure  19) . 

Factors  Affecting  Data  Interpretation  -  Several  other  factors 
need  to  be  considered  in  temperature  data  interpretation.  First, 
mixing  and  entrainment  are  not  gradual  continuous  processes  as  indi¬ 
cated  by  the  seasonal  temperature  structure  in  Figure  16  and  by  the 
eddy  diffusion  approach.  They  are  actually  dynamic  processes  resulting 
from  the  interactions  and  variations  of  solar  radiation,  strong  winds, 
inflows,  power  operations,  etc.  Examples  showing  the  extent  of  mixing 
resulting  from  several  hours  of  strong  winds  are  given  in  Figure  20. 
Entrainment  of  nutrients  from  the  metalimnion  into  the  euphotic  zone 
can  therefore  be  considered  as  an  instantaneous  process  rather  than  a 
gradual  continuous  process. 

In  addition  to  the  seasonal  cycle  in  temperature,  there  are  also 
diel  and  synoptic  cycles  in  temperature  (Figures  21  and  22) .  These 
result  from  variations  in  solar  radiation  and  mixing.  These  temporal 
variations  can  be  Important  in  the  verification  of  mathematical  models 
and  in  data  Interpretation.  Since  most  mathematical  models  use  a 
dally  time  step,  predictions  represent  a  full  day  of  heating  and 
cooling.  In  contrast,  field  measurements  are  usually  taken  during 
daylight  hours  and  frequently  in  mid-afternoon  when  diurnal  heating 
is  maximum.  Diel  variations  can  be  as  large  as  7°C  or  more,  but 
are  typically  the  order  of  1°  or  2°C.  The  magnitude  of  the  variation 
will  depend  on  the  size  and  depth  of  the  lake,  among  other  factors. 

Horizontal  variations  also  occur.  These  can  result  from  differ¬ 
ential  heating,  inflow,  or  mixing.  Differential  heating  occurs  when 
the  smaller  volume  of  water  in  the  littoral  zones  and  headwaters  of 
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Impoundment  warms  or  cools  more  rapidly  than  the  open  water  regions. 

In  large  lakes,  this  phenomenon  Is  significant  and  results  In  the 
formation  of  thermal  bars.  Similarly,  rivers  flowing  Into  a  reservoir 
may  be  of  different  temperatures  creating  longitudinal  variations. 
Examples  of  horizontal  variations  are  shown  In  Figure  23.  Horizontal 
variations  typically  create  temperaure  differences  of  1°  or  20C  or  more. 

Consideration  must  be  given  to  the  above  mentioned  variations 
and  to  the  physical  processes  responsible  for  them  In  data  Interpre¬ 
tation.  If  discrepancies  In  data  are  not  explainable,  then  the  data 
could  be  erroneous.  An  example  of  erroneous  data  Is  shown  In  Figures 
18  and  24  for  Cave  Run  Lake.  The  hypolimnetic  temperatures  for  Julian 
day  213  are  unrealistic  with  respect  to  days  184  and  219.  Examination 
of  the  meteorological  data  for  this  period  confirmed  that  there  were  no 
extreme  events  that  could  have  possibly  caused  such  a  temperature 
deviation. 


DISSOLVED  OXYGEN  DATA 


Two  basic  reasons  for  measuring  dissolved  oxygen  in  lakes  and 
reservoirs  are  aesthetics  and  aquatic  life.  DO  is  required  in  reser¬ 
voirs  to  prevent  the  onset  of  septic  conditions  and  the  accompanying 
malodorous  emissions.  The  development  of  carbon  dioxide,  hydrogen 
sulfide,  and  methane  in  the  sediments  under  anaerobic  conditions  can 
loosen  sludge  which  can  float  to  the  water  surface  and  form  an  unsightly 
scum.  DO  is  also  essential  to  the  metabolism  of  aquatic  life.  EPA  (1976) 
recommends  a  minimum  of  5  mg/f,  DO  to  maintain  good  fish  populations. 

DO  is  an  important  water-quality  parameter  that  Influences  many  other 
chemical  and  biological  processes  in  the  reservoir. 

The  sources  of  DO  to  a  reservoir  are  photosynthesis,  atmospheric 
reaeratlon,  and  hydromechanical.  DO  sinks  are  community  respiration 
and  oxidation  of  organic  matter.  The  magnitude  of  these  sources  and 
sinks  will  depend  on  reservoir  morphology,  water  temperature,  meteorology, 
the  amount  and  characteristics  of  the  oxldlzable  substance,  etc.  For 
lakes  and  reservoirs,  the  most  important  sources  and  sinks  are  probably 
photosynthesis,  community  respiration,  and  autochthonous  organic  matter 
(see  Table  I  and  Ruttner,  1963). 

Before  discussing  DO  distributions  in  reservoirs,  the  solubility 
properties  of  DO  in  water  need  to  be  reviewed.  DO  Is  poorly  soluble 
In  water,  because  It  does  not  react  chemically  with  water  and  because 
air  contains  only  about  21  percent  oxygen  by  volume.  The  solubility 
of  DO  Is  proportional  to  its  partial  pressure.  It  therefore,  varies 
with  pressure  and  temperature  (Figure  25).  Potential  problems  with 
low  DO  can  occur  at  low  pressure  (i.e.,  high  altitudes)  and  at  high 
temperatures.  Other  parameters  such  as  salinity  can  reduce  the  solu¬ 
bility  of  DO. 
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TABLE  I 


Average  DO  Budget  In  Lake  Keystone  on  CLOUDLESS  DAYS 
(after  Eley,  1970) 


X  Gains  on  Clear  Days  (metric  tons  O2) 


Production 

+ 

976.1  day  ^ 

Diffusion  in  (X  k=1.49) 

+ 

42.0  day~^ 

Cimarron  R.  Inflow 

+ 

6.8  day  ^ 

+ 

1024.9  .  .  . 

+1024.9  day 


X  Losses  on  Clear  Days  (metric  tons  O2) 


Respiration 

-  787.8  day~^ 

-1 

Diffusion  out  (X  k  =  1.49) 

-  194.8  day  ^ 

-1 

Reservoir  discharge 

1.4  day 

-  984.0  . -  984.0  day 

X  Calculated  Net  Gain  +  40.9  day  ^ 
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DO  Distributions  -  Ollgotrophlc  lakes  are  typically  characterized 
by  low  nutrient  concentrations  and  low  production  of  organic  matter. 

The  DO  distribution  is  controlled  by  physical  processes  such  as  tempera¬ 
ture,  pressure,  and  stratification  instead  of  biological  processes. 

The  resulting  DO  profile  is  called  an  orthograde  profile.  In  an  idealized 
orthograde  profile,  DO  is  100  percent  saturated  and  its  distribution 
will  depend  only  on  pressure  and  temperature  (Figure  26).  Profiles 
of  this  type  are  very  seldom  found  in  nature.  Usually  the  DO  concentrations 
in  an  oligotrophic  lake  are  constant  throughout  the  water  column. 

In  contrast  to  oligotrophic  lakes,  eutrophic  lakes  are  character¬ 
ized  by  high  nutrient  concentrations  and  high  organic  production. 

The  DO  in  the  hypolimnion  is  depleted  by  the  oxidation  of  the  settling 
organic  matter  and  by  oxidation  at  the  sediment-water  Interface.  The 
density  gradients  in  the  metallmnion  prevent  renewal  of  DO  by  mixing 
and  light  is  usually  not  sufficient  for  renewal  by  photosynthesis. 

The  resulting  DO  profile  is  called  cllnograde  and  it  generally  follows 
the  shape  of  the  temperature  profile  (Figure  27).  All  things  being 
equal,  insight  into  the  productivity  of  a  lake  can  be  obtained  by 
examining  a  DO  profile  to  determine  if  it  is  orthograde  or  clinograde. 

Power  operations  can  also  affect  DO  distributions.  DO  might  be 
expected  to  Improve  with  the  Increased  mixing  from  power  operations, 
but  this  is  not  necessarily  true.  In  Figure  19,  power  operations  did 
increase  the  temperatures  in  the  hypolimnion,  but  DO  did  not  improve. 

This  example  reaffirms  the  point  that  estimates  of  mixing  from  tempera¬ 
ture  profiles  alone  are  not  sufficient  to  answer  all  water-quality 
questions. 

DO  Profile  Variations  -  Not  all  DO  profiles  are  purely  orthograde 
or  clinograde.  Some  exhibit  maxima  or  minima  in  the  metallmnion. 

The  maxima  are  due  to  photosynthesis  by  algae  and  aquatic  plants. 

The  high  stability  of  the  metallmnion  prevents  the  DO  from  diffusing. 

DO  profiles  containing  a  maxima  are  called  positive  heteorograde 
curves. 

Negative  heteorograde  curves  denote  metalimnetic  DO  minima.  The 
minima  can  result  from  the  accumulation  of  oxldlzable  material,  basin 
morphology,  zooplankton,  and  the  temperature  gradient  alone.  The 
sources  of  the  oxldlzable  material  accumulating  in  the  metallmnion  can 
be  settling  from  above  (autochthonous)  or  density  currents  (allochthonous) . 
In  reservoirs,  both  are  possible.  The  exact  source  can  only  be  deter¬ 
mined  by  field  studies.  Attempts  at  such  studies  are  described  by 
Gordon  and  Skelton  (1977)  and  Drury  and  Gearheart  (1975) .  It  is  possible 
that  allochthonous  demands  will  dominate  in  the  headwaters  and 
autochthonous  demands  will  dominate  near  the  dam.  The  basin  morphology 
is  important  if  there  are  large  flats  of  sediment  in  the  metallmnion. 
Oxidation  at  the  sediment-water  interface  will  create  the  minima. 

Temperature  is  Important  because  the  solubility  of  DO  varies  with 
temperature.  Ruttner  (1963)  argued  that  the  temperature  gradient 
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alone  can  account  for  the  metallmnetlc  oxygen  minima  because  oxidation 
will  precede  much  faster  in  the  epllimnion  and  metallmnlon  than  in 
the  hypolimnion  due  to  the  warmer  temperatures.  In  addition,  DO  will 
Increase  in  the  epllimnion  from  transfer  through  the  air-water  Interface, 
but  not  in  the  metalimnion  and  hypolimnion  because  mixing  is  inhibited. 

The  net  result  is  a  metallmnetlc  minima. 

Wetzel  (1975)  stated  metallmnetlc  maxima  are  much  more  prevalent 
in  lakes  than  metallmnetlc  minima.  Yet,  reservoir  data  (e.g..  Figures 
19  and  26),  continually  show  metallmnetlc  minima.  There  are  two  possible 
reasons  for  this  difference.  First,  in  a  reservoir,  sufficient  quantities 
of  allochthonous  material  may  be  entering  the  metalimnion  via  density 
currents  causing  the  minima.  Second,  in  natural  lakes  there  appears 
to  be  a  relationship  between  the  depth  of  the  thermocline  and  light 
penetration.  Photosynthesis  can  still  occur  in  the  metalimnia  causing 
DO  maxima.  However,  in  a  reservoir  this  may  not  be  the  case.  The  in¬ 
creased  turbidity  from  inflows  and  the  increased  mixing  from  advection  will 
result  in  a  deeper  epllimnion  and  may  keep  light  from  penetrating 
into  the  metalimnion.  The  result  would  be  a  DO  minima  Instead  of  a 
maxima. 

Seasonal  Variations  -  The  seasonal  variation  in  DO  as  measured  in 
DeGray  Lake  for  three  consecutive  years  is  shown  in  Figure  28.  During 
periods  of  turnover,  DO  was  distributed  uniformly  throughout  the  water 
column.  Both  metallmnetlc  minima  and  maxima  occurred.  The  minima 
were  prevalent  during  the  summer  for  the  three  years  while  the  maxima 
were  not  so  prevalent  and  occurred  only  during  the  spring.  The  general 
increase  in  DO  from  1974  to  1976  illustrates  the  improvement  in  water 
quality  that  can  be  expected  as  a  new  impoundment  progresses  through 
its  initial  transitory  period. 

Temporal  and  Horizontal  DO  Variations  -  DO,  as  temperature,  exhibits 
diel  and  synoptic  variations  (Figures  21  and  29) .  Both  examples  are 
related  to  solar  radiation  (i.e.,  clear  weather)  and  the  resultant 
biological  activity  (photosynthesis  and  respiration) .  If  temperature 
was  the  dominant  factor,  then  the  diel  DO  variation  in  Figure  21  would 
be  inversely  proportional  to  temperature  as  implied  in  Figure  25. 

According  to  Wetzel  (1975),  the  oxygen  regime  in  the  littoral 
zone  is  totally  different  from  the  pelagic  (open  water)  zone.  Two 
reasons  for  this  difference  are  the  large  number  of  aquatic  plants  and 
the  smaller  volume  of  water  found  in  the  littoral  zone.  For  these  two 
reasons,  the  littoral  zone  is  usually  characterized  by  a  more  pronounced 
diel  cycle,  which  in  turn,  creates  horizontal  DO  variations  between 
the  littoral  and  pelagic  regions  of  a  lake.  Horizontal  DO  variations 
are  also  created  at  the  end  of  the  growing  season  when  large  popu¬ 
lations  of  aquatic  plants  die  and  decompose.  Since  reservoirs  charac¬ 
teristically  have  large  shoreline  development  ratios  (i.e.,  extensive 
littoral  zones),  the  oxygen  regime  in  the  littoral  zone  can  play  a 
significant  role  in  the  total  oxygen  budget  of  a  reservoir. 
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Allochthonous  inflows  to  a  reservoir  can  also  act  as  a  source  or 
sink  of  DO  and  thereby  create  horizontal  variations.  The  oxygen 
demand  will  be  greater  in  the  headwater  regions  resulting  in  lower  DO. 


OXYGEN  DEMANDS 


Three  general  classes  of  matter  that  exhibit  oxygen  demands  in 
reservoirs  are  carbonaceous  material,  oxidizable  nitrogen,  and  chemical 
reducing  compounds.  The  impact  of  each  of  these  on  the  DO  distribution 
will  depend  on  many  things  including  the  amount  of  oxidizable  material, 
the  settling  rate  of  the  oxidizable  material,  the  water  temperature, 
and  the  bottom  profile  and  depth  of  the  lake.  Methods  to  measure  these 
oxygen  demands  in  reservoirs  include  biochemical  oxygen  demand  (BOD) , 
chemical  oxygen  demand  (COD),  and  total  organic  carbon  (TOC),  among 
others. 

BOD  -  BOD  is  defined  as  the  amount  of  DO  required  by  bacteria 
to  stabilize  decomposable  organic  matter  under  aerobic  conditions. 

It  has  been  widely  used  to  determine  the  pollution  strength  of  domestic 
and  industrial  wastes. 

The  BOD  test  is  a  procedure  that  requires  dilute  wastes,  no 
toxic  substances,  accessory  nutrients,  and  a  seed  consisting  of  a 
group  of  diverse  organisms.  The  reaction  is  assumed  to  be  complete 
in  20  days  although  it  theoretically  goes  to  infinity.  It  is  usually 
measured,  however,  in  5  days  when  only  a  portion  of  the  total  BOD  is 
exerted.  The  precise  portion  depends  on  the  character  of  the  seed  and 
the  nature  of  the  organic  matter  but  it  is  typically  70  to  80  percent 
of  the  total  BOD.  The  total  or  ultimate  BOD  is  estimated  from  the 
5-day  BOD  assuming  a  first  order  reaction  with  a  constant  rate  coefficient 
For  more  specific  information  on  the  BOD  test,  see  Standard  Methods 
and  Sawyer  and  McCarty  (1967). 

An  idealized  BOD  reaction  is  shown  in  Figure  30.  The  initial 
demand  is  due  to  carbonaceous  material.  After  about  8  days  significant 
nitrification  begins.  The  5-day  BOD  was  selected  to  avoid  the  oxygen 
demand  from  nitrification.  The  nitrification  demand  can  be  measured 
by  measuring  the  various  forms  of  nitrogen  and  using  stoichiometry 
to  calculate  the  demand. 

Since  water  samples  from  rivers,  lakes,  and  reservoirs  contain 
significant  populations  of  nitrifying  organisms,  nitrification  may 
begin  early  and  thereby  affect  the  interpretation  of  the  5-day  BOD. 
Although  the  action  of  nitrifying  bacteria  can  be  arrested  by  inhibiting 
agents  (e.g.,  methylene  blue)  or  by  reducing  populations  through 
pasteurization,  chlorination,  or  acid  treatment,  there  is  not  a  standard 
procedure  for  measuring  the  carbonaceous  BOD  in  natural  waters  (Sawyer 
and  McCarty,  1967).  In  addition,  any  type  of  algae  growth  in  such 
waters  makes  BOD  data  difficult  to  interpret. 
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In  interpreting  BOD  data,  the  assumptions  of  the  BOD  test  need  to 
be  considered.  The  first  assumption  is  the  reaction  is  first  order. 

Both  nitrification  and  protozoa  that  eat  bacteria  can  cause  deviations 
from  a  first  order  reaction.  Second,  it  is  assumed  that  the  rate 
constant  is  truly  a  constant.  Table  II  shows  that  this  is  not  always 
true.  The  constant  will  depend  on  the  nature  of  the  organic  material 
and  the  ability  of  the  organisms  to  utilize  organic  matter.  Since 
these  two  assumptions  are  not  always  met,  extrapolating  5-day  BOD  to 
total  BOD  may  be  subject  to  error. 

COD  -  The  COD  is  the  oxygen  required  for  oxidation  of  organic 
matter  to  carbon  dioxide  (CO2)  and  water.  It  is  based  on  the  fact 
that  most  organic  compounds  can  be  oxidized  by  the  action  of  strong 
oxidizing  agents  under  acid  conditions.  A  complete  description  of 
the  COD  test  is  found  in  Sawyer  and  McCarty  (1967)  and  Standard  Methods. 

The  disadvantages  of  the  COD  are: 

(1)  Organic  matter  is  converted  to  CO2  and  water  regardless 
of  the  biological  assimilability  of  the  substance.  The  COD  is  usually 
greater  than  the  BOD  (Table  II). 

(2)  The  procedure  can  reduce  some  Inorganic  ions  causing 
erroneously  high  results. 

(3)  There  is  no  evidence  of  the  rate  at  which  biological 
active  material  will  stabilize. 

The  advantages  of  the  COD  are  that  it  requires  a  short  time 
for  analysis  and  it  is  not  subject  to  toxic  conditions. 

TOC  -  The  amount  of  carbonaceous  material  in  a  water  sample  can 
be  obtained  from  a  total  carbon  analyzer  (i.e.,  infrared  spectro¬ 
photometer).  Through  appropriate  procedures  (see  Standard  Jlethods) ,  the 
organic  carbon  can  be  separated  from  the  Inorganic  carbon  to  obtain 
TOC.  It  is  also  possible  to  separate  TOC  into  dissolved  organic 
carbon  (DOC)  and  particulate  organic  carbon  (POC)  by  filtration, 
settling,  or  centrifugation. 

The  advantage  of  TOC  measurements  is  that  only  carbonaceous 
material  is  considered.  It  is  also  possible  to  relate  TOC  and  POC 
to  plankton  and  zooplankton  (Figure  31) . 

The  disadvantages  of  using  TOC  measurements  for  estimating 
oxygen  demand  are  that  not  all  of  the  carbon  may  be  biologically 
oxldizable  and  that  there  is  no  evidence  on  the  rate  of  oxidation. 

Discussion  -  Typical  BOD  and  COD  temporal  and  spatial  variations 
are  given  in  Table  II  for  Lake  Keystone.  The  values  decrease  from 
the  headwaters  (Station  I)  to  the  dam  (Station  IV)  indicating  the 
decay  of  allochthonous  material.  In  addition,  the  COD  is  always 
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TABLE  II 


SPATIAL  AND  TEMPORAL  VARIATION  OF  MEAN  CONCENTRATIONS  OF 
BOD  AND  COD  AND  MEAN  BOD  REACTION  CONSTANTS 

IN  LAKE  KEYSTONE  DURING  1966-67  (after  Eley,  1970) 


5-day 

BOD 

20-day 

BOD 

BOD 

k 

COD 

Station  I 

3.8 

11.2 

0.035 

176 

Station  II 

2.5 

6.6 

0.041 

46 

Station  III 

1.4 

5.2 

0.029 

28 

Station  IV 

1.1 

4.5 

0.026 

18 

8/1/66 

1.4 

9/24/66 

2.1 

4.7 

33 

11/19/66 

3.1 

7.9 

0.042 

42 

12/17/66 

2.8 

8.0 

0.038 

48 

1/28/67 

3.6 

11.6 

0.031 

163 

3/29/67 

2.8 

9.2 

0.028 

103 

6/4/67 

1.6 

7.5 

0.019 

36 

7/23/67 

1.2 

4.7 

0.026 

22 

8/24/67 

1,1 

3.1 

0.042 

29 
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greater  than  the  BOD  values  because  It  Is  a  measure  of  all  oxldlzable 
material,  not  just  the  biologically  oxldizable  material. 


In  the  section  on  DO,  it  was  stated  that  the  DO  budget  in  lakes 
and  reservoirs  was  dominated  by  photosynthesis  and  community  respiration. 
If  this  statement  is  correct,  then  these  processes  must  be  considered 
when  evaluating  the  oxygen  demand  in  reservoirs.  Since  the  BOD,  COD, 
and  TOC  do  not  include  the  effects  of  photosynthesis  and  community 
respiration,  their  value  for  estimating  the  oxygen  demands  in  reservoirs 
is  questionable.  It  should  also  be  remembered  that  the  BOD  and  COD 
tests  are  primarily  Intended  to  measure  the  pollution  strength  of 
wastes,  not  the  oxygen  demands  in  the  reservoirs. 

There  is  not  single  method  available  to  estimate  the  oxygen 
demand  in  a  reservoir.  A  combination  of  BOD  tests  and  TOC  measure¬ 
ments  provides  useful  Information.  The  BOD  will  give  a  measure  of 
the  biologically  oxldizable  organic  matter,  provided  nitrification  is 
arrested.  TOC  measurements  estimate  the  total  quantity  of  dissolved 
and  particulate  organic  matter.  Together  they  can  be  used  to  estimate 
biologically  resistant  organic  matter.  COD  measurements  are  totally 
worthless  in  studying  reservoirs  because  they  contribute  nothing  to 
the  understanding  of  oxygen  dynamics. 

Perhaps  the  best  procedure  for  evaluating  oxygen  budgets  and 
demands  in  reservoirs  is  a  time  series  of  DO  measurements.  All 
naturally  occurring  phenomena  are  therefore  taken  into  consideration 
in  the  proper  proportions.  Measurements  of  this  type  are  sometimes 
referred  to  as  oxygen  deficits. 


SUMMARY 


The  subject  of  temperature  and  DO  data  interpret  ition  was  discussed 
from  the  point  of  view  that  all  data  are  explainable.  This  approach 
requires  a  complete  understanding  of  the  numerous  interactions  between 
the  morphometry,  physics,  chemistry,  and  biology  of  the  system.  It  is 
only  with  this  complete  understanding  that  a  reservoir  can  be  managed 
effectively  with  respect  to  water  quality. 
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Figure  1.  Relationship  between  surface  area  and  volume  for  187  CE  impoundments. 
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Figure  2.  Relationship  between  surface  area  and  drainage  area  for  187  CE  impoundments. 
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Figure  3.  Relationship  between  total  annual  flow  and  drainage  area  for 
187  CE  Impoundments. 


Thermocline  depth  as  a  function  of  surface  area.  (after  Ford,  1976) 


BLANTON  (1973) 

MORTIMER  (HUTCHINSON,  1957) 
OTHERS 


S/jll  ‘AilAISnjJIfl  Aa03 


Paper  7 


SURFACE  AREA,  KM* 

Figure  5.  Increase  in  hypolimnetic  eddy  dif fusivities  with  lake  surface  area 
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Figure  6.  Response  of  surface  temperatures  and  mixed  layer  depths 
In  McCarrons  Lake,  Minnesota,  to  dally  changes  In  solar 
radiation  and  wind  speed. 
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Figure  8.  Basic  terminology  for  thermal 
stratification. 
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Figure  11.  Relationship  between  the  percent  light  absorbed  in  the  top  two  feet  and  the 
extinction  coefficient  (G.  Drummond,  ORD,  unpublished  data). 
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Figure  12.  Relationship  between  the  extinction  coefficient  and  the 
Secchl  disk  depth  (G.  Drummond,  ORD,  unpublished  data). 
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Typical  seasonal  variations  in  Secchl  disk  depths  (after 
Culver,  1975). 
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Figure  15.  Spatial  and  temporal  variations  In  back 
calculated  eddy  dlffuslvltles,  McCarrons 
Lake,  Minnesota. 
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Temporal  variations  In  Secchl  disk  depths 
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DISSOLVED  OXYGEN ,  m^/jL 

Figure  19.  Effect  of  power  operations  on  temperature  and  dissolved 
oxygen  profiles. 
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Figure  20.  Examples  of  wind  mixing. 
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Figure  22.  Comparison  of  daily  (synoptic) 
and  weekly  isotherms. 
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CAVE  RUN, 1975 

Figure  24.  Example  of  erroneous  temperature  profile  for 
day  213. 
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DISSOLVED  OXYGEN  PROFILES,  8/23/74 

Figure  26.  Comparison  of  measured  and  saturated  DO  profiles. 
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Lgure  27.  Example  of  cllnograde  oxygen  profile. 


Figure  30.  Classical  BOD  reaction  where  L  is 
the  20-day  BOD  minus  nitrification. 
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COLIFORMS  AND  WATER  QUALITY: 

USE  OF  DATA  IN  PROJECT  DESIGN  AND  OPERATION 


By 

Kent  W.  Thornton* 


INTRODUCTION 


Conform  bacteria  are  Important  in  many  CE  Impoundments  because 
of  their  Impact  on  project  purposes  such  as  recreation  and  water  supply. 
Knowledge  of  coliform  transport  and  die-off  may  provide  information 
for  project  design  and  operation  that  minimizes  potential  problems. 

The  selection  of  appropriate  recreation  sites  can  be  improved  by  knowing 
the  path  of  coliform  transport  into  and  through  an  impoundment.  This 
same  information  may  be  useful  in  determining  withdrawal  depths  for 
water  supply. 

The  coliform  bacteria  are  defined  to  Include  all  the  aerobic  and 
facultative  anaerobic,  gram-negative,  nonsporulatlng  bacilli  that 
produce  acid  and  gas  from  the  fermentation  of  lactose  (APHA,  1971) . 

These  organisms  are  primarily,  but  not  exclusively,  of  human  enteric  or 
excretal  orgin.  Two  of  the  most  common  members  of  the  coliform  group 
are  Escherichia  coll  and  Enterobacter  aerogenes . 

Another  group  of  bacteria  of  interest  in  a  water  quality  sampling 
program  are  the  fecal  streptococci  bacilli.  Fecal  streptococci  are 
also  found  in  the  gastrointestinal  tract  of  warm  blooded  animals, 
including  man,  but  are  typically  nonpathogenlc .  Several  species  of 
fecal  streptococci  and  their  hosts  include  Streptococcus  faecalis  -  man, 
bovis  -  cattle,  and  equinus  -  horses.  Fecal  streptococci  are  also 
found  in  hogs  and  fowl  such  as  chickens,  ducks,  turkeys,  etc. 

The  primary  purpose  of  determining  the  coliform  bacteria  concen¬ 
tration  in  water  is  to  indicate  the  possible  presence  of  water-borne 
pathogenic  organisms.  Pathogenic  organisms  that  cause  typhoid, 
dysentery,  and  infectious  hepatitis  are  also  of  human  origin  and  may 
be  transmitted  through  water.  However,  even  during  epidemics,  these 
pathogenic  organisms  may  be  present  in  low  concentrations,  reducing 
the  possibility  of  detection.  In  addition,  culturing  the  organisms  for 
Identification  is  tedious  and  time  consuming.  For  these  reasons,  a 
surrogate  indicator  such  as  fecal  collforms  is  desirable.  To  be  an 
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effective  indicator,  a  microorganism  should  possess  several  attributes: 
(1)  it  should  be  applicable  to  all  types  of  waters  such  as  lakes, 
streams,  reservoirs,  estuaries,  etc.;  (2)  It  should  be  present  whenever 
pathogens  are  present  with  a  survival  time  comparable  to  that  of  the 
hardiest  pathogen;  and  (3)  It  should  not  reproduce  In  contaminated 
waters  thereby  resulting  in  Inflated  values  (Scarplno,  1974).  Unfortu¬ 
nately,  no  indicator  organism  Is  known  that  satisfies  all  these  criteria. 
Fecal  conform  bacteria  do,  however,  meet  some  of  the  requirements  -  they 
are  of  enteric  origin,  have  similar  die-off  times  as  some  pathogens,  and 
have  been  correlated  with  an  Increase  or  decrease  in  some  pathogens 
such  as  Salmonella  typhl  (typhoid) .  In  addition,  the  collform  bacteria 
are  relatively  easy  to  culture. 

This  paper  emphasizes  the  water  quality  relations  of  collform 
distributions  rather  than  public  health  aspects.  While  the  two  are 
Inexorably  related,  the  discussion  centers  on  the  interpretation  of 
collform  data  with  respect  to  water  quality.  Sampling  and  culturing 
techniques  and  procedures  are  outlined  In  AFHA  Standard  Methods  (1975) 
and  Mlllipore  Application  Manual  AM302  (1973) • 


SOURCE  IDENTIFICATION 


While  fecal  collform  (FC)  and  fecal  streptococci  (FS)  bacteria 
are  present  in  both  man  and  other  livestock,  their  abundance  varies  from 
species  to  species  (Table  1) .  The  preponderance  of  bacteria  In  man 
consists  of  FCs.  The  reverse  is  true  for  livestock;  that  is,  FS 
are  the  most  abundant.  This  Information  can  be  used  to  assist  In  the 
Identification  of  possible  origins  of  waste  Input. 

Since  one  group  of  organisms  originate  predominantly  from  man  and 
the  other  group  from  livestock,  a  ratio  of  the  two  may  provide  infor¬ 
mation  on  the  probable  source  of  the  waste  discharge  (Table  2) .  A  FC 
to  FS  ratio  greater  than  4.0  would  Indicate  human  origin  since,  as 
Indicated  In  Table  1,  man  is  the  only  animal  with  a  FC/FS  ratio  greater 
than  4.0.  Similarly,  the  various  species  of  livestock  have  FC/FS  ratios 
less  than  0.7.  A  range  of  FC/FS  ratios  from  2. 0-4.0  or  0. 7-1.0 
Indicates  origins  primarily  from  human  or  livestock  origin,  respectively. 
Ratios  between  1.0  and  2.0  are  Inconclusive  and  are  not  really  indica¬ 
tive  of  either  source.  In  this  Instance,  the  sampling  station  should 
probably  be  moved  upstream  closer  to  the  source  of  the  waste  discharge. 

The  above  ratios  were  computed  based  on  three  underlying  limi¬ 
tations:  (1)  the  samples  were  collected  no  more  than  24  hr  travel  time 
downstream  from  the  source  of  the  waste;  (2)  the  pH  of  the  water  was 
within  the  range  of  4. 0-9.0;  and  (3)  only  FCs  are  compared  with  FS 
(Mlllipore,  1973).  Differential  die-off  of  the  two  groups  after  24  hr 
obfuscates  the  conclusions  that  can  be  derived  from  the  ratio  values. 

This  differential  die-off  minimizes  the  usefullness  of  the  ratio. 

Extreme  pH  values  also  affect  the  survival  of  the  bacterial  groups, 
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again,  resulting  In  differential  die-off  rates.  Finally,  consideration 
of  total  conforms  rather  than  FCs  may  result  In  erroneous  conclusions 
since  the  proportion  of  total  conforms  originating  from  human  sources 
Is  highly  variable.  The  total  collform  group  Includes  Enterobacter 
aerogenes  that  are  ubiquitous  In  soil  and  may  be  washed  off  during  storm 
events. 

One  use  of  collform  data,  then,  would  be  to  determine  the  possible 
origin  of  the  waste  Input  and  assess  the  seasonality  of  violations 
of  the  standard  for  body  contact  recreation  (log  mean  of  200  FC  cells/ 
100  ml;  EPA,  1976).  Increased  violations  during  the  summer  months,  due 
to  an  Influx  of  people  may  require  additional  package  treatment  plants 
at  the  project  site  or  Increased  monitoring  of  Inflows. 


STORM  INFLOWS 


Total  and  fecal  collform  concentrations  generally  Increase  In  the 
Inflow  to  a  project  during  storm  events.  This  Increase  Is  a  function  of 
many  factors  such  as  the  size  of  the  drainage  area;  magnitude  of  the 
storm  event;  land  use  In  the  basin;  duration.  Intensity,  and  location 
of  the  rainfall;  and  other  variables.  In  general,  collforms  exhibit 
the  "first-flush"  phenomena.  This  Implies  the  loading  of  collforms 
to  a  stream  occurs  during  the  early  stages  of  a  storm.  Perrier  et  al. 
(1977)  studied  the  storm  event  loading  of  collforms  and  other  water 
quality  variables  to  the  Caddo  River  above  DeGray  Reservoir  In  south 
central  Arkansas.  Although  there  was  a  lot  of  variation  In  the  data, 
the  first-flush  phenomena  was  evident  (Figure  1). 

Regression  equations  were  derived  to  predict  collform  loadings  to 
the  Caddo  River  as  a  function  of  streamflow  during  storm  events  (Perrier 
et  al.,  1977).  Using  data  from  the  entire  hydrograph  to  predict  collform 
concentrations  resulted  in  an  R^  value  of  0.19  or  flow  alone  accounted 
for  only  19  percent  of  the  variance  In  the  FC  concentrations.  Using 
data  from  only  the  rising  limb  of  the  hydrograph  in  the  equation  resulted 
In  an  r2  value  of  0.A2  or  flow  alone  now  accounted  for  A2  percent  of 
the  variance  In  FC  concentrations.  The  regression  equation  using  data 
from  the  rising  limb  of  the  hydrograph 

was 


log  C  «  -0.00263  +  0.8177  log  Q  S.E.  =  0.6736 

where 

C  =  FC  concentration,  cells/100  ml 

3 

Q  *  streamflow,  m  /s 
S.E.  >  standard  error 
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By  monitoring  several  storm  events  a  year  for  several  years,  it  may  be 
possible  to  quantify  the  relation  between  coliform  concentrations  and 
streamflow.  This  knowledge,  coupled  with  the  following  discussion, 
may  permit  reservoir  operations  personnel  to  predict  when  coliform 
violations  may  occur  and  for  how  long  they  may  persist.  In  addition, 
planners  may  be  able  to  determine  the  most  appropriate  sites  for 
recreation  areas. 

Given  the  magnitude  of  the  streamflow  and  its  associated  coliform 
load,  analyses  can  be  performed  to  determine  if  coliform  concentrations 
exceed  standards  and  how  far  into  the  pool  these  violations  may  persist. 


IN-POOL  WATER  QUALITY 


Storm  events  have  been  monitored  by  the  USAE  Waterways  Experiment 
Station  (WES)  in  both  the  Caddo  River  and  as  the  events  proceeded  through 
DeGray  Reservoir  (Figures  2-4).  These  storm  events  were  monitored 
using  a  transmissometer  to  measure  the  change  in  the  percent  transmission. 
A  percent  transmission  of  100  percent  would  indicate  pure  water,  for 
example,  whereas  0-percent  transmission  Indicates  very  turbid  water 
that  permits  no  transmission  of  light  through  the  light  path.  By 
using  spatial  and  temporal  variations  in  percent  transmission,  it  was 
possible  to  trace  the  progress  of  stormflow  into  and  through  the  reser¬ 
voir.  Knowledge  of  where  the  flow  goes  also  permits  an  assessment  of 
the  fate  of  quality  constituents. 

The  storm  flow  that  occurred  on  24  October  1976  entered  as  an 
underflow  (Figure  5).  As  it  proceeded  into  the  reservoir,  it  appeared 
to  become  an  interflow  near  Station  12  on  27  October  1976  (Figure  6). 

The  inflowing  streamwater  had  an  initial  temperature  of  ca.  15°C. 

A  comparable  density  strata  existed  at  around  11  or  12  m  below  the 
surface.  Comparison  with  the  temperature  profile  for  this  month 
(Figure  7)  indicated  15°C  temperature  ca.  12  m  below  the  surface. 

On  3  March,  the  inflowing  stream  waters  were  approximately  10°C 
while  the  reservoir  was  nearly  Isothermal  at  approximately  8  or  9°C 
(Figure  8).  The  streamflow  initially  followed  the  old  stream  channel 
as  it  entered  the  reservoir  and  stayed  near  the  bottom  CFigure  9) . 

However,  as  the  stormflow  proceeded  further  into  the  reservoir,  the 
buoyant  forces  of  the  slightly  less  dense  water  started  to  dominate 
and  the  flow  rose  in  the  pool  until  it  moved  primarily  as  an  overflow. 

On  16  June,  the  inflowing  stream  water  was  approximately  19-20°C.  It 
is  evident  from  Figure  10  that  the  storm  water  proceeded  into  the  pool 
as  a  rather  narrow  interflow  at  approximately  8m.  A  review  of  Figure  11 
will  indicate  a  strong  thermal  stratification  with  the  20°C  water 
located  at  approximately  8  m  below  the  surface. 

This  information  is  quite  useful  in  interpreting  other  water 
quality  data.  A  perusal  of  data  from  Red  Rock  Reservoir  (Figure  12) 
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will  indicate  a  reasonable  correspondence  of  suspended  solids  concen¬ 
trations  and  FC  concentrations.  Coliforms  may  be  washed  into  the 
stream  unattached  or  adsorbed  to  solid  particles.  There  is  also 
evidence  that  some  coliforms  may  enter  the  water  column  through  bed 
scour  (Grimes,  1975;  Matson  et  al.,  1978;  Van  Donsel  and  Geldreich, 

1971).  The  use  of  percent  transmission,  then,  may  be  an  excellent  means 
of  tracking  the  initial  movement  of  coliforms  into  the  pool.  If  this  is 
true,  then  there  should  be  some  correspondence  between  peak  coliform 
concentrations  in  the  pool  and  the  expected  movement  of  stormflows 
through  the  reservoir. 

A  comparison  of  coliform  profiles  at  various  sites  in  the  reser¬ 
voir  following  storm  events  demonstrated  the  correlation  among  coliforms 
and  storm  events.  The  storm  event  of  24  October  1976  entered  initially 
as  an  underflow  and  then  proceeded  as  an  interflow  at  approximately 
the  10-m  depth  at  Station  12.  This  same  pattern  is  illustrated  in  the 
FC  concentrations  in  the  pool  (Figure  13).  Samples  were  taken  at 
1,  3,  5,  10  m  and  5-m  Increments  thereafter,  so  it  is  not  possible  to 
state  the  FC  concentrations  at  intermediate  depths  that  correspond 
to  the  depth  of  the  maxima  flow  regime.  By  27  October  1976,  the 
Inflow  had  proceeded  to  Station  12  in  the  reservoir.  The  peak  con¬ 
centrations  were  located  around  10  m.  Station  13  had  a  peak  at  5  m, 
which  corresponded  roughly  to  the  peak  increase  in  turbidity  on 
25  October  1976.  By  3  November  1976,  the  coliform  concentration 
at  both  stations  had  diminished  below  the  standard  of  200  cells/100  ml 
due  to  both  die-off  and  settling. 

The  storm  event  of  3  March  1977  provided  additional  insight  into 
coliform  transport.  On  3  March  the  streamflow  plunged  into  the  pool 
following  the  thalweg.  As  it  proceeded  down  the  reservoir,  the  slightly 
warmer  stream  water  began  to  rise  toward  the  surface.  On  3  March,  the 
coliform  concentration  at  the  upper  Station  13  was  greatest  near  the 
bottom  (Figure  14).  On  5  March,  the  storm  inflow  had  proceeded  to 
Station  10  with  the  leading  edge  located  at  a  depth  of  approximately  5  m. 
It  is  evident  from  Figure  14  that  this  same  pattern  occurred  with 
respect  to  coliforms.  The  streamflow  and  associated  FC  concentration 
were  distributed  throughout  more  of  the  pool  at  Station  12.  This  was 
a  function  both  of  the  settling  of  suspended  solids  and  coliforms  as 
well  as  the  stormflow  rising  in  the  pool  due  to  the  buoyant  density 
differences. 

The  settling  process  was  demonstrated  even  more  vividly  by  the 
storm  event  of  16  June  1977.  The  storm  inflow  was  narrowly  defined 
as  an  interflow  at  a  depth  of  approximately  7-8  m  on  18  June.  The 
coliform  concentration  peaked  sharply  at  5  m  on  18  June  and  diminished 
rapidly  on  either  side  at  Station  12  (Figure  15).  By  21  June,  the  inflow 
had  proceeded  to  Station  8  (Figure  15).  The  FC  concentrations  were  now 
distributed  more  evenly  between  5  and  10  m  at  Station  12.  With  these 
flow  patterns  and  coliform  profiles,  it  is  possible  to  compute  infor¬ 
mation  that  is  useful  in  locating  recreational  sites  and  in  determining 
the  possible  Impact  of  storm  events  on  existing  recreational  sites  or 
water  supply. 
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APPLICATION  OF  STORM  EVENT  DATA 


The  information  from  the  movement  of  the  turbidity  plumes  can  be 
used  to  compute  estimates,  albeit  rough,  of  the  speeds  of  various 
storm  events  as  they  proceed  through  the  reservoir  (Table  3) .  Since 
the  distance  among  stations  is  known  and  the  measurement  times  were 
recorded,  it  is  possible  to  compute  the  speed  in  terms  of  miles 
per  hour  (m/s  '  0.5  x  mph) .  In  addition,  recurrence  intervals  on  a 
seasonal  basis  can  be  computed  from  long-term  flow  records  to  determine 
the  frequency  of  occurrence  of  a  given  storm.  In  this  example, 
based  on  climatological  similarities,  winter  was  considered  December- 
February,  spring  represented  March-May,  and  fall  was  represented  by 
September-November .  The  seasonal  basis  is  important  since  erroneous 
conclusions  can  be  reached  for  operational  purposes  if  an  annual 
basis  is  used.  The  peak  discharge  for  the  3  March  1977  storm  was 
approximately  22,000  cfs  while  the  peak  discharge  for  the  16  June  1977 
storm  was  ca.  13,000  cfs.  Although  the  June  storm  had  a  lower  peak 
discharge,  on  a  seasonal  basis,  this  was  equivalent  to  a  1  in  100  yr 
event  while  the  March  storm  was  equivalent  to  a  1  in  44  yr  event. 

If  enough  storm  events  are  measured,  it  may  be  possible  to  roughly 
correlate  the  speed  or  change  in  speed  with  which  a  storm  proceeds 
through  a  reservoir  with  the  Inflow  discharge.  While  actual  measure¬ 
ments  are  obviously  the  most  desirable,  it  also  may  be  possible 
to  compute  some  of  these  speeds  using  infrared  aerial  or  land- 
satellite  photos  of  turbidity  plume  movement  in  a  project,  especially 
if  the  event  is  an  overflow. 

An  assumption  generally  made  for  coliform  die-off  is  that  it  is 
a  first  order  process;  that  is,  die-off  is  exponential.  Gunnison  (un¬ 
published  data)  conducted  a  series  of  experiments  to  predict  the  effect 
of  temperature  on  coliform  die-off  using  natural  lake  water  as  the  media 
rather  than  an  artificial  media.  Gunnison  (personal  communication) 
obtained  significantly  greater  die-off  of  conforms  using  natural 
water  than  artificial  media.  The  results  of  his  study  for  three 
temperatures  are  shown  in  Table  4.  Figure  16  was  constructed  for 
illustrative  purposes  using  an  initial  concentration  of  100,000  organisms 
and  the  equations  of  Gunnison.  At  5°C,  it  takes  approximately  15  days 
for  the  FC  concentration  to  decrease  from  100,000  cells/100  ml  to 
200  cells/100  ml.  This  time  is  decreased  to  ca.  4  days  at  20°C  with 
only  a  slight  decrease  in  time  (3.5  days)  at  35°C.  The  three  die-off 
rates  may  be  plotted  versus  temperature  and  estimates  of  die-off 
rates  at  other  temperatures  Interpolated  (Figure  17).  While  the 
scientific  validity  of  this  interpolation  is  questionable,  the  values 
obtained  appear  to  be  reasonable.  The  die-off  rate  at  20°C  is  approxi¬ 
mately  double  that  at  10°C,  which  could  be  expected  from  a  Q^q  assump¬ 
tion  of  a  doubling  of  a  chemical  reaction  rate  for  each  10°C  increase 
in  temperature.  Whether  biological  reaction  rates  obey  the  assump¬ 
tion  is  a  point  of  contention,  but  it  may  be  possible,  then,  to 
obtain  estimates  for  various  coliform  die-off  rates  as  a  function  of 
temperature  under  controlled  conditions. 


6 


Paper  8 


The  storm  event  Information  on  collform  concentrations  collected 
in  the  field  has  already  been  presented.  The  peak  FC  concentrations 
have  been  tabulated  as  a  function  of  reservoir  station  and  storm 
event  (Table  5).  With  several  assumptions,  this  laboratory  and  field 
information  can  be  used  to  assist  in  project  operation  and  design. 


If  we  assume  there  is  an  exponential  die-off  of  conforms  as  a 
function  of  temperature,  then 


C  = 


where 


C  =  conform  concentration  at  time  t  and  some 
specified  temperature 

C  =  initial  conform  concentration 
o 

X  =  die-off  rate  as  a  function  of  temperature 
t  =  time  in  days. 

Then, 


t  = 


X  =  - 


From  the  laboratory  data,  Cq  and  x  are  known,  t  is  an  independent 
variable  so  C  may  be  computed.  From  the  field  data,  C  and  Cq  are 
known,  and  using  t  as  the  travel  time,  x  may  be  computed,  or  x  may  be 
assumed  to  be  identical  to  the  laboratory  rate  at  a  given  temperature 
and  t  computed. 

Using  this  approach,  it  is  then  possible  to  compare  field  rates 
and  times  with  laboratory  rates  and  times.  For  example,  if  we  ^'sumed 
plug  flow  (Lagranglan  reference  frame)  with  no  dilution  or  entrainment, 
etc.,  based  on  a  given  temperature  and  time  of  travel,  we  could 
compute  a  die-off  rate.  Conversely,  if  we  compute  or  assume  a  die-off 
rate  at  a  given  temperature,  we  can  compute  the  expected  time  required 
to  reach  the  concentration  of  C  cells/100  ml.  Again,  this  assumes 
no  dilution  or  losses  due  to  settling.  A  comparison  of  these  values  is 
shown  in  Table  6. 
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For  the  storm  event  of  24  October  1976,  there  was  excellent  agreement 
between  the  two  approaches  indicating  that,  for  this  particular  event, 
plug  flow  may  have  been  a  representative  model.  In  general,  however, 
the  "actual"  die-off  and  the  time  required  were  much  greater  and 
shorter,  respectively,  than  predicted  under  laboratory  conditions. 

This  is  expected  due  to  the  effects  of  dilution  and  settling  in 
addition  to  die-off.  The  laboratory  equations  would,  in  general, 
therefore,  represent  a  conservative  case  for  collform  die-off  and  travel 
times. 

This  Information  can  be  used  to  assist  in  the  location  of  recrea¬ 
tional  sites  and  in  anticipating  possible  body  contact  standard 
violations  (Table  7).  By  installing  or  using  recording  gages  and  record¬ 
ing  thermographs  on  major  tributaries  into  CE  projects,  both  flow  and 
temperature  can  be  continuously  monitored.  Thermal  model  predictions 
and  actual  temperature  profiles  in  the  reservoir  over  several  years 
provide  good  information  on  what  the  average  temperature  profile  will 
be  at  a  given  time.  Information  from  equations  of  flow  versus  FC 
concentrations  provides  estimates  of  the  initial  inflowing  collform 
concentrations.  Using,  first,  the  laboratory  equations  and  appropriate 
rates,  the  time  required  to  reach  the  200  FC  cells/100  ml  can  be  computed. 
Once  this  time  is  obtained,  information  from  a  table  such  as  Table  3 
giving  average  speeds  as  a  function  of  storm  events  can  be  used 
to  compute  how  many  miles  into  the  reservoir  the  inflow  will  proceed 
before  the  collform  concentration  decreases  below  the  body  contact 
standard.  The  seasonal  frequency  of  the  storm  event  can  be  used  to 
assess  the  significance  of  this  degree  of  penetration.  For  example, 
if  a  storm  event  that  proceeds  6  miles  into  the  project  has  a  recurrence 
interval  of  1  in  100  yr,  the  decision  to  locate  sites  further  up  the 
Impoundment  for  greater  recreational  benefits  may  be  reasonable. 

This  same  approach  may  be  used  to  determine  when  potential  problems 
or  violations  may  occur  at  an  existing  recreational  site.  This  infor¬ 
mation  may  help  determine  where  and  when  additional  sampling  and 
monitoring  of  existing  sites  is  desirable  to  minimize  health  risks. 


DISCUSSION  AND  LIMITATIONS 


The  decrease  in  collform  concentrations  within  the  water  column 
is,  as  suggested  by  Zanonl  et  al.  (1978)  perhaps,  best  described  by 
a  disappearance  rate  rather  than  die-off  rate.  Studies  by  Matson  et  al. 
(1978)  indicate  the  sediments  may  serve  as  a  reservoir  for  collform 
bacteria.  The  effects  of  settling  and  entrainment  result  in  a  decrease 
in  collform  concentrations  but  do  not,  necessarily.  Increase  the 
actual  die-off  rate. 

The  magnitude  of  the  sampling  error  and  resolution  in  the  field 
data  must  constantly  be  kept  in  mind.  An  example  of  this  magnitude 
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in  sampling  error  may  be  obtained  by  considering  the  average  values  of 
conforms  collected  at  discrete  depths  versus  that  of  composite  samples 
from  DeGray  Reservoir  (Table  8) .  It  is  known  that  bacteria  may  exist 
in  layers  only  a  few  centimeters  thick  and  can  easily  be  missed  by 
sampling  only  at  a  fixed  depth.  The  error  in  one  Instance  is  almost 
three  orders  of  magnitude  in  estimating  concentrations. 

The  site  specificity  of  these  conclusions  is  also  important.  The 
DeGray  Reservoir  has  stable  land  use  with  minor  point-source  discharges. 

A  system  with  changing  land  use  and  multiple,  significant  point  sources 
will  respond  quite  differently.  The  site  specificity  must  be  incorporated 
in  all  sampling  programs  as  well  as  in  any  conclusions  formed  from  the 
data  set. 

All  conclusions  and  approaches  are  predicated  on  obtaining  good 
field  data  from  CE  projects.  The  importance  of  good  field  data  should 
never  be  underestimated.  It  is  extremely  difficult  to  refute  sound 
conclusions  drawn  from  a  well-designed  and  implemented  field  data 
collection  program. 
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TABLE  1 

AVERAGE  INDICATOR* 


10®  CELLS/g  FECES 


RATIO 

ANIMAL 

F  COLI 

F  STREP 

FC/FS 

MAN 

13.0 

3.0 

4.4 

DUCK 

33.0 

54.0 

0.6 

SHEEP 

16.0 

38.0 

0.4 

CHICKEN 

1.3 

3.4 

0.4 

COW 

0.2 

1.3 

0.2 

PIG 

3.3 

84.0 

0.04 

*  GELDREICH  &  KENNER,  1969. 
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TABLE  2 

INTERPRETATION 
FC/FS  RATIO 


FC/FS  >  4.0 
FC/FS  <  0.7 
2  <  FC/FS  <  4.0 
0.7  <  FC/FS  <  1.0 
1.0  <  FC/FS  <  2.0 


HUMAN  WASTE 
LIVESTOCK  WASTE 
PRIMARILY  HUMAN 
PRIMARILY  LIVESTOCK 
?  MOVE  NEARER  SOURCE 
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TABLE  3 

STORM  EVENT  MOVEMENT 


STORM 

EVENT 

aT 

(DAY) 

MILES 

MPH 

iNIT 

STA 

TERM 

STA 

RECURRENCE 
INTERVAL  (YR) 

10/24/76 

1.12 

8.9 

0J3 

16 

13 

1  IN  25 

3/3/77 

0.58 

10.3 

0.74 

16 

13,12 

1  IN  44 

3/5/77 

1.91 

6.1 

0.13 

13,12 

10 

6/17/77 

1.41 

8.9 

0.26 

16 

13 

1  IN  98 

6/18/77 

1.12 

8.6 

0.32 

13 
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TABLE  4 


Y  =  3.1  X 

Y  =  11.4 

Y  =  5.7  X 


EQUATIONS  DESCRIBING 
COLIFORM  DIE-OFF 

10«  EXP(-0.42T)  5“C  R 

10*  EXP(-1.64T)  20^  C  R 

10®  EXP(-1.86T)  35  C  R 


^  =  0.84 
^  =  0.99 
^  =  0.97 
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TABLE  5 

STORM  EVENT  DATA 


STORM 

EVENT 

STATION 

16 

13 

10 

10/24/76 

14,600* 

3600 

3/3/77 

5,700 

1250 

3/5/77 

1250 

590 

6/17/77 

60,000 

560 

6/18/77 

560 

26 

*  F  COLI  CELLS/100  ML. 
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TABLE  6 

COMPARISON  OF  COLiFORM  DIE-OFF  RATES  AND  TIMES 


STORM 

EVENT 

TEMP 

’C 

RATE  (DAV^) 

TIME  (DAY) 

LAGRANGIAN 

ACTUAL 

THEOR 

ACTUAL 

10/24/76 

15 

-1.23 

-1.25 

1.13 

1.12 

3/3/77 

10 

-0.82 

-2.61 

1.85 

0.58 

3/5/77 

10 

-0.82 

-0.39 

0.91 

1.91 

6/17/77 

20 

-1.64 

-3.31 

2.85 

1.41 

6/18/77 

20 

-1.64 

-2.74 

1.87 

1.12 
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TABLE  7 


COMPUTATION 

GIVEN:  STORM  DISCHARGE.  Q 
TEMPERATURE 
INITIAL  CONC 

LOG  C  =  0.00263  +  0.8177  LOG  Q 
TIME  CALCULATION 

,  S(i) 

X 

Y  =  200  CELLS/100  ML 
A  =  C 

X  =  UB  OR  KNOWN  FROM  FIELD  CALC 
DISTANCE 

MILES  =  MPH  X  TIME 
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TABLE  8 
STATION  12 


MEAN  FOR  DEPTHS  COMPOSITE  OF 

0.  3,  5  METERS  UPPER  5  METERS 

DATE  TOTAL  COLIFORMS  TOTAL  COLIFORMS 


1977  CELLS/100  ML  CELLS/100  ML 


BEGINNING  OCTOBER  2 


FIGURE  2a 


FIGUKE  2b 


STORM  EVENT  BEGINNING  OCTOBER  24,  1976 


FIGURE  2c 


SURFACE  OBSERVATION  -  MUD  PLUME 
14:17  HR  MARCH  3,  1977 


FIGURE  3a 


CM 


PRIMARY  SAMPLING 
LOCATION 


STORM  EVENT  BEGINNING  MARCH  3,  MH 


FIGURE  3b 
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STORM  EVENT  BEGINNING  MARCH  3.  1977 
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FIGURE  3d 


STORM  EVENT  BEGINNING  JUNE  16,  19 


SURFACE  -  PERCENT  TIME  OBSERVATIONS 
13:40  HR  JUNE  23.  1977 


DEPTH 


FIGURE  7 


10/9/76 


TEMP,  ®C 
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STORM  EVENT  BEGINNG  MARCH  3,  1977 


ATI  ON  S 


FIGURE  10b 
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FIGURE  11 
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WEEK  OF  SAMPLING 


10/27/76 


FIGURE  13a 


f 

O 

g 

\ 

(/) 

-J 

-J 

LJ 

o 


1^ 

-J 

o 

o 

_1 

< 

o 

UJ 


o 


I 


Paper  8 


11/3/76 


FECAL  COLIF.,  CELLS/100  m£ 


AQUATIC  BIOLOGY  -  DATA  INTERPRETATION 
(PART  I):  PRIMARY  PRODUCER  ORGANISMS 
(ALGAE  AND  MACROPHYTES) 


By 

John  W.  Barko,  Research  Biologist 
Environmental  Effects  Laboratory 
Waterways  Experiment  Station 
Vicksburg,  MS  U.S.A. 


INTRODUCTION 


The  process  of  photosynthesis  utilizes  solar  energy  to  convert 
carbon  into  organic  matter,  and  thereby  satisfies  the  nutritional  and 
energetic  requirements  of  most  organisms  within  the  biosphere.  Organisms 
capable  of  photosynthesis  are  referred  to  as  primary  producers.  These 
organisms  are  autotrophic  -  meaning  that  they  require  only  a  source  of 
kinetic  energy  (light)  and  certain  essential  nutrients  for  self- 
maintenance.  Through  their  metabolism,  primary  producers  considerably 
Influence  the  quality  of  the  environment. 

In  aquatic  ecosystems,  the  primary  producers  are  the  algae  and 
macrophytes  ("large  plants").  Because  of  their  significant  interactions 
with  other  components  cf  the  aquatic  ecosystem,  these  organisms  markedly 
affect  water  quality.  Therefore,  an  understanding  of  the  role  of 
aquatic  primary  producers  is  appropriate  for  the  proper  evaluation  and 
interpretation  of  water  quality  data. 

Interactions  among  aquatic  primary  producer  organisms,  other 
aquatic  organisms,  and  the  physical  and  chemical  aquatic  environment 
are  discussed  herein.  This  article  provides  basic  information  on  these 
interactions  and  elucidates  some  of  the  limitations  and  constraints 
associated  with  related  methods  of  parameterization. 


ECOLOGICAL  CLASSIFICATION 
OF  AQUATIC  PRIMARY  PRODUCERS 


Aquatic  macrophytes  are  simply  large  aquatic  plants.  Although 
the  term  "macrophyte"  is  most  commonly  used  in  reference  to  the  more 
structurally  complex  aquatic  flowering  plants,  certain  large  algae 
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and  also  aquatic  bryophytes  (mosses)  are  also  referred  to  as  macrophytes. 

For  example,  the  macroalga,  Chara ,  and  the  moss,  Fontlnalis  are  macrophytes. 
Aquatic  macrophytes  can  be  classified  Into  four  categories  on  the  basis 
of  their  life  form. 


Life  Form _ Example 


emergent 

cattail 

floating- leaved 

water  Illy 

submergent 

pondweed 

free-floating 

duckweed 

Members  of  the  first  three  categories  (emergent,  floating-leaved, 
and  submergent  plants)  are  rooted  and  attached  to  a  substratum.  Free- 
floating  plants  often  have  roots,  but  these  do  not  function  In  attachment. 
Emergent  aquatic  plants,  although  adapted  to  constant  or  periodic  standing 
water,  are  most  similar  to  terrestrial  plants.  These  have  easily 
distinguishable  and  functionally  separable  aerial  and  rooted  portions. 
Submergent  plants  are  the  moat  highly  adapted  to  the  aquatic  environment. 
These  plants  remain  totally  submersed  throughout  their  lifespan,  lack 
structural  rigidity,  and  therefore  are  dependent  upon  the  buoyancy  of 
the  water  for  their  support.  Floating-leaved  plants  are  Intermediate 
In  form  between  emergent  and  submergent  life  forms,  and  have  character¬ 
istics  of  both.  Similar  to  the  emergent  aquatic  plants,  floating¬ 
leaved  plants  possess  separate  aerial  and  rooted  portions.  In  addition 
to  the  large  floating  leaves,  smaller  leaves  similar  to  those  of  submerged 
plants  are  often  produced  at  the  sediment  surface. 

Algae  are  structurally  simplified  aquatic  plants.  In  general, 
they  lack  functionally  differentiated  structures  (l.e.,  roots,  leaves, 
etc.)  and  are  relatively  (compared  to  macrophytes)  small  in  size.  This 
criterion  of  size  necessitates  the  classification  of  macroalgae  as 
macrophytes.  Because  of  the  obvious  ambiguities  Involved,  the  term 
microphyte,  meaning  small  aquatic  plant.  Is  occasionally  substituted 
for  the  term,  algae. 

Algae  can  be  ecologically  classified  on  the  basis  of  their  mode  of 
occurrence.  In  this  respect,  there  are  two  major  categories:  phyto¬ 
plankton  and  attached  algae.  Phytoplankton  are  free-floating,  and 
Inhabit  the  pelagic  region  (l.e.,  open  water)  of  bodies  of  water. 

Attached  algae,  often  called  periphyton,  are  physically  associated  with 
a  substratum.  Their  attachment  can  be  simply  mechanical,  or  It  may 
Involve  a  form  of  mutualism.  For  example,  there  is  some  evidence 
indicating  that  epiphytic  algae  (l.e.,  algae  attached  to  other  plants) 
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may  derive  some  of  their  nutrition  from  the  host  plant.  In  contrast, 
the  attachment  of  algae  to  rock  surfaces  is  usually  simply  mechanical. 
Attached  algae  are  easily  dislodged  by  water  movements;  thus,  they  are 
frequently  found  in  the  phytoplankton,  but  normally  so  in  relatively 
low  abundance.  Also,  some  algal  species  exist  abundantly  In  both  the 
planktonic  and  the  attached  modes  of  occurrence.  Attached  algae  are 
frequently  further  classified  according  to  the  type  of  substratum 
with  which  they  are  associated.  Commonly  used  examples  of  this  type 
of  classification  are  given  below. 


!  of  Attachment 

Substratum 

eplpellc 

fine  sediments 

epilithic 

rocks 

eplsammlc 

sand  grains 

epiphytic 

plants 

eplzolc 

animals 

INTERACTIONS  BETWEEN  AQUATIC  PRIMARY 
PRODUCERS  AND  THE  PHYSICAL  ENVIRONMENT 


Sunlight  is  used  by  primary  producers  in  the  process  of  photo¬ 
synthesis,  presented  below  In  a  simplified  manner. 

sunlight 

and 

pigments 

Inorganic  +  Water - »  Organic  +  Oxygen 

Carbon  Matter 

In  photosynthesis,  sunlight  Is  absorbed  by  various  pigments.  Among  the 
many  pigments  used  In  this  process,  chlorophyll  a  is  the  most  Important 
because  It  Is  directly  Involved  In  the  conversion  of  solar  energy  to 
chemical  energy  as  organic  matter  is  synthesized  from  Inorganic  carbon. 
Only  chlorophyll  a  Is  found  In  all  primary  producers.  Other  pigments 
Include  carotenoids,  phycobollns,  and  other  chlorophylls.  These  other 
pigments  are  termed  "accessory  pigments"  because  they  transfer  absorbed 
light  energy  to  chlorophyll  a.  On  a  dry  weight  basis,  the  mass  of 
chlorophyll  ^  In  plant  cells  Is  usually  much  greater  than  the  mass  of 
other  pigments.  Changes  In  the  relative  proportions  of  different 
pigments  (i.e.,  chromatic  adaptation)  occur  In  response  to  changes 
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in  the  quality  of  light  as  explained  below.  Thus,  both  the  color  and 
light  absorption  capability  of  a  single  species  can  differ  among 
different  habitats. 

Light  of  wavelengths  between  approximately  390  and  760  nannometers 
is  used  in  photosynthesis.  This  range  represents  the  violet  to  red 
portion  of  the  visible  light  spectrum.  In  general,  the  penetration  of 
different  wavelengths  of  light  Into  natural  water  depends  upon  the 
quantity  and  kinds  of  dissolved  materials  and  suspended  particulate 
materials  that  differentially  absorb  and  reflect  portions  of  the  light 
spectrum.  In  many  reservoirs,  light  penetration  is  significantly 
limited  by  turbidity  associated  with  suspended  sediment.  The  same 
effect  can  occur  during  phytoplankton  blooms  that  may  be  dense  enough 
to  also  cause  turbidity.  Turbidity  limits  primary  productivity  by 
decreasing  light  penetration  and  consequently  limiting  photosynthesis. 

Water  temperature  affects  the  rates  of  various  biochemical 
processes.  In  general,  overall  cellular  metabolism  is  temperature- 
dependent.  Different  macrophyte  and  algal  species  exhibit  different 
temperature  optima.  Temperature  and  light  interact  to  define  four 
hypothetical  categories  of  primary  producer  organisms,  each  adapted 
to  a  specific  season  (see  below) . 


low  light 

high  light 

low  temp. 

Winter 

Spring 

high  temp. 

Autumn 

Summer 

Ignoring  for  the  moment  other  important  factors  (l.e.,  chemical, 
biological,  and  other),  temperature  and  light  collectively  account  for 
the  seasonality  of  occurrence  of  primary  producer  organisms. 

Indirect  effects  of  temperature  on  primary  producers  are  complex, 
but  nonetheless  important.  Changes  in  the  density  of  water  with  water 
temperature  change  affect  phytoplankton  buoyancy,  and  ultimately  the 
vertical  distribution  of  phytoplankton.  Temperature  also  affects  the 
availability  of  nutrients.  Rates  of  microbial  mineralization  of 
organic  matter  and  subsequent  releases  of  nutrients  are  greater  at 
higher  than  at  lower  temperatures.  However,  the  solubility  of  gases 
including  carbon  dioxide,  an  important  substrate  in  photosynthesis, 
is  inversely  related  to  water  temperature. 

The  majority  of  macrophytes  are  rooted,  and  their  growth  and 
distribution  are  markedly  affected  by  sediment  type.  The  role  of 
sediments  in  the  nutrition' of  macrophytes  is  variable,  because  nutrients 
may  be  absorbed  either  from  the  water  by  the  shoots  or  from  the  sediment 
by  the  roots.  For  the  most  part,  regions  of  nutrient  absorption 
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depend  upon  the  relative  availability  of  nutrients  with  respect  to  the 
distribution  of  plant  tissues  having  an  absorptive  capability.  Nutrient 
uptake  by  emergent  macrophytes,  which  have  extensive  root  systems  and 
little  absorptive  surface  exposed  to  the  water,  probably  occurs  almost 
exclusively  from  the  sediment.  In  the  case  of  submergent  macrophytes, 
recent  evidence  has  suggested  that  their  roots  are  functionally  similar 
to  those  of  emergent  life  forms.  Therefore,  because  of  the  usually 
greater  availability  of  nutrients  in  sediments  than  In  water,  sediments 
represent  a  potentially  Important  source  of  nutrition  to  these  macrophytes 
Primarily  for  nutritional  reasons,  macrophyte  communities  achieve 
maximum  growth  on  fine-textured  sediments. 

By  virtue  of  their  ability  to  colonize  fine-grained  materials, 
rooted  macrophytes  promote  sediment  stabilization.  Dense  beds  of 
macrophytes  are  also  capable  of  reducing  turbidity  by  decreasing  the 
velocity  of  water  currents  that  pass  through  them;  this  consequently 
Increases  rates  of  sedimentation.  Macrophytes  significantly  affect 
the  shoreline  regions  of  hard  water  lakes  by  promoting  the  development 
of  marl  (calcium  and  magnesltim  carbonate)  formations  that  gradually 
expand  the  area  of  shallow  regions  at  the  expense  of  the  open  water  zone. 


INTERACTIONS  BETl-fEEN  AQUATIC  PRIMARY 
PRODUCERS  AND  THE  CHEMICAL  ENVIRONMENT 


Various  nutrients  are  essential  to  primary  producers.  Carbon, 
nitrogen,  and  phosphorus  are  among  those  nutrients  that  are  required 
In  large  concentrations,  and  are  referred  to  collectively  as  macro- 
nutrients.  Mlcronutrlents  Include  inorganic  elements  and  organic 
compounds  such  as  metals  and  vitamins,  respectively,  that  are  required 
In  small  concentrations.  Nutrients  required  by  aquatic  primary  producers 
are  basically  the  same  as  those  required  by  their  terrestrial  counterparts 

As  Indicated  earlier,  macrophytes,  depending  upon  their  life  form 
among  other  factors,  can  absorb  nutrients  from  either  the  sediment  or 
the  water.  Algae  lack  specialized  structures  for  nutrient  absorption 
and  therefore  derive  their  nutrition  from  the  water.  However,  attached 
algal  species  are  capable  of  supplementing  their  nutrition  through 
nutrient  absorption  from  the  substratum  to  which  they  are  attached. 

This  mode  of  nutrition  Is  particularly  common  among  eplpellc  and 
epiphytic  forms  that  obtain  some  of  their  nutrients  from  sediments 
and  plant  hosts,  respectively.  Summarily,  the  possible  modes  of 
nutrition  for  aquatic  primary  producers  are  given  below. 
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Primary  Producer 
Macrophyte 
emergent 
floating-leaved 
submergent 
free-floating 
Algae 

phytoplankton 
attached  algae 


Nutrition 


sediment 

sediment  and  water 
sediment  and  water 
water 

water 

water  and  substratum 


Both  the  biomass  and  the  productivity  of  primary  producers  are 
controlled  to  a  large  extent  by  nutrient  supply.  Accordingly,  Investi¬ 
gators  for  many  years  have  sought  to  predict  the  growth  of  primary 
producers  from  measurements  of  essential  nutrients.  In  this  regard, 
nutrient  supply  has  been  assessed  from  measurements  of  nutrient  concen¬ 
trations  In  both  the  environment  and  In  the  tissues  of  primary  producer 
organisms.  Additionally,  nutrient  bloassays  have  been  developed  to 
enable  evaluations  of  nutrient  supply  from  fertilization  experiments. 

All  three  of  these  techniques  are  laden  with  assumptions  and  weaknesses; 
these  are  discussed  briefly  below. 

In  evaluations  of  water  chemistry,  the  investigator  frequently 
assumes  that  there  Is  a  relationship  between  concentrations  of  nutrients 
in  water  and  the  growth  of  primary  producers.  Although  one  can  generally 
expect  greater  growth  In  a  system  where  nutrients  are  abundant,  there 
actually  Is  no  simple  relationship  between  the  concentration  of  nutrients 
present  in  the  system  and  their  availability  to  the  primary  producers. 

The  availability  of  a  particular  nutrient  depends  for  the  most  part 
upon  Its  chemical  form,  but  also  upon  the  ease  with  which  It  can  be 
removed  from  the  medium.  For  example,  most  of  the  phosphorus  In  a 
sediment  Is  unavailable  to  rooted  plants  because  It  Is  present  In  an 
Insoluble  form  as  a  component  of  the  sediment  matrix.  The  question  of 
availability  has  been  circumvented  to  some  extent  by  the  development  of 
chemical  techniques  for  determining  soluble  and  extractable  Ions,  which 
are  often  assumed  to  be  available  to  primary  producers.  Superimposed 
on  the  ambiguity  associated  with  determinations  of  "available  nutrients" 
In  the  environment.  Is  the  question  of  nutrient  supply  rate. 
Theoretically,  the  supply  of  a  particular  nutrient  could  be  adequate 
at  concentrations  even  near  zero,  as  long  as  the  rate  of  Its  replenish¬ 
ment  remains  equal  to  the  rate  of  its  removal  from  the  environment  by 
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the  plant.  Certainly,  different  nutrients  are  recycled  within  the 
environment  at  different  rates.  Similarly,  primary  producers  are 
capable  of  selective  removal  of  nutrients  at  different  rates.  Generally, 
measurable  quantities  of  nutrients  in  the  environment  represent  an 
excess.  For  these  and  other  reasons,  correlations  between  environ¬ 
mental  nutrients  and  the  growth  of  primary  producers  are  extremely 
tenuous  and  need  to  be  interpreted  with  great  caution. 

The  determination  of  nutrient  concentrations  in  the  biomass 
(l.e.,  tissue  concentrations)  of  primary  producers  provides  a  more 
direct  means  of  evaluating  nutrient  supply.  Living  tissues  have  fairly 
fixed  requirements  for  specific  nutrients.  At  concentrations  below 
the  optimal  tissue  concentration,  termed  the  critical  concentration  of 
a  particular  nutrient,  growth  of  the  organism  is  limited  by  that  nutrient. 
Conversely,  at  concentrations  above  the  critical  tissue  concentration 
of  a  particular  nutrient,  growth  of  the  organism  is  not  limited  by 
that  nutrient.  Before  this  technique  can  be  used  to  evaluate  nutrient 
supply,  the  critical  concentrations  of  the  relevant  nutrients  must  be 
known.  Unfortunately,  the  process  of  obtaining  this  information 
requires  extensive  experimentation  under  carefully  controlled  con¬ 
ditions.  From  the  literature  available  on  this  subject,  it  appears 
that  critical  concentrations  vary  considerably  among  different  nutrients, 
plant  species,  plant  parts,  and  seasons  of  the  year.  Because  of  the 
great  difficulties  associated  with  the  separation  of  algal  tissue  from 
other  associated  organic  materials,  tissue  nutrient  analyses  in  most 
cases,  are  not  suitable  for  assessing  algal-nutrient  supply.  In 
contrast,  these  analyses  are  generally  quite  suitable  for  assessing 
macrophyte-nutrient  supply;  their  utility  in  this  regard  should 
greatly  increase  as  more  data  on  macrophyte  critical  tissue-nutrient 
concentrations  become  available. 

Nutrient  bioassays  are  essentially  fertilization  experiments. 

This  technique  is  based  on  the  premise  that  nutrient (s)  limiting 
growth  will  support  growth  when  their  concentration  in  the  environment 
is  Increased.  On  one  extreme,  "whole  system"  fertilization  studies, 
both  in  lakes  and  wetlands,  have  provided  verifiable  and  Important 
information  on  the  eutrophication  process.  On  the  other  extreme, 
short-term  "bottle  bioassays"  have  provided  confusing  and  often  mis¬ 
leading  Information  when  used  to  predict  eutrophication.  Short-term 
nutrient  bloassays  provide  Information  on  potential  growth  limiting 
nutrients  only  during  usually  very  brief  and  discrete  periods  of 
experimentation.  Unfortunately,  these  data  are  commonly  extrapolated 
to  other  portions  of  the  year,  and  in  some  cases  even  to  reservoirs 
still  in  the  planning  stage.  Potentially  very  serious  problems  arise 
from  the  possibility  that  nutrients,  other  than  those  implicated  at 
the  time  of  the  bioassay,  may  actually  exert  greater  control  over  the 
annual  growth  of  the  primary  producers.  Another  serious  problem 
related  to  bottle  bioassays  is  the  restriction  on  the  extent  of 
mineral  recycling.  In  nature,  recycling  occurs  among  components  of 
the  entire  ecosystem  and  not  just  within  the  confines  of  a  small 
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volume  of  water.  Other  problems  with  bottle  bloassays  can  be  associated 
with  the  use  of  "test"  organisms  that  often  are  foreign  to  the  environ¬ 
ment  being  examined.  It  Is  extremely  Important  to  realize  that 
different  primary  producers  have  different  nutrient  requirements,  and 
thus  may  actually  demonstrate  a  growth  response  to  different  nutrients 
under  the  same  conditions. 

There  Is  a  dynamic  equilibrium  between  inorganic  and  organic 
components  of  water.  This  equllibrivim  Is  partially  mediated  by  the 
activities  of  primary  producers.  During  the  production  of  biomass, 
primary  producers  assimilate  inorganic  compounds  and  produce  organic 
matter.  When  this  biomass  (organic  matter)  undergoes  degradation.  It 
Is  reduced  to  Its  Inorganic  components.  Thus,  the  inorganic-organic 
cycle  completes  Itself.  Following  periods  of  active  biomass  production, 
nutrients  become  depleted  from  the  environment  and  Incorporated  In  the 
biomass.  Contrastingly,  following  a  period  of  active  biomass  degradation, 
the  converse  occurs.  These  two  processes,  of  course,  occur  simultaneously. 
In  addition  to  changes  In  concentrations  of  nutrients  and  biomass, 
there  are  also  quite  significant  changes  In  the  concentrations  of 
dissolved  gases  associated  with  these  processes.  Oxygen  supersaturation 
and  carbon  dioxide  depletion  result  from  high  rates  of  production; 
rapid  degradation  results  In  increased  carbon  dioxide  and  diminished 
oxygen  concentrations.  The  processes  of  production  and  degradation 
respectively  promote  basic  and  acidic  conditions  particularly  In  poorly 
buffered  and  unmlxed  environments. 


BIOLOGICAL  INTERACTIONS  WITH 
AND  AMONG  AQUATIC  PRIMARY  PRODUCERS 


One  of  the  key  Interactions  between  algae  and  macrophytes  Is 
their  competition  for  light.  For  example,  submergent  macrophytes  are 
disadvantaged  because  of  light  limitation  in  systems  where  phytoplankton 
blooms  are  dense.  During  blooms,  light  penetration  into  the  water  Is 
severely  curtailed;  consequently,  among  the  primary  producers,  only 
those  organisms  at  or  near  the  water  surface  can  survive.  Owing  to 
the  marked  buoyancy  of  many  species  of  blue-green  algae,  these  algae 
are  capable  of  shading  other  algal  species  and  macrophytes  existing 
at  greater  depths.  In  dense  macrophyte  beds,  self-shading  by  component 
plants  has  been  demonstrated  to  affect  not  only  the  productivity  of 
the  macrophytes,  but  also  the  species  composition  of  the  macrophyte 
community. 

Competition  for  nutrients  Is  considered  to  be  an  Important  factor 
affecting  the  seasonal  distribution  and  species  composition  of  algae, 
particularly  the  phytoplankton.  Among  the  attached  algae,  competition 
Is  also  Important.  For  example,  it  seems  reasonable  to  assume  that 
epiphytic  algae  compete  to  some  extent  with  their  host  for  certain 
nutrients.  Similarly,  macrophytes  may  at  times  compete  with  the 
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phytoplankton  for  nutrients  In  the  water.  The  extent  of  this  compe¬ 
tition  would  of  course,  depend  on  the  degree  to  which  the  macrophytes 
derive  their  nutrients  from  the  water. 

Primary  producers  are  affected  by  various  growth-enhancing  and 
growth-inhibiting  substances.  In  general,  these  substances  are  dissolved 
organic  compounds  resulting  from  plant  excretory  and  secretory  processes. 
Some  types  of  algae  are  unable  to  produce  sufficient  quantities  of 
certain  vitamins  for  their  own  use.  Algae,  primarily  those  with 
"animal-like"  characteristics,  are  often  dependent  upon  other  primary 
producers  for  these  materials.  Other  substances  produced  by  both 
macrophytes  and  algae  have  specific  antibiotic  effects  and  thereby 
can  exert  some  control  over  community  composition.  For  example,  there 
is  limited  evidence  indicating  that  certain  macrophytes  are  capable 
of  producing  antibiotics  having  a  rather  broad  range  of  effectiveness 
against  phytoplankton.  Other  macrophytes  have  been  demonstrated  to 
produce  substances  that  inhibit  the  growth  of  zooplankton,  and  in  some 
cases,  the  larvae  of  mosquitoes  as  well. 

The  most  obvious  interaction  between  primary  producers  and  animals 
is  grazing.  However,  the  quantity  of  primary  producer  tissue  consumed 
in  the  living  state  by  grazing  animals  (herbivores)  is  small,  and 
usually  less  than  10%  is  produced.  Due  to  the  selective  feeding  habits 
of  many  grazers,  the  species  composition  of  phytoplankton  can  be 
affected.  Also,  by  reducing  the  overall  biomass  of  phytoplankton  and 
thereby  the  effects  of  competition  under  crowded  conditions,  grazers 
may  often  enhance  phytoplankton  productivity.  The  foliage  and  certain 
reproductive  structures  of  macrophytes  are  used  as  a  source  of  food 
by  a  variety  of  waterfowl  and  mammals.  Although  very  little  macrophyte 
tissue  is  directly  consumed  by  grazers,  many  grazers  feed  on  the 
epiphytically  associated  bacteria  and  algae.  The  tissue  of  primary 
producers  is  ultimately  consumed  in  the  dead  state  by  detritus-feeding 
organisms.  Organisms  feeding  on  dead  particulate  material  in  combi¬ 
nation  with  microbial  agents  of  mineralization  are  ultimately  responsible 
for  converting  organic  matter  to  its  inorganic  components. 


PRIMARY  PRODUCER 
BIOMASS  AND  DENSITY 


The  term  biomass  refers  to  the  mass  of  living  tissue  of  an 
organism  or  a  composite  of  organisms  at  a  given  time.  Biomass  is 
commonly  presented  in  units  of  mass  per  unit  surface  area  or  volume. 
For  example,  the  biomass  of  a  phytoplankton  community  can  be  presented 
as  g/m^  of  water  when  data  on  a  volume  basis  are  Integrated  with 
respect  to  depth  or  as  mg/f.  when  concentration  per  unit  volume  is 
known.  In  the  strictest  sense,  and  as  implicit  in  the  term  Itself, 
biomass  estimates  should  be  corrected  for  dead  mass.  Obviously,  this 
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Is  often  difficult,  and  In  the  case  of  phytoplankton  for  example,  biomass 
estimates  are  frequently  too  large  because  of  the  Inclusion  of  suspended 
dead  particulate  mass.  With  macrophyte  communities,  dead  mass  Is 
relatively  more  easily  separated  from  living  mass;  thus  biomass  and 
litter  mass  categories  are  routinely  distinguished. 

Biomass  can  be  determined  directly  as  fresh,  oven  dry,  or  organic 
weight.  Fresh  weight  Is  not  recommended,  because  It  Is  too  variable; 
but.  It  Is  often  obtained  when  large  quantities  of  plant  material 
need  to  be  processed.  Oven  dry  weight  Is  obtained  by  drying  materials 
In  an  oven  to  a  constant  weight  at  a  temperature  of  100°-105°  Celsius. 
Biomass  data  are  most  commonly  presented  In  this  form.  Where  the  mass 
of  primary  producers  Is  affected  by  Inorganic  contamination,  the 
determination  of  weight  loss  after  Ignition  at  550°  Celsius  allows  the 
estimation  of  organic  weight  In  a  manner  that  Is  unaffected  by  contami¬ 
nation.  In  general,  the  ash  content  (l.e.,  mineral  portion)  of  organic 
matter  Is  less  than  10%.  When  this  percentage  after  combustion  Is 
found  to  greatly  exceed  10%,  contamination  by  nonvolatile  substances 
should  be  suspected.  Such  a  situation  commonly  arises  In  hard  water 
lakes  where  macrophytes  become  encrusted  with  carbonate  deposits.  Also, 
the  belowground  biomass  of  rooted  macrophytes  Is  notoriously  difficult 
to  wash  free  of  adherent  sediment  particles.  These  problems  can  be 
circumvented  by  routinely  combusting  the  plant  materials,  and  expressing 
the  biomass  data  as  organic  weight  (l.e.,  ash-free  dry  weight). 

Estimates  of  algal  biomass  from  chlorophyll  or  cell  volume  data 
are  Indirect,  and  should  be  considered  only  as  approximations.  However, 
because  of  the  previously  noted  difficulties  associated  with  the  separation 
of  algae  from  detritus,  these  approximations  of  algal  biomass  are 
commonly  used.  A  few  precautions  need  to  be  mentioned.  Calculating 
mass  from  cellular  volume  assumes  constant  Intracellular  density,  but 
In  reality  Intracellular  density  Is  quite  variable.  Also,  calculations 
of  cellular  volume  assume  the  conformity  of  algal  cells  to  simple 
geometric  shapes.  These  shapes  need  to  be  carefully  matched  to  different 
algal  species.  In  regard  to  estimates  of  algal  biomass  from  chlorophyll  £ 
data.  It  Is  Important  to  determine  whether  or  not  data  have  been 
corrected  for  chlorophyll  degradation  products  that  Interfere  with 
actual  chlorophyll  a  measurements.  Uncorrected  chlorophyll  £  data 
are  Inappropriate  for  calculating  algal  biomass. 

The  density  of  primary  producers  Is  a  measure  of  biomass  concen¬ 
tric  ion.  More  common  measures  of  density  are  based  on  number  or  volume. 

:n  turn,  these  parameters  are  presented  on  the  basis  of  either  surface 
iroi  ur  volume.  Numerical  density  determinations  are  the  most  common. 

-•'h  phytoplankton,  numbers  of  algal  cells  are  most  frequently  expressed 
i  mimber  per  unit  volume  of  water.  In  the  case  of  attached  algae 
K.*  r  iphytes,  numerical  density  estimates  are  usually  given  on  a 
■  .  -  ir«>a  of  substrate  basis.  The  less  commonly  determined  volumetric 

'  different  algal  species  can  be  useful  In  determining  the 
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relative  Importance  of  each  Individual  species  within  a  phytoplankton 
community. 


PRIMARY  PRODUCTIVITY 


Primary  productivity  is  the  rate  of  production  of  primary  producer 
biomass.  Productivity  is  expressed  as  a  rate  in  units  of  mass  (usually 
dry  weight  or  carbon)  per  unit  area  or  volume,  both  per  unit  time. 
Sometimes  units  of  energy  are  substituted  for  mass,  because  in  the 
process  of  primary  production,  quantities  of  energy  are  fixed  along 
with  the  production  of  mass.  Primary  productivity  is  not  measured  as 
a  single  process.  Indeed,  the  processes  of  gross  primary  production, 
net  primary  production,  and  net  ecosystem  (or  net  community)  production 
are  functionally  distinguishable  from  one  another.  These  processes 
and  their  interactions  with  others  are  depicted  below. 


Solar  Energy 

'i 

Gross  Primary  Production 

^ - >  Autotrophic  (l.e..  Plant)  Respiration 

Net  Primary  Production 

- -  Heterotrophic  (i.e..  Animal  and  Decomposer) 

Respiration 

T 

Net  Ecosystem  Production 

Net  Primary  Production  =  Gross  Primary  Production  -  Autotrophic 

Respiration 

Net  Ecosystem  Production  •  Gross  Primary  Production  -  Ecosystem 

Respiration 

Where:  Ecosystem  Respiration  ~  Autotrophic  Respiration  + 

Heterotrophic  Respiration 


The  diagram  presented  abcve  simplistlcally  demonstrates  the 
relationships  of  energy  and  matter  in  an  ecosystem.  Gross  primary 
production  represents  the  capacity  of  primary  producers  to  utilize 
solar  energy  in  the  fixation  of  inorganic  carbon.  The  organic  matter 
formed  by  this  process  provides  biological  structure  and  a  potential 
energy  base  for  virtually  all  other  biological  processes.  Net  primary 
production  represents  energy  and  matter  produced  in  excess  of  the 
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requirements  of  the  primary  producers  themselves.  These  requirements 
are  measured  as  plant  respiration.  Net  ecosystem  production  represents 
energy  and  matter  produced  in  excess  of  the  requirements  of  the 
ecosystems'  animals  and  decomposers.  Net  ecosystem  production  can  be 
measured  as  newly  accrued  biomass  in  growing  systems,  as  accumulated 
detritus  in  stagnant  systems,  or  as  export  to  adjacent  ecosystems. 

There  are  various  methods  of  measuring  primary  productivity. 

Among  these,  the  biomass,  oxygen,  and  carbon-14  techniques  are  the  most 
common.  Each  of  these  has  inherent  strengths  and  weaknesses,  and 
measures  different  aspects  of  the  process  of  primary  production. 

With  the  biomass  technique,  net  primary  production  is  determined 
as  the  difference  between  minimum  and  maximum  estimates  of  biomass. 

This  method  is  appropriate  only  when  interim  losses  of  biomass  due 
to  grazing,  mortality,  or  other  sources  are  negligible  or  can  other¬ 
wise  be  accounted  for.  For  this  reason,  the  technique  cannot  be  applied 
to  algae.  It  is  best  applied  to  macrophyte  populations  that  demonstrate 
distinct  annual  periods  of  growth.  Determinations  of  the  productivity 
of  perennial  macrophyte  populations  are  often  problematical,  because 
biomass  often  remains  unchanged,  and  also  since  annual  patterns  of 
growth  are  confused  by  variable  rates  of  mortality  and  Internal 
translocations  of  mass  above  and  below  ground. 

The  oxygen  technique  Involves  measurements  of  community  respiration 
in  the  dark  (O2  loss)  and  net  community  production  (O2  gain)  in  the 
light.  The  sum  of  these  two  quantities  provides  an  estimate  of  gross 
primary  production.  This  technique  does  not  measure  net  primary 
production,  because  plant's  own  respiration  is  Inseparable  from  other 
sources  of  respiration  such  as  microbial  decomposition.  This  technique 
is  applicable  to  the  phytoplankton  for  determinations  of  gross  pro¬ 
duction,  net  community  production,  and  pelagic  (l.e.,  open  water) 
respiration.  In  general,  the  oxygen  technique  is  not  applicable  to 
macrophytes,  because  of  their  ability  to  internally  store  and  recycle 
gases  including  oxygen.  Thus  macrophytes  are  capable  of  obscuring 
changes  in  oxygen  concentration  associated  with  internal  processes 
of  production  and  respiration. 

Among  the  various  methods  of  determining  phytoplankton  productivity, 
the  carbon-14  technique  is  becoming  the  most  commonly  used.  It  is 
advantageous  because  of  its  simplicity  and  directness.  Rather  than 
measuring  changes  in  concentrations  of  substrates  or  products  of 
photosynthesis  as  is  done  with  other  techniques,  the  C-14  technique 
directly  measures  carbon  fixation.  This  measure  is  achieved  by  monitoring 
the  photosynthetic  Incorporation  of  radioactively  labeled  inorganic 
carbon  (l.e.,  carbon-14).  Because  some  of  the  fixed  carbon  is  lost 
to  the  plant's  own  respiration,  the  C-14  technique  is  thought  to 
measure  net  primary  production.  However,  because  this  loss  appears  to 
be  minimal  under  certain  conditions,  this  technique  can  also  be  con¬ 
sidered  to  measure  gross  primary  production.  For  practical  purposes, 
investigators  usually  assume  that  the  C-14  technique  measures  something 
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between  gross  and  net  primary  production.  Whatever  Is  the  case,  this 
technique  affords  excellent  replicability  and  Is  extremely  sensitive. 
Carbon-14  has  been  used  with  some  success  to  measure  the  productivity 
of  attached  algae  as  well  as  phytoplankton.  Although  the  technique 
can  also  be  adapted  to  measurements  of  macrophyte  productivity,  such 
applications  are  methodologically  cumbersome,  and  the  results  are 
difficult  to  Interpret. 


CONCLUSIONS 


In  this  article  I  have  briefly  reviewed  and  summarized  the  more 
Important  aspects  of  the  role  of  primary  producer  organisms  In  aquatic 
ecosystems.  Much  of  the  Information  contained  herein  Is  of  quite 
recent  origin,  and  consequently  subject  to  change  as  new  data  become 
available. 

Primary  producer  organisms  are  an  Integral  part  of  the  aquatic 
environment  and  exert  a  significant  Influence  on  the  quality  of  the 
same.  If  Federal  agencies  And  the  public  In  general  are  seriously 
concerned  about  the  degradation  of  water  quality  In  this  country,  then 
the  research  efforts  of  qualified  laboratories  and  Individuals  must 
continue  to  be  supported.  The  Importance  of  sustained  research  related 
to  the  biology  and  role  of  primary  producers  In  aquatic  ecosystems 
cannot  be  overemphasized. 
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WATER  QUALITY  EFFECTS  OF 
MULTI-PURPOSE  RESERVOIRS  AND  NAVIGATION  DAMS 
ON  THE  MONONGAHELA  RIVER 


by 

Michael  Koryak^ 


I.  INTRODUCTION 

This  is  a  case  study  of  a  basin-wide  system  of  Corps  of  Engineers  water 
resource  projects.  The  discussion  emphasizes  functional  aspects  of  the  data  and 
how  it  is  being  utilized  to  establish  design  and  operational  criteria. 

II.  DESCRIPTION  OF  THE  AREA 

The  Monongahela  River  is  formed  by  the  confluence  of  the  West  Fork  and 
Tygart  Rivers  at  Fairmont  West  Virginia  (Plate  1).  From  Fairmont,  the  Mononga¬ 
hela  River  flows  in  a  northerly  direction  for  128.7  miles  before  joining  the  Alle¬ 
gheny  River  at  Pittsburgh,  Pennsylvania  to  form  the  Ohio  River. 

The  river  drains  7,386  square  miles  of  northern  West  Virginia,  western  Mary¬ 
land,  and  southwestern  Pennsylvania. 

Slack  water  navigation  is  maintained  year-round  on  the  entire  128.7  mile 
length  of  the  Monongahela  River  by  a  series  of  Corps  of  Engineers  low-head  navi¬ 
gation  dams.  There  are  nine  navigation  Locks  and  Dams  (L/D)  on  the  mainstem 
of  the  Monongahela  River  while  the  lower  11.2  miles  are  in  the  pool  of  Emsworth 
L/D,  located  at  Ohio  River  mile  6.2. 

There  are  two  large  tributary  storage  reservoirs  in  the  basin,  Tygart  River 
Lake  cuid  Youghiogheny  River  Lake.  Also,  one  additional  reservoir.  Stonewall 
Jackson  Lake,  is  authorized  on  the  West  Fork  River. 

The  navigation  system  links  the  coal  fields  of  northern  West  Virginia  and 
southwestern  Pennsylvania  with  the  basic  metal  industries  of  the  Pittsburgh  region. 
Consequently,  acid  mine  drainage  from  bituminous  coal  mining  activities  is  the 
princi[^  water  quality  problem  in  the  upper  basin,  and  industrial,  thermal,  and 
domestic  pollution  degrade  the  lower  30  miles. 

The  Monongahela  River  is  very  low  yielding  at  base  flow  and  most  water 
quality  problems  are  low-flow  related. 


1  Aquatic  Uiologist,  Hydrology  and  Hydraulics  Branch,  U.S.  Army  Engineer 
District  Pittsburgh. 
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ni.A.  TYGART  RIVER  LAKE 

Tygart  Dam  is  located  in  northern  West  Virginia,  and  controls  1,184  square 
miles  .of  the  Tygart  River  Basin.  At  summer  pool,  the  lake  is  10  miles  lor^  with 
maximum  and  mean  depths  of  134  feet  and  63  feet  respectively.  By  May  of  each 
year  there  is  usually  110,000  acre  feet  of  water  stored,  and  all  but  9,700  acre 
feet  of  this  total  is  available  for  augmentation. 

There  is  no  authorized  water  quality  storage  in  Tygart  River  Lake.  The 
official  project  purposes  are  flood-control  and  low-flow  augmentation  to  maintain 
uninterrupted  navigation  on  the  Monongahela  and  upper  Ohio  Rivers.  We  can 
and  do  however,  utilize  excess  storage  from  Tygart  River  Lake  for  water  quality 
purposes. 

III.B.  AUGMENTATION  FROM  TYGART  DAM  TO  MITIGATE  ACID  MINE  DRAINAGE 

Tygart  River  Lake  provides  an  assured  minimum  flow  of  340  cfs  at  Opekiska 
L/D  near  the  head  of  the  Monongahela  River  mainstem.  We  also  attempt  to  main¬ 
tain  a  2  to  1  flow  ratio  between  the  discharge  from  Tygart  Dam  and  the  flow 
of  the  mine  drainage  degradated  West  Fork  River. 

Tygart  Dam  has  been  operational  since  1939  and  it's  somewhat  difficult 
to  evaluate  the  effectiveness  of  an  operation  that  has  been  going  on  continuou^y 
for  almost  forty  years.  Plate  2  shows  the  reduction  in  summer  acidity  values 
at  McKeesport,  Pennsylvania  that  occurred  when  operations  to  augment  Mononga¬ 
hela  River  flows  were  initiated  in  1939.  These  values  were  measured  at  a  point 
16  miles  above  the  mouth  of  the  Monongahela  River  and  135  miles  downstream 
of  Tygart  Dam. 

Another  factor  to  be  considered  in  the  evaluation  of  the  data  presented 
in  Plate  2  is  the  fact  that  during  the  depression  years,  prior  to  the  construction 
of  Tygart  Dam,  coal  production  in  the  basin  was  very  low.  With  the  beginning 
of  World  War  R,  shortly  after  completion  of  the  project,  coal  production  in  the 
Monongahela  River  Basin  increased  substantially.  Therefore  the  reduction  in 
summer  acidity  values  shown  in  Plate  2  occurred  during  a  period  when  increased 
acidity  would  have  been  expected. 

III.C.  INFLOW  QUALITY 

Bituminous  coal  mining  is  also  an  important  economic  activity  upstream 
of  Tygart  Dam,  and  periodic  intrusions  of  acid  mine  drainage  enter  Tygart  River 
Lake  during  the  summer  and  autumn  months. 

A  pattern  of  flow  related,  depressed  inflow  pH  values  during  the  summer 
of  1973  can  be  seen  in  Plate  3.  The  outflow  pH  is  also  plotted  in  Plate  3  for  the 
same  period.  In  spite  of  the  extreme,  periodic  degradation  of  the  inflow  water 
a  high  quality  discharge  is  consistently  maintained  at  Tygart  Dam. 
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Biological  samples  were  collected  that  verify  this  observation.  Very  low 
macroinvertebrate  productivity  at  the  inflow  station  indicated  stress  from  acid 
mine  drainage.  At  the  outflow  station  a  5,000%  increase  in  the  avei^e  macroin-  « 
vertebrate  dry  weight  per  square  foot  was  observed  (from  1.2  mg/ft^  to  62.2  mg/ft^). 
The  inflow  samples  were  numerically  dominated  by  acid  tolerant  Chironomidae 
(53.8%).  At  the  outflow,  the  percent  composition  of  Chironomidae  dropped  to 
32.3%,  while  the  more  sensitive  Ephemeroptera  increased  from  0  to  66.1%. 

Also  consistent  with  these  observations  there  is  an  excellent  outflow  sports 
fishery  at  Tygart,  and  no  fishery  at  all  at  the  inflow. 

m.D.  INTERNAL  HYDRODYNAMICS  AND  WITHDRAWAL  ELEVATION 

The  depth  of  penetration  and  resulting  mixing  patterns  of  the  acid  inflows 
into  Tygart  River  Lake  are  influenced  by  the  design  and  operation  of  the  dam. 

The  internal  hydrodynamics  of  the  reservoir  in  turn  influence  the  chemistry  and 
biology  of  both  the  impoundment  and  the  outflow. 

Tygart  River  Lake  has  bottom  withdrawal,  and  a  relatively  short  hydraulic 
retention  time.  The  volume  of  water  normally  discharged  between  June  and  Augfust 
is  equivalent  to  1.65  times  the  storage  at  summer  pool.  In  spite  of  only  moderate 
vertical  thermal  gradients,  inflows  can  penetrate  the  impoundment  as  well-defined 
temperature-density  currents.  Interflow  occurs  during  most  of  the  summer  months 
when  acid  loading  problems  are  most  severe. 

A  seasonal  cycle  of  temperature-density  current  variation  can  be  seen  in 
Plate  4.  Conductivity  values  are  used  in  this  plate  to  follow  the  path  of  the  den¬ 
sity  flows  that  originate  from  the  Tygart  River  inflow.  Beginning  in  May,  the 
impoundmnet  is  still  cool,  nearly  isothermic,  and  conductivity  values  are  low  and 
relatively  uniform  throughout  the  reservoir.  During  this  time  of  year,  the  more 
easily  warmed  inflow  can  be  expected  to  overflow  the  cooler  denser  reservoir 
as  seen  on  the  26  June  1973  survey. 

With  summer  stratification,  interflow  conditions  gradually  develop.  Inter¬ 
flow  is  first  clearly  detectable  on  the  basis  of  conductivity  values  in  early  July 
in  the  upper  reservoir.  This  interflow  penetrated  Tygart  River  Lake  in  a  horizontal 
layer  concentrated  at  approximately  30-40  feet  below  the  surface. 

There  was  constant  mixing  and  entrainment,  especially  in  the  shallow  upper 
portion  of  the  impoundment.  But  for  the  most  part  the  current  retained  its  iden¬ 
tity  and  could  still  be  recognized  when  it  reached  the  dam,  10  miles  downstream 
and  roughly  30  days  after  first  entering  the  impoundment. 

The  conductance  models  were  generally  useful  in  discribing  mixing  patterns 
of  less  conservative  acid  mine  drainage  related  parameters.  For  example,  Plate 
5  demonstrates  a  close  correlation  between  pH  and  conductance  distribution  pat¬ 
terns  for  a  15  August  74  survey. 
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Because  of  bottom  withdrawal,  sharp  slugs  of  acid  that  move  through  Tygart 
River  Lake  as  either  overflow  or  interflow  are  will  mixed  and  diluted  before  they 
are  discharged  from  the  Dam.  This  accounts  for  the  relative  consistency  in  the 
quality  of  the  outflow.  In  the  autumn,  the  inflows  are  likely  to  cool  before  the 
reservoir,  and  underflow  conditions  can  develop,  as  was  noted  in  late  September 
1973  (Plate  4).  Underflows  short  circuit  the  reservoir.  With  bottom  withdrawal 
the  retention  time  of  underflows  is  short,  and  mixing  and  dilution  is  not  as  thor¬ 
ough  as  with  overflow  and  interflow. 

Selective  withdrawal  options  are  not  available  at  Tygart  River  Lake,  but 
if  they  were,  a  higher  elevation  discharge  would  probably  only  be  desirable  during 
underflow  periods.  During  these  periods  the  higher  level  withdrawal  would  protect 
discharge  water  quality. 

Considering  the  nature  of  the  problem  at  Tygart.  the  low  level  outlets  of 
the  existing  structure  generally  provide  an  acceptable  means  for  the  partial  miti¬ 
gation  of  acid  slugs  in  the  outflow  without  sacrificing  reservoir  water  quality 
and  the  reservoir  fishery.  Withdrawal  from  a  higher  elevation  during  the  summer 
months  would  not  only  create  downstream  pH  depressions,  but  would  result  in 
concentrated  acid  interflows  penetratir^  the  reservoir  at  a  higher  elevation. 

This  would  have  a  more  negative  impact  on  the  biologically  sensitive  surface 
strata. 

IV.A,  YOUGHIOGHENY  RIVER  LAKE 

Youghiogheny  Dam  is  located  74  miles  above  the  mouth  of  the  Youghiogheny 
River.  The  Youghiogheny  River  makes  its  confluence  with  the  Monongahela  River 
16  miles  upstream  of  Pittsburgh.  Pennsylvania.  Youghiogheny  Dam  controls  a 
434  square  mile  drainage  area.  Total  storage  at  summer  pool  is  154.000  acre- 
feet  and  the  project  has  been  totally  operational  since  1948. 

Unlike  Tygart  Reservoir,  there  is  authorized  water  quality  storage  in  Youghio¬ 
gheny  River  Lake,  and  there  are  no  significant  sources  of  acid  mine  drainage  up¬ 
stream  of  this  project.  Similar  to  Tygar^however,  the  dilution  and  neutralization 
of  downstream  acid  mine  drainage  is  an  important  aspect  of  the  operation  of  Yough¬ 
iogheny  River  Lake. 

IV.B.  AUGMENTATION  FROM  YOUGHIOGHENY  RIVER  LAKE 

Plate  6  shows  the  pH  of  the  Monongahela  River  above  and  below  the  conflu¬ 
ence  of  the  Youghiogheny  River  in  1965.  Except  for  a  very  brief  period,  down¬ 
stream  pH  was  higher  than  at  the  upstream  station,  indicating  the  year  round 
positive  influence  of  the  Youghiogheny  River.  During  the  summer  months,  the 
pH  of  the  station  above  the  confluence  of  the  Youghiogheny  River  dropped  dra¬ 
matically  while  the  pH  of  the  station  downstream  of  the  Youghiogheny  River 
actually  climbed  during  this  same  period.  No  major  sources  of  alkalinity  other 
than  augmented  Youghiogheny  River  flows  could  be  identified  in  the  reach  between 
the  two  stations. 
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As  illustrated  by  the  data  presented  in  Plate  6,  the  Youghiogheny  River 
Lake  water  quality  augmentation  program  has  been  highly  successful.  One  limita¬ 
tion  to  effective  year  round  water  quality  control  however,  and  a  potential  con¬ 
straint  in  any  water  quality  augmentation  operation,  is  also  evident  in  this  plate. 

In  April  1965  the  pH  of  both  Monongahela  River  stations  was  mildly  depressed. 
Similar  spring  pH  depressions  also  occur  regularly  in  both  the  Tygart  and  Youghio- 
gheny  Rivers.  The  timing  of  these  pH  depressions  corresponds  with  the  spring 
filling  cycles  of  the  tributary  storage  reservoirs.  Excess  run-of  must  be  stored 
during  the  spring  for  low-flow  augmentation  during  the  dry  summer  and  autumn 
months  when  severely  acidic  conditions  are  most  likely  to  develop  in  the  down¬ 
stream  target  areas.  During  the  filling  cycles,  and  when  flood  control  operations 
are  made,  however,  the  percentage  flow  contributions  of  uncontrolled  polluted 
tributaries  can  be  quite  high. 

Acidity  values  in  the  acid  mine  drainage  effected  streams  of  the  upper  Ohio 
River  Basin  are  generally  inversely  related  to  flow.  Therefore,  the  filling  cycle 
problems  are  always  relatively  mild  and  of  short  duration  compared  to  probable 
summer  conditions  without  augmentation. 

Thermal  pollution  and  anoxia  also  exhibit  a  strong  inverse  relationship  to 
flow.  Plates  7  and  8  show  the  mitigating  effect  of  the  cool,  well  aerated  Yough- 
iogheny  River  on  thermal  pollution  and  anoxia  respectively  in  the  lower  Monon¬ 
gahela  River  during  a  1975  summer  low-flow  survey.  The  flow  in  the  Monongahela 
River  upstream  of  the  Youghiogheny  confluence  at  this  time  was  700  cfs.  The 
flow  in  the  Youghiogheny  was  1000  cfs,  75%  of  which  was  augmented  flow  from 
Youghiogheny  River  Lake. 

V.  MONONGAHELA  RIVER  NAVIGATION  STRUCTURES 

Navigation  structures  on  the  Monoi^ahela  River  directly  influence  a  number 
of  physical,  chemical,  and  biological  parameters.  Dissolved  oxygen  is  the  most 
significantly  affected  water  quality  parameter.  The  degree  and  nature  of  the 
effect  is  dependent  on  the  design  and  operation  of  the  individual  structures,  local¬ 
ized  biochemical  oxygen  demands,  and  flow  and  temperature  conditions. 

At  high  flows  turbulent  reaeration  can  be  expected  below  all  of  the  struc¬ 
tures.  Problems  can  occur  however  at  low  summer  flows,  and  because  of  some 
design  features,  particularly  downstream  of  Opekiska  L/D. 

Dissolved  oxygen  stratification  can  develop  in  the  larger  navigation  pools 
at  low  summer  flows.  The  dissolved  oxygen  concentration  of  Opekiska  Pool  during 
the  July  1975  survey  (Q  650  cfs)  presented  in  Plate  8  decreased  vertically  from 
a  maximum  of  10.2  mgA  at  the  surface  to  a  minimum  of  0.1  mg/1  near  the  bottom. 
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The  invert  elevation  of  the  outlets  of  Opekiska  Dam  (fixed  sill  elevation 
831  msl),  is  26  feet  below  the  surface  of  Opekiska  Pool  (normal  pool  elevation 
857  msl).  At  this  depth,  during  low  summer  flow  conditions,  the  concentration 
of  dissolved  oxygen  is  low.  The  outlets  of  Opekiska  Dam  are  submerged  in  Hilde¬ 
brand  pool  (normal  pool  elevation  835  msl)  and  the  low-flow  reaeration  potential 
of  such  an  arrangement  is  limited.  As  illustrated  in  Plate  8,  the  Opekiska  dis¬ 
charge  can  depress  dissolved  oxygen  concentration  over  the  entire  length  of  Hilde¬ 
brand  Pool.  Similar  D.O.  depressions  have  also  been  observed  below  other  sub¬ 
merged  sill,  gated  dams  on  the  mainstem  Ohio  River  at  low  flow. 

A  stratified  pattern  of  dissolved  oxygen  distribution  was  also  apparent  in 
Hildebrand  Pool  during  the  July  1975  low-flow  survey.  The  vertical  dissolved 
oxygen  profile  varied  from  a  maximum  of  6.4  mg/1  at  the  surface  to  a  minimum 
of  3.9  mg/1  at  a  depth  of  29  feet.  The  elevation  of  the  "hypolimnetic”  strata  of 
Hildebrand  Pool  appeared  to  be  influenced  by  the  invert  elevation  of  the  outlets 
of  Hildebrand  Dam  (816  msl  or  19  feet  below  the  normal  pool  elevation). 

The  outlets  of  Hildebrand  Dam  are  elevated  above  the  next  receiving  pool, 
Morgantown  Pool,  and  turbulent  reaeration  from  Hildebrand  Dam  was  sufficient 
to  mitigate  the  dissolved  oxygen  depression  that  was  originally  created  by  the 
Opekiska  Dam  discharge.  Reaeration  at  Hildebrand  Dam  may  also  be  increased 
by  bucket  sills  which  flip  the  discharged  water. 

Similar  to  what  was  observed  in  Hildebrand  Pool,  Morgantown  Pool  appeared 
to  trend  towards  low-flow  dissolved  oxygen  stratification  below  the  invert  eleva¬ 
tion  of  the  outlets  of  Morgantown  Dam.  However,  this  pattern  in  Morgantown 
Pool  was  weakly  developed. 

Some  low-flow  reaeration  was  noted  in  L/D  8  Pool  from  the  discharge  of 
Morgantown  Dam.  Morgantown  Dam  also  has  flip  bucket  sills.  Likewise  below 
L/D  8  some  reaeration  was  measureable  at  low-flow.  The  sills  at  L/D  8  are  ele¬ 
vated  12  feet  above  normal  L/D  7  Pool. 

In  general,  with  the  exception  of  Opekiska  L/D,  either  limited  reaeration 
or  no  significant  change  was  observed  below  all  mainstem  navigation  structures 
on  the  Monongahela  River. 

Existing  gated  navigation  projects  are  being  re-examined  to  determine  the 
feasibility  of  increasing  aeration  at  low-flow  by  the  adjustment  of  the  current 
operation  of  the  gates. 

Since  it  has  been  established  that  these  structures  aerate  when  discharging 
high  volumes  of  water,  it  is  likely  that  further  increasing  the  discharge  per  gate 
at  low  flows  by  using  fewer  gates  with  larger  openings  may  increase  low  flow 
turbulence  and  aeration. 
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We  have  already  had  some  success  with  operationally  induced  turbulence 
on  the  mainstem  Ohio  River.  Hydraulic  studies  at  Montgomery  L/D  for  instance, 
demonstrated  that  1  gate  opened  2  feet  provided  li  times  the  amount  of  D.O. 
as  2  gates  opened  1  foot. 

Even  where  it  may  be  possible  to  successfully  operate  these  structures  for 
reaeration,  the  restrictions  imposed  by  stilling  basin  action,  downstream  riprap 
stability,  currents,  channel  scour,  available  tailwater,  etc.,  may  limit  operations 
to  create  turbulence  for  the  sake  of  reaeration  alone. 

The  problem  can  be  approached  more  effectively  in  the  design  of  new  facil¬ 
ities.  Currently  major  portions  of  a  modernization  program  for  the  Monongahela 
River  navigation  system  are  in  the  planning  stage.  The  program  will  most  likely 
consist  of  replacement  structures  for  L/D  7  and  L/D  3,  and  structural  modifica¬ 
tions  at  L/D  2.  L/D  4,  and  L/D  8.  Low-flow  dissolved  oxygen  conditions  in  the 
lower  river  indicate  that  significant  water  quality  benefits  could  be  realized  at 
low  additional  costs  by  designing  and  incorporating  aerating  hydraulic  features 
into  L/D  2  and  L/D  3. 
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MONONGAHELA  RIVER  ABOVE  McKEESPORT.PA. 

AVERAGE  METHYL  ORANGE  ACIDITY  (  JUNE  -  NOVEMBER  ) 


An-AlOl  a02  CORPS  OF  ENGINEERS  HA5HXNGT0N  DC  F/6  13/2 

AD  AlOl  602  OP  *  SEMINAR  ON  MATER  QUALITY  DATA  INTERPRETATION#  — ETC(U) 


KltOMETCRS  ABOVE  0AM 


KILOMETERS  ABOVE  0AM 


INTERNAL  FLOWS  THROUGH  TYGART  RIVER  LAKE 
AS  INDICATED  BY  CONDUCTIVITY  VARIATION 
MAY  TO  OCTOBER  73 

EACH  ISOPLETH  REPRESENTS  A  lOijmhos/cm 
CONDUCTIVITY  INCREMENT. 
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DEPTH  IN  FEET 
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DEPTH  IN  METERS 


EACH  ISOPLETH  REPRESENTS  A  I.O*C  WT.  INCREMENT 
THE  LOCATIONS  OF  INDUSTRIES  WITH  AN  AVERAGE 
WITHDRAWL  OVER  200  M6D  ARE  INDICATED. 
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